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1 .  Contract  Title;  High-Speed  Strained  Quantum-Well  Lasers  and 

Optoelectronic  Devices 

Number;  NOOO 14-96- 1-0303 

PI;  Professor  S.  L.  Chuang 

University  of  Illinois  at  Urbana-Champaign 
Department  of  Electrical  and  Computer  Engineering 
1406  W.  Green  Street,  Urbana,  IL  61801 
Phone:  (217)  333-3359  Fax:  (217)  333-5701 
e-mail:  s-chuang@uiuc.edu 

Program  Manager;  Dr.  Y.  S.  Park,  Code:  ONR  312 

2.  Technical  Objectives; 

We  investigate  research  issues  on  strained  quantum-well  optoelectronic  devices  and 
their  applications  for  navy  needs.  This  work  studies  the  fundamental  physics  in  strained 
quantum-well  devices  such  as  Fabry-Perot  (FP)  lasers,  distributed-feedback  (DFB)  lasers, 
and  vertical  cavity  surface-emitting  lasers  (VCSELs)  through  the  development  of  theory 
and  experiment.  Integrated  optoelectronic  devices  are  also  investigated  for  applications  in 
wavelength-division-multiplexed  (WDM)  fiber-optic  communications  systems  such  as 
wavelength  tunability  and  wavelength  conversion.  In  this  work,  we  investigate  the 
following  areas: 

a.  Comprehensive  characterization  of  fundamental  parameters  of  strained  quantum- 
well  lasers 

b.  Modeling  and  experiment  of  integrated  laser/modulator  devices 

c.  Modeling  and  experiment  of  VCSELs 

d.  Wavelength  conversion  using  four-wave  mixing  and  cross-phase/gain/absorption 
modulation 

e.  Optical  gain  theory  for  GaN-based  blue-green  quantum-well  lasers 


3 .  Approach 

a.  Use  of  our  amplified  spontaneous  emission  (ASE)  measurements  to  study  the 
intrinsic  parameters  of  high-speed  laser  systems  and  to  assess  the  temperature 
sensitivity  of  these  parameters. 

b.  Use  of  our  high-speed  modulation  and  measurement  system  to  experimentally 
assess  the  carrier  dynamics  and  the  temperature  effects  on  the  high  speed  response 
of  quantum-well  structures  in  both  lasers  and  integrated  devices. 

c.  Comparison  of  experimental  data  and  theory  using  a  comprehensive  model  for  the 
emission  spectra  of  integrated  optoelectronic  devices  such  as  three-electrode  curved- 
waveguide  DFB  lasers  with  electroabsorption  modulators,  taking  into  account  the 
longitudinal  dependence  of  the  carrier  and  photon  densities. 

d.  Theoretical  modeling  of  curved  waveguide  DFB  lasers  with  integrated 
electroabsorption  modulators  for  use  as  wavelength-tunable  laser  sources. 

e.  Comparison  of  experiment  and  theory  for  the  high-speed  modulation  response  and 
relative  intensity  noise  in  selectively  oxidized  VCSELs. 

f.  Conduct  experiments  to  measure  the  effects  of  polarization  dependence  and 
multiple  pump  beams  on  the  conversion  efficiency  of  four-wave  mixing  in 
conventional  and  novel  three-electrode  DFB  lasers. 
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g.  Demonstrate  a  novel  wavelength  conversion  experiment  demonstrating  cross¬ 
absorption  modulation  in  an  integrated  electroabsorption  modulator  /  laser. 

4.  Accomplishments 

A.  Comprehensive  characterization  of  strained  quantum-well  lasers 

I.  Fundamental  parameters  of  DFB  lasers.  We  have  extended  our  model  for  the 
amplified  spontaneous  emission  (ASE)  light  from  the  facet  of  a  laser  diode  to  a  distributed 
feedback  structure.  By  comparing  measured  spectra  to  the  theory  with  excellent  agreement 
we  are  able  to  obtain  important  cavity  parameters  such  as  the  coupling  coefficient  and  the 
facet  phases  which  are  otherwise  very  difficult  to  determine.  Also,  this  measurement  and  fit 
enables  us  to  extract  critical  material  parameters  such  as  the  gain  and  refractive  index  and 
the  corresponding  differential  gain,  differential  index  shift,  and  linewidth  enhancement 
factor  which  are  critical  measures  of  laser  modulation  bandwidth  and  chirp  performance. 
[Publications  1  and  9] 

II.  Temperature  dependence  of  InGaAsP  quantum-well  lasers.  We  have  also 
used  the  facet  light  emission  to  study  the  temperature  dependence  of  bulk  InGaAsP 
semiconductor  lasers.  The  side  light  spontaneous  emission  is  measured  and  compared  with 
the  gain  obtained  from  the  facet  ASE  spectrum  through  the  fundamental  relationship 
between  these  two  parameters.  Through  this  process  we  can  calibrate  the  spontaneous 
emission  measurement  and  obtain  the  actual  radiative  recombination  rate.  From 
measurements  at  different  temperatures  we  have  shown  that  the  characteristic  temperature 
for  the  non-radiative  current  density  is  44  K  for  our  lasers  which  causes  the  high- 
temperature  sensitivity  of  these  lasers.  [Publication  6] 

III.  Strained  InGaAlAs  quantum-well  lasers.  We  have  measured  intrinsic 
parameters  for  both  strained  InGaAsP  and  new  strained  and  unstrained  InGaAlAs  material 
systems.  For  the  InGaAlAs  system  we  have  found  that  compressive  strain  gives  a  larger 
differential  gain  and  refractive  index  shifts  than  no  strain  and  slight  tensile  strain.  Also,  the 
linewidth  enhancement  factor  is  greatly  reduced  for  the  compressively  strained  active 
regions  when  compared  to  that  of  the  unstrained  ones.  When  compared  to  an  InGaAsP 
sample  of  similar  compressive  strain,  the  InGaAlAs  samples  have  a  much  larger 
differential  gain  as  well  as  a  much  larger  refractive  index  shift.  InGaAlAs  active  regions 
should  yield  much  higher  modulation  bandwidths  and  should  provide  good  candidates  for 
devices  such  as  modulators  and  wavelength  converters  which  require  a  large  index  change. 
The  linewidth  enhancement  factor  for  the  two  systems  seem  to  be  similar.  [Publications  16 
and  18] 

IV.  High-speed  microwave  modulation  of  quantum-well  lasers  by  electrical 
and  optical  injection.  We  investigated  optical  injection  in  the  absorption  region  of  a 
quantum-well  laser  and  compared  with  direct  electrical  injection.  We  show  experimentally 
and  theoretically  that  differences  between  the  modulation  responses  for  the  two  injection 
methods  come  from  carrier  transport  and  circuit  parasitic  effects.  We  also  obtained  more 
accurate  measurements  of  several  important  dynamic  laser  parameters  by  using  both 
injection  methods.  [Publication  21] 


B.  Modeling  of  integrated  laser/modulator  devices 

I.  Spatial  hole-burning  and  ASE  of  DFB  lasers.  We  have  developed  a  transfer- 
matrix  based  model  for  distributed-feedback  laser  structures.  With  this  model  we  are  able 
to  treat  the  spatial  dependence  of  the  photon  density  and  carrier  density  self-consistently. 
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We  have  also  achieved  excellent  agreement  with  the  measured  lasing  spectrum  as  well  as  a 
measurement  of  the  longitudinal  carrier  density  profile  obtained  from  the  side  light 
spontaneous  emission  from  a  typical  index -coupled  DFB  laser.  [Publication  12] 

II.  Wavelength  tuning  in  novel  integrated  laser/modulator  devices.  Another 
advantage  of  this  technique  is  that  complicated  structures  with  longitudinal  dependence 
such  as  a  chirped  grating  pitch  as  well  as  integrated  devices  can  easily  be  modeled.  We 
have  applied  this  technique  to  explain  the  behavior  of  integrated  DFB  lasers  integrated 
with  electroabsorption  modulators  experimentally  characterized  in  our  lab.  We  have  also 
had  success  explaining  the  tuning  mechanisms  for  a  tunable  DFB  laser  integrated  with  a 
modulator  which  serves  as  a  phase  tuning  section.  We  show  that  the  wavelength-tuning 
range  can  be  enhanced  by  increasing  the  bending  of  the  curved  waveguide.  We  also  show 
that  the  use  of  the  curved  waveguide  can  suppress  detrimental  spatial  hole  burning  effects 
so  higher  power  operation  can  be  achieved.  [Publication  20] 

III.  High-speed  modulation  of  integrated  devices.  We  have  developed  an 
experimental  setup  to  measure  the  frequency  response  of  integrated  laser  and  modulator 
devices  and  have  compared  the  direct  and  external  modulation  on  the  same  laser. 


C.  Modeling  and  experiment  of  VCSELs 

I.  Modal  properties  of  selectively  oxidized  VCSELs.  We  have  developed  a 
very  efficient  model  to  calculate  the  electromagnetic  properties  of  VCSEL  cavities  based  on 
the  beam  propagation  method  and  a  discrete  Bessel  transform.  This  method  compares 
favorably  to  more  computationally  intensive  models  such  as  the  Finite-Difference  Time 
Domain  model.  We  have  used  this  model  to  study  threshold  properties  and  mode 
competition  in  pillar  structures  and  selectively  oxidized  aperture  devices.  By  determining 
the  cavity  quality  factor  and  resonant  wavelength  of  each  mode  we  can  design  optimized 
structures  for  low  threshold  and  single  mode  operation.  Experimentally,  we  have 
developed  a  fiber-tip  near-field  measurement  system  which  allows  us  to  measure  the  near¬ 
field  pattern  of  the  VCSEL  modes.  We  are  currently  using  these  techniques  to  measure  the 
spatial  and  spectral  properties  of  VCSEL  laser  modes  including  beam  size  and  modal 
competition  of  selectively  oxidized  VCSELs.  We  have  also  measured  the  relative  intensity 
noise  (RIN)  spectrum  and  studied  the  effects  of  competition  between  multiple  transverse 
modes.  [Publication  10] 

II.  Bistability  in  selectively-oxidized  VCSELs.  We  have  performed  extensive 
experimental  characterization  of  selectively  oxidized  VCSELs  provided  by  our  research 
collaborators  at  Sandia  National  Laboratoiy.  We  have  observed  a  strong  correlation 
between  electrical  and  optical  characteristics.  The  onset  of  lasing  threshold  for  both  the 
fundamental  mode  and  for  each  higher  order  mode  is  indicated  by  a  dip  in  I-dV/dl  curve. 
We  have  also  observed  new  phenomena  including  both  electrical  and  optical  bistability. 
[Publication  8] 

III.  High-speed  modulation  of  VCSELs.  We  have  developed  an  experiment 
to  study  the  high-speed  modulation  of  VCSELs  in  both  the  small  signal  and  large  signal 
regimes.  Initial  experiments  for  the  small  signal  modulation  give  a  bandwidth  of  about  8 
GHz.  The  spectral  behavior  is  markedly  different  from  that  of  single-mode  edge  emitting 
lasers.  For  the  large  signal  domain,  we  observe  several  higher  order  harmonics  in  the 
modulation  spectrum  whose  power  dies  down  as  the  current  is  increased.  We  have  also 
observed  interesting  effects  such  as  period  doubling  for  some  operating  conditions.  In 
addition,  the  effects  of  lateral  diffusion,  carrier  transport,  and  multiple  transverse  modes  in 
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VCSELs  were  studied  by  modeling  the  high-speed  modulation  response,  showing  good 
agreement  with  experiment. 


D.  Wavelength  conversion  by  four-wave  mixing,  cross-phase  modulation, 
and  cross-absorption  modulation 

I.  Four-wave  mixing  in  a  DFB  laser.  We  have  performed  a  detailed  study  of 
wavelength  conversion  by  four-wave  mixing  in  distributed-feedback  (DFB)  lasers  where 
the  lasing  modes  serve  as  pump  waves  to  facilitate  the  wavelength  conversion. 

In  a  DFB  laser  with  a  single  longitudinal  mode,  we  have  characterized  the  conversion 
efficiency  up  to  detuning  frequencies  of  500  GHz.  For  small  detunings,  multiple  conjugate 
signals  are  observed  and  measured.  The  dependence  of  each  of  these  on  the  probe  and 
pump  powers  were  measured  and  agree  quite  well  with  a  simple  theoretical  analysis.  For 
larger  detunings,  the  laser  cavity  has  a  pronounced  effect  on  the  efficiency  spectrum  which 
correlates  very  well  to  a  simple  analysis  of  the  experimental  ASE  spectrum  of  the  laser. 
[Publications  9  and  11] 

II.  Four-wave  mixing  in  two-lasing-mode  DFB  laser.  We  have  also  measured 
four-wave  mixing  in  DFB  lasers  with  two  simultaneous  lasing  modes.  We  have 
demonstrated  that  the  shifted  conjugate  conversion  efficiency  is  virtually  independent  of 
pump  separation  even  though  the  frequency  shift  is  increased  by  that  amount.  We  have 
achieved  a  conversion  efficiency  of  2%  for  a  frequency  shift  of  900  GHz  even  though  the 
gain  is  clamped  at  threshold  in  these  lasers.  When  compared  to  the  single  pump  case  for  the 
same  overall  frequency  shift,  this  conversion  efficiency  is  almost  1000  times  greater. 

III.  Four- wave  mixing  in  a  novel  tunable  DFB  laser.  We  have  characterized 
four-wave  mixing  in  a  tunable  three-electrode  DFB  laser  with  a  curved  waveguide.  This 
device  gives  a  tuning  range  of  over  4  nm  and  allows  us  to  tune  to  the  desired  output 
wavelength.  Initial  studies  indicate  that  the  current  injection  profile  and  thus  the  spatial 
profile  of  the  pump  wave  can  have  a  large  effect  on  the  observed  conversion  efficiencies. 

IV.  Polarization  dependence  of  four-wave  mixing.  We  measured  the  spectral 
dependence  of  the  conversion  efficiency  for  polarization  resolved  four-wave  mixing  in  a 
single  mode  distributed-feedback  laser.  We  then  used  our  measurements  of  the  refractive 
index  dispersion  to  calculate  the  polarization  dependent  phase  mismatch  for  four-wave 
mixing.  In  addition,  we  measured  the  polarization  resolved  four-wave  mixing  in  a 
distributed-feedback  laser  with  simultaneous  lasing  TE  and  TM  modes.  The  experimental 
observations  agree  very  well  with  a  qualitative  analysis  of  the  theoretical  mechanism  for 
four-wave  mixing. 

V.  Fundamental  parameters  for  cross-gain  and  cross  phase  modulation.  We 

have  performed  measurements  of  the  gain  and  refractive  index  in  the  presence  of  an  optical 
pump  light.  Depending  on  the  location  of  the  pump  (above  or  below  transparency)  we  see 
enhancement  or  compression  of  the  gain  and  index  spectra.  For  the  refractive  index  spectra 
we  have  demonstrated  that  the  induced  index  shift  is  almost  constant  over  a  100  nm 
wavelength  span.  This  would  indicate  that  efficient  two-way  cross-phase  modulation  could 
be  achieved  over  the  entire  gain  bandwidth  of  the  semiconductor  active  region.  [Publication 
13] 


VI.  Novel  wavelength  conversion  by  cross-absorption  modulation.  We 

demonstrated  wavelength  conversion  by  a  novel  method,  cross-abso^tion  modulation,  in 
an  integrated  electroabsorption  modulator  /  laser  (EML)  for  the  first  time.  The  EML  serves 
as  both  the  conversion  medium  and  the  converted  wavelength  light  source,  a  technique 
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which  requires  less  components  than  other  methods  of  wavelength  conversion  using 
semiconductor  optical  amplifiers.  Cross-absorption  modulation  as  a  function  of  input  and 
converted  signal  power  was  measured  to  study  absorption  saturation  effects  in  the 
modulator.  Dependencies  on  input  signal  wavelength  were  also  studied.  The  modulator 
voltage  can  be  chosen  to  optimize  cross-absorption  modulation  at  a  specific  wavelength. 


E.  Optical  gain  theory  for  GaN  based  blue-green  quantum-well  lasers.  We 

have  published  several  fundamental  papers  on  the  band-structure  and  gain  of  wurtzite 
semiconductors.  The  k.p  method  was  used  to  derive  the  proper  Hamiltonian  for  strained 
nitride-based  semiconductors  with  the  wurtzite  crystal  symmetry.  We  have  subsequently 
applied  this  Hamiltonian  to  the  study  of  the  band-structure  of  strained  GaN  and  InGaN 
quantum  wells.  In  addition,  we  have  added  many-body  effects  which  includes  bandgap 
renormalization  and  Coulomb  enhancement  of  the  gain  model  to  our  calculations  for  these 
material  systems.  We  have  found  that  these  effects  significantly  contribute  to  the  blue-shift 
of  the  peak  gain  in  addition  to  its  magnitude.  Our  model  agrees  very  well  with  the  first  set 
of  experimental  data  published  recently  by  Nakamura.  Recently,  we  have  studied  the 
effects  of  the  strong  piezoelectric  fields  in  GaN  based  systems  on  the  electronic  band 
structure  and  optical  gain.  We  have  found  that  at  lower  carrier  densities  these  induced 
fields  dramatically  lowers  the  calculated  gain  and  can  shift  the  peak  gain  position.  The 
impact  of  the  piezoelectric  field  on  the  laser  performance  has  also  been  investigated,  and  we 
find  that  the  laser  threshold  can  be  reduced  for  a  certain  crystal  orientation  direction  such 
that  the  piezoelectric  field  is  reduced  and  the  band  edge  effective  mass  is  smaller  than  the 
conventional  growth  direction  along  the  c  axis.  [Publications  2-5,  14,  15,  17] 


5.  Significance  of  this  work: 

The  fundamental  research  conducted  here,  including  theory  on  gain,  high-speed 
response,  and  noise  in  lasers,  allows  for  a  better  understanding  of  current  device  limitations 
and  will  enable  higher  performance  device  design.  Our  investigation  of  applications  for 
strained  quantum-well  optoelectronics,  such  as  wavelength-tunable  sources  for 
wavelength-division-multiplexing  (WDM)  fiber-optic  communications  systems, 
wavelength  conversion  in  order  to  enable  wavelength  switching  in  dense  WDM  networks, 
and  high-speed  modulation  in  distributed-feedback  (DFB)  lasers  and  vertical -cavity 
surface-emitting  lasers  (VCSELs)  is  essential  for  enabling  the  growth  of  future  optical 
communications  systems. 
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Theory  and  Experiment  of  Ini_2;Gaa;AsyPi_j/ 
and  Ini_2;_j^Gaa;Al2/As  Long-Wavelength 
Strained  Quantum- Well  Lasers 

J.  Minch,  S.  H.  Park,  T.  Keating,  and  S.  L.  Chuang 


Abstract— We  present  a  comprehensive  model  for  the  calcu¬ 
lation  of  the  bandedge  profile  of  both  the  Ini-xGa^AsyPi-^ 
and  Ini-x-yGaiAlyAs  quantum-well  systems  with  an  arbitrary 
composition.  Using  a  many-body  optical  gain  model,  we  compare 
the  measured  net  modal  gain  for  both  material  systems  with 
calculations  from  the  realistic  band  structure  including  valence 
band  mixing  effects.  Calibrated  measurements  of  the  side  light 
spontaneous  emission  spectrum  based  on  its  fundamental  relation 
to  the  optical  gain  spectrum  give  values  for  the  radiative  current 
density.  These  measurements  allow  us  to  extract  the  relationship 
between  total  current  density  and  carrier  density.  A  fit  of  this 
relation  yields  values  for  the  Auger  coefficient  for  each  material 
system. 

Index  Terms—  Auger  recombination,  gain,  InGaAsP, 
InGaAlAs,  semiconductor  lasers,  strained  quantum-weU 
lasers. 


I.  Introduction 

There  are  two  main  material  systems  used  to  fab¬ 
ricate  long-wavelength  semiconductor  lasers.  While 
the  Ini-.:,Ga^ASyPi^y  system  has  been  firmly  established 
as  a  candidate  for  high-performance  lasers  [1],  [2],  the 
Ini«ar-yGaxAlyAs  system  has  only  recently  begun  to  yield 
quality  results  [3]-[5],  The  interest  in  the 
system  stems  from  its  larger  conduction  band  offset  of 
AEc/AEg  =  0.7  compared  to  a  value  ofAEc/AEg  =  0.4 
for  Ini_xGaxASyPi_y.  This  larger  conduction  band  offset  had 
been  predicted  to  result  in  better  electron  confinement  in  the 
conduction  band  and,  therefore,  a  higher  temperature  stability. 

Previous  reports  have  been  given  on  the  model  for  the 
band  structure  of  these  two  systems  [6],  theoretical  calcu¬ 
lations  comparing  the  systems  in  terms  of  threshold  current 
density  and  differential  gain  [7],  as  well  as  the  tempera¬ 
ture  dependence  of  the  threshold  current  [8].  There  have 
also  been  measurements  of  the  contributions  of  the  various 
recombination  mechanisms  for  each  system  [9].  However, 
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little  has  been  published  directly  comparing  the  gain  and 
recombination  mechanisms  of  these  two  systems,  or  comparing 
the  existing  many-body  optical  gain  models  with  experimental 
measurements. 

In  this  paper,  we  present  a  comprehensive  model  for  the 
calculation  of  the  band  structure  of  samples  made  of  either 
of  these  systems  with  an  arbitrary  composition.  Then,  using 
a  many-body  model  for  the  optical  gain  and  spontaneous 
emission,  we  compare  the  model  to  experimental  data  for 
each  material  system.  Also,  from  simultaneous  measurements 
of  both  the  spontaneous  emission  spectrum  from  the  side  of 
the  lasers  and  the  gain  spectrum  from  the  facet  emission,  we 
are  able  to  compare  the  relative  contributions  of  the  different 
recombination  mechanisms  in  each  laser  system. 

In  Section  II,  we  present  our  model  for  the  calculation  of 
the  bandedge  of  each  material  system  and  the  quantum-well 
(QW)  band  edge  discontinuities,  followed  by  the  model  for 
the  optical  gain  and  spontaneous  emission.  In  Section  III,  we 
describe  the  structure  of  the  lasers  as  well  as  our  experimental 
measurements.  Section  IV  describes  the  comparison  of  the 
optical  gain  data  with  calculations  using  our  model,  and  the 
extraction  of  the  relative  contributions  of  the  recombination 
mechanisms  in  each  system.  Finally,  a  brief  conclusion  is 
given  in  Section  V. 

II.  Theoretical  Model 

A.  Calculation  of  Bulk  Bandgap  with  Strain 

To  obtain  most  parameters  for  both  the  Ini^^^Gdi^AsyPi-y 
and  Ini..a._yGax  AlyAs  material  systems,  a  linear  interpolation 
between  the  parameters  of  the  relevant  binary  semiconductors 
is  used.  The  interpolation  formulas  for  all  physical  parameters 
P  used  in  the  calculation  of  the  band  edge,  except  for  the 
bandgap,  are  given  as  [10],  [11] 

i^(Ini  ^j-Ga^  ASyP  i  -.y) 

=  P(GaAs)xj/  -f  P(GaP)x(l  -  y)  -I-  P(InAs)(l  —  x)y 
+  P(InP)(l  -x)(l  -y)  (]) 

P(Ini  .j.yGa^c  Aly  As) 

=  P{InAs)(l  ~  X  -  y)  4-  P(GaAs)x  +  P(AlAs)y.  (2) 

The  material  parameters  of  the  binary  semiconductors  can  be 
found  in  Table  I. 
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TABLE  I 


Parameters  for  the  Calculation  of  Strain  a.nd  Bandgap  Energy  for  the  Ini  -x 

GSx  ASy  Pi 

—  y  AND  In]_x*-*y 

GaxAIyAs  Material  Systems  [10],  [11] 

Parameter 

Symbol(unit) 

GaAs 

InAs 

InP 

GaP 

AlAs 

Lattice  Constant 

«(A) 

5.6533 

6.0584 

5.8688 

5.4505 

5.660 

Elastic  Stiffness  Constant 

Cn  ( 1 0”  dy  n/cm^ ) 

11.879 

8.329 

10.11 

14.05 

12.5 

Elastic  Stiffness  Constant 

Ci2(10*^dyn/cm’) 

5.376 

4.526 

5.61 

6.203 

5.34 

Hydrostatic  deformation  potential 

for  conduction  band 

«c(eV) 

-7.17 

-5.08 

-5.04 

-7.14 

-5.64 

for  valence  band 

1.16 

1.00 

1.27 

1.70 

2.47 

Shear  deformation  potential 

for  valence  band 

6(eV) 

-1.7 

-1.8 

-1.7 

-1.8 

-1.5 

Valence  band  parameter 

7l 

6.8 

20.4 

4.95 

4.05 

3.45 

7j 

1.9 

8.3 

1.65 

0.49 

0.68 

73 

2.73 

9.1 

2.35 

1.25 

1.29 

Electron  effective  mass 

m,/m. 

0.067 

0.023 

0.077 

0.25 

0.15 

Heavy-hole  effective  mass 

0.50 

0.40 

0.60 

0.67 

0.79 

The  one  exception  to  the  linear  interpolation  is  the  formula 
for  the  unstrained  bandgap.  For  each  material  system,  this 
quantity  is  given  as 

Ini-.,GaxASyPi_y£p(x,  y) 

=  1.35  +  0.642X  ^  l.lOly  +  0.758x^  4-  O.lOly^ 

-  0.159xy  -  0.28x2y  +  O.lOQxy^  (eV)  (3) 

Ini-j«yGaxAlj,AsEp(x,  y) 

=  0.36  +  2.093y  +  0.629x  +  0.577y^  -f  0.436x^ 

4-  1.013xy  -  2.0xy(l  -  x  -  y)(eT^).  (4) 

The  effects  of  strain  are  calculated  in  the  following  way. 
First,  the  strain  in  the  plane  of  the  epitaxial  growth  is 

Co  ~  a  ' 

e  —  6x1  —  6j,y - - — :  (5) 

a 

where  a  is  the  lattice  constant  of  the  quaternary  epitaxial  layer 
and  a©  is  the  lattice  constant  of  the  substrate  which  is  assumed 
to  be  InP.  The  strain  in  the  perpendicular  direction  can  be 
expressed  as 

=  -2^£  (6) 

where  Cn  and  C12  are  the  elastic  stiffness  constants.  The 
conduction  band  is  shifted  by  the  energy 

and  the  valence  bands  are  shifted  by 

6Ehh[x-  y)=  -P(-Qc  (8) 

6E,h{x.y)=  -Pc  +  Q,  (9) 

where 

Pf  =  —  (6x1  +  €yy  4-  Czz)  (10) 


6Ec{x,  y)  =  ac(£„  +  £yy  +  e. 


»)  =  2ac^ 


where  and  are  the  conduction-band  and  valence-band 
hydrostatic  deformation  potentials,  and  b  is  the  valence-band 
shear  deformation  potential. 

The  strained  bandgaps  can  then  be  expressed  as 


Ec-hkix,  y)  =Eg{x,  y)  +  SEc{x,  y)  -  SEhhix,  y)  (14) 
and 

Ec-th{x,  y)  =  Eg{x,  y)  +  6Ec{x,  y)  -  6Eih{x,  y).  (15) 


Fig.  1(a)  shows  a  contour  plot  relating  x  and  y  to  the  strain 
and  energy  gap  of  the  Inj^xGaxASyPi^y  system  grown  on 
an  InP  substrate.  Fig.  1(b)  shows  the  same  information,  but 
inverted,  so  that,  given  a  strain  and  a  bandgap,  the  Ga  and 
As  mole  fractions  can  be  obtained.  Fig.  2(a)  and  (b)  shows 
the  results  for  the  Ini-x-yGaxAIyAs  system  with  an  InP 
substrate.  In  all  of  these  figures,  the  bandgap  with  the  smallest 
energy  is  presented  (i.e.,  c-hh  for  compressive  strain,  c-lh  for 
tensile  strain)  and  the  bandgap  shift  due  to  the  strain  effects  is 
taken  into  account.  Note  that  these  results  are  different  from 
a  conventional  plot  where  the  “unstrained’’  bandgap  energy  is 
usually  presented. 

B.  Calculation  of  Band  Offset 

When  modeling  the  bandstruciure  for  a  QW,  the  relative 
alignment  of  the  band  edges  of  the  well  and  barrier  material 
is  very  important.  Unfortunately,  there  are  few  data  for  arbi- 
trar>’  compositions  of  these  two  quaternary'  material  systems, 
especially  for  the  strained  cases.  Here  we  present  and  compare 
two  models  for  the  band  alignment  of  these  systems. 
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(a) 


(b) 

Fig.  1.  Strain  and  strained  bandgap  energy  versus  Ga  (i)  and  As  (i/)  mole 
fractions  for  the  Ini— xGax Asy Pi  — y  material  system  with  (a)  contours  of 
constant  strain  and  bandgap  energy  and  (b)  contours  of  constant  x  and  y. 


AIAs 


Fig.  2.  Strain  and  strained  bandgap  energy  versus  Ga  (jr)  and  A1  (y)  mole 
fractions  for  the  Ini-x-yGaxAlyAs  material  system  with  (a)  contours  of 
constant  strain  and  bandgap  energy  and  (b)  contours  of  constant  x  and  y. 


J)  Model-Solid  Theory^  [JO],  [12]:  The  valence  band  posi¬ 
tion  of  a  quaternary  is  given  by 


£v(x,  y)  = 


-^v,av(x,  y)  H  +  6Eh,k{^^  y), 

for  hh  (compressive  strain) 

*^v,av(x,  y)  -f  — “2’ y), 

.  for //i  (tensile  strain)  ; 


(16) 


where  Ey^av{x,  y)  is  the  average  valence  subband  energy  and 
A  is  the  spin-orbit  spiil-off  band  energy.  These  values  are 
obtained  by  a  linear  interpolation  [see  (1)  or  (2)]  of  the  binary 
values  listed  in  Table  II,  The  conduction  band  position  may 
be  calculated  by  simply  adding  the  strained  bandgap  energy 
to  the  valence  band  position 


E'vix,  y)  4-  Ec~-kh{x^  y)j 
for  hh  (compressive  strain) 

■^v{x.  y)  +  Ecr-iK{,X)  y), 

.  for  Ih  (tensile  strain). 


(17) 


The  conduction  band  offset  is  given  by 


A^C  ^  -  El 

AE,  El  ~  £- 


(18) 


where  E'^;'  and  E^  are  the  valence  band  positions  in  the  well 
and  barrier  materials,  respectively,  and  E^'  and  E^  are  the 
strain  adjusted  band  gaps  (Er^hh  for  compressive  strain  and 
Ec^ih  for  tensile  strain)  for  the  well  and  barrier  materials. 


2)  Harrison’s  Model  [13]:  The  position  of  both  the  con¬ 
duction  and  valence  bands  are  determined  by 


r  (^>  y)  +  SEhh{x,  y), 

E  (x  7/)  —  J  (compressive  strain) 

\  (x,  y)  +  SEihix,  y), 

\  for  Ih  (tensile  strain) 

Ecix,  y)=E»  +  6Ecix,  y) 


(19) 

(20) 


where  E^(x,  y)  and  E^{x,  y)  are  obtained  by  a  linear 
interpolation  of  the  binary  parameters  found  in  Table  II,  and 
^E)hh{x^  y),  6Eih{x^  y),  and  6Ec{x^  y)  are  the  strain-induced 
energy  shifts  given  in  (7)-(9).  The  superscript  “i/”  refers 
to  Hamson’s  model.  The  conduction  band-edge  discontinuity 
may  then  be  calculated  as 


AE,  _ 

^^9  ~  -  Elf’’’)  +  -  E^''^) 

where  the  superscripts  w  and  6  indicate  the  well  and  barrier 
materials,  respectively. 

It  should  be  noted  that  this  method  simply  is  meant  to  give 
the  parameter  AEdAEg  which  may  be  used  to  determine 
the  alignment  of  the  well  and  barrier  materials.  The  difference 
between  Ey{x^  y)  and  Ec{x^  y)  from  (19)  and  (20)  should  not 
be  used  to  calculate  the  bandgap  of  the  quaternary  material. 
Rather,  (3)  and  (4)  should  be  used. 

Fig.  3  shows  a  comparison  of  the  model-solid  theory 
and  Harrison’s  model  for  the  calculation  of  the  conduction 
band  offset  for  both  lattice-matched  Ini-xGaxASj/Pi.y  and 
Ifii-i-yCaiAlyAs.  Also  shown  are  empirical  curves  (solid) 
based  on  e.xperimental  measurements  of  these  quantities 
[14],  [15).  For  the  Inj-a-CaxASyPi-y  system,  Harrison’s 
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TABLE  II 

Parameters  for  the  Calculation  of  Band  Alignment  of  the  Ini-^CaxAsyPi-y  and 
Ini  ^x-»yGax  AlyAs  Material  Systems  Using  the  Model-Solid  Theory  [12J  and  Harrison's  Model  [13] 


Parameter 

Symbol(unit) 

GaAs 

InAs 

InP 

GaP 

AlAs 

Model-solid  theory: 

Average  valence  band  position 

£v,,.(eK) 

-6.92 

-6.67 

-7.04 

-7.40 

-7.49 

Spin-orbit  spIit-ofT  energy 

A(«V) 

0.34 

0.38 

0.11 

0.08 

0.28 

Harrison’s  model: 

Conduction  band  position 

E»{eV) 

1.53 

0.801 

1.35 

2.352 

2.5255 

Valence  band  position 

E!/(eV) 

0.111 

0.441 

0.00 

-0.388 

-0.4245 

1,0 

0.8 

K 

0.2 


0.0 


'  Empiric*! 

^  -  Htrriion^  model 
* . Model -Solid  Theory 


In.  Ga  As  P 

l  x  *  y  l-y 


0.0  0.2 


0.4  0.6  0.8 

y 

(a) 


1.0 


y 

(b) 


Fig.  3.  The  conduction  band  offset  for  (a)  Inj  ,»xCaxAsyPi_y  and  (b) 
Ini-x-yGaxAlyAs  lattice-matched  to  InP.  The  solid  lines  show  empiri¬ 
cal  results  [14],  [15]  obtained  from  experimental  data.  Dashed  lines  show 
calculations  using  the  model-solid  theory  and  Harrison’s  model. 


model  is  excellent  while  the  model-solid  theory  gives  the 
opposite  trend,  only  being  good  near  the  InGaAs  boundary. 
For  Ini^aj.yGajcAlyAs,  both  models  have  some  error,  but 
Hamson’s  model  is  again  closer  to  the  empirical  curve.  For 
the  calculations  of  the  band  structure  in  this  paper,  we  therefore 
have  chosen  to  use  Harrison’s  model. 


C.  Calculation  of  Band  Structure  and  Gain 

The  optical  gain  with  many- body  effects  is  calculated  by  a 
non-Markovian  model  [16].  The  plasma  screening,  bandgap 
renormalization  (BGR),  and  the  excitonic  or  the  Coulomb 
enhancement  (CE)  of  the  interband  transition  probability  are 
included  in  the  model.  The  valence  band  structure  is  calculated 
by  using  the  block-diagonalized  3x3  Hamiltonian  based  on  the 
k'p  method  [17].  Aside  from  the  conduction  band  effective 
mass  for  the  Ini-xCa^ASyPi-y  system,  parameters  used  in 
the  calculation  are  interpolated  using  (I)  and  (2)  and  the 


binary  parameters  summarized  in  Table  I.  The  conduction  band 
effective  mass  for  Ini«xGaxASyPi_y  is  given  as 

7ne(Ini^a.GaxASyPi-.y) 

=  0.08  -  OAiex  -h  0.026y  -  0.059xy 
-h  (0.064  -  0.02x)y2  -f  (0,06  -h  0.032y)x^  (22) 


The  optical  gain  spectra  are  related  to  the  spontaneous 
emission  spectra  from  the  detailed  balance  [18],  [19]  between 
absorption  and  emission  of  photons.  This  means  that  there  is 
a  transparency  point  in  the  gain  spectra  that  is  determined  by 
the  Fermi  level  separation  AF  that  suggests  the  carriers  are 
in  a  quasi-equilibrium  distribution.  The  following  fundamental 
relationship  [18],  [19]  is  then  used  to  calculate  the  optical  gain: 


9{^)  = 

^sp(^)  = 


1  -  exp 


TT^C^ 


'hcj  -  AF 


kT 


(u) 


(23) 

(24) 


where  AF  is  the  quasi-Fermi  level  separation,  k  is  Boltz¬ 
mann’s  constant,  T  is  the  device  temperature,  c  is  the  speed 
of  light,  n  is  the  effective  refractive  index  of  the  optical 
waveguide  mode,  and  rspon{uj)  is  the  spontaneous  emission 
rate.  The  spontaneous  emission  coefficient  9sp{u^)  taking  into 
account  the  non-Markovian  relaxation  and  the  many-body 
effects  is  given  by  [16] 


x{l-/m(fc||))(l-Re9*„) 


ReI(E,.n(A:||))  -Img,„ImI(£,„.(fell)) 

{l-Reqk,)^  +  {lmgk,r 

(25) 

where  uj  is  the  angular  frequency,  fig  is  the  vacuum  permeabil¬ 
ity,  €  is  the  dielectric  constant,  fc||  is  the  in-plane  wavevector, 
L,  is  the  well  thickness,  |M/„p  is  the  momentum  matrix 
element  in  the  strained  QW,  //  and  are  the  Fermi  function 
for  the  conduction  band  states  and  the  valence  band  states, 
and  h  is  the  reduced  Planck’s  constant.  The  indices  /  and  m 
denote  the  electron  states  in  the  conduction  band  and  the  heavy 
hole  (light  hole)  subband  states  in  the  valence  band.  Also, 
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TABLE  HI 

A  Summary  of  the  Structure  of  the  Three  Laser  Systems  Stl^died 


Sample  A  Well: 

Barrier: 

Number  of  wells 

Material 

Strain 

Width 

PL  wavelength 

Material 

Strain 

Width 

PL  wavelength 

7 

InGa.^sP 

—0.9%  compressive 

80A 

1.58;im 

InGaAsP 

lattice  matched 

80A 

1.255pm 

SCH  width:  700A 

Stripe  width:  1,2pm 
Cavity  length:  350pm 

Sample  B  Well: 

Number  of  wells 

5 

SCH  width:  lOOOA 

Material 

InGaAlAs 

Stripe  width:  1.5pm 

Strain 

—0.78%  compressive 

Cavity  length:  890pm 

Width 

54A 

PL  wavelength 

1.56pm 

Barrier: 

Material 

InGaAlAs 

Strain 

lattice  matched 

Width 

57A 

PL  wavelength 

1.21/im 

Sample  C  Well: 

Number  of  wells 

5 

SCH  width:  600A 

Material 

InGaAlAs 

Stripe  width:  1.46pm 

Strain 

lattice  matched 

Cavity  length:  638pm 

Width 

86A 

PL  wavelength 

1.56pm 

Barrier: 

Material 

InGa.MAs 

Strain 

lattice  matched 

Width 

50A 

PL  wavelength 

1.21pm 

^/m(fc||)  —  Ef{k\i)  —  £?^(A:||)  +  Eg  +  AEsx  +  AjEch 
is  the  renormalized  transition  energy  between  electrons  and 
holes,  where  Eg  is  the  bandgap  of  the  material,  and  A^sx 
and  AEch  are  the  screened  exchange  and  Coulomb-hole 
contributions  [21]  to  the  bandgap  renormalization.  The  factor 
accounts  for  the  excitonic  or  Coulomb  enhancement  of 
the  interband  transition  probability  [20],  [21].  The  lineshape 
function  is  Gaussian  for  the  simplest  non-Markovian  quantum 
kinetics  and  is  given  by  [16] 


(26) 


and 


i  ]dt. 


(21) 


The  iniraband  relaxation  time  Tin  and  correlation  time  are 
assumed  to  be  constant  and  used  as  fitting  parameters. 


III.  Experiments 


A.  Laser  Structures 

In  our  experiments,  we  perform  measurements  using  lasers 
with  three  types  of  active  regions.  First,  we  tested  a  buried 
heterostructure  laser  with  five  lattice-matched  QW*s  made  of 
the  Ini_a.^yGaa.AIyAs  material  system,  as  well  as  a  buried 
heterostructure  laser  with  five  —0.78%  compressively  strained 
QWs  of  the  same  system.  Details  on  the  growth  and  fab¬ 
rication  of  these  devices  can  be  found  elsewhere  [5].  Also 
tested  was  a  buried  heterosiructure  laser  with  seven  -0.9% 
compressively  strained  QW’s  made  of  the  Ini_*GaxASj,Pi_y 
material  system.  Table  IE  contains  detailed  information  for  the 
structure  of  each  active  region. 

B.  Experimental  Procedure 

The  facet  emission  spectrum  of  each  laser  w'as  measured  for 
various  currents  below  lasing  threshold.  The  device  under  test 
is  mounted  on  a  heat  sink  and  is  held  at  a  constant  temperature 
of  25  using  a  thermoelectric  cooling  unit.  The  current 
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10  the  device  was  supplied  by  an  ILX  Lightwave  LDX-3412 
current  source  and  monitored  to  the  hundredth  of  a  milliamp 
using  a  Fluke  87  III  digital  multimeter.  The  facet  output 
was  collected  using  an  E-TEK  Dynamics  laser  optical  fiber 
interface  with  a  built-in  polarizer  to  control  the  polarization 
of  the  measured  emission.  An  optical  isolator  was  used  to 
prevent  optical  feedback  effects  and  the  light  was  monitored 
using  an  HP  7095  IB  optical  spectrum  analyzer  with  a  spectral 
resolution  of  0.08  nm.  The  net  modal  gain  for  these  devices 
was  obtained  using  the  well  known  Hakki-Paoli  method  [22] 
in  which  the  ratio  of  the  maximum  /max  and  minimum  /min 
of  the  amplified  spontaneous  emission  spectrum  yields  the  net 
modal  gain  Guet  by  the  following  relations: 


where  g  is  the  material  gain  for  one  well,  T{—Nu,T^)  is  the 
optical  confinement  factor  per  well  (T^^)  multiplied  by 
the  number  of  wells,  a,-  is  the  intrinsic  loss,  L  is  the  cavity 
length,  and  Ri  and  R2  are  the  facet  power  reflectivities. 

The  side  light  spontaneous  emission  [24]  of  the  same  device 
was  also  monitored  for  currents  up  to  lasing  threshold.  In  this 
measurement,  a  cleaved  multimode  fiber  is  brought  up  close 
to  the  side  of  the  device,  usually  near  the  middle  in  order  to 
avoid  scattered  light  from  the  laser  facets.  First,  using  the  same 
equipment  as  for  the  gain,  the  spontaneous  emission  spectrum 
is  taken  using  a  resolution  bandwidth  of  10  nm  for  each  of 
the  currents  at  which  the  facet  emission  was  measured.  Next, 
the  current  is  controlled  using  an  HP  4145  semiconductor 
parameter  analyzer  in  order  to  achieve  currents  as  low  as 
20  //A.  The  current  is  swept  and  the  light  captured  by  the 
fiber  is  then  fed  into  an  HP  8 153 A  Lightwave  Multimeter  to 
obtain  the  integrated  spontaneous  emission  power.  For  these 
measurements,  the  width  of  the  active  region  is  sufficiently 
small  (1.2-1.46  ^rm)  so  that  amplification  of  the  SE  spectra 
is  negligible. 


IV.  Analysis  and  Discussion 

The  measured  net  modal  gain  obtained  from  the  amplified 
spontaneous  emission  (ASE)  spectrum  using  (28)  for  each  of 
the  three  material  systems  studied  is  presented  in  Fig.  4  for 
several  currents  up  to  threshold.  Typical  mode  spacings  in 
the  ASE  spectrum  for  samples  A,  B,  and  C  are  about  0.93, 
0.38,  and  0.51  nm,  respectively.  To  check  for  the  effects  of 
spectrometer  resolution  on  the  measured  net  modal  gain,  the 
min/sum  method  [23],  which  is  less  sensitive  to  the  resolution 
limits  of  the  measurement  (although  more  sensitive  to  noise), 
was  also  used  to  extract  the  gain  spectra.  The  resulting  spectra 
were  nearly  identical  except  for  values  closest  to  the  threshold 
gain  value  where  the  Hakki-Paoli  method  underestimates  the 
gain,  but  these  differences  give  less  than  a  59c  error  on  the 
value  of  the  modal  gain. 

Also  presented  in  Fig.  4  are  calculated  spectra  using  the 
model  presented  above  for  the  material  gain  of  the  three 
systems.  The  iniraband  relaxation  time  Tin  and  the  correlation 


1480  1520  1560  1600  1640 

Wavelength  (nm) 


(c) 

Fig.  4.  The  measured  modal  gain  (solid  lines)  and  calculated  (dashed  lines) 
values  from  the  electronic  band  structure  for  (a)  lni-*Gax  ASyPi-y  with 
-0.9%  compressive  strain  (Sample  A),  (b)  Ini-r-yGar  AlyAs  with  -0.78% 
compressive  strain  (Sample  B).  and  (c)  Im-x-yGax  AlyAs  lattice  matched 
(Sample  C).  The  device  parameters  arc  listed  in  Table  III 


lime  Tc  used  for  the  fits  are  given  in  Table  IV.  For  the  screened 
exchange  AEsx  and  Coulomb-hole  AEch  contributions  to 
the  bandgap  renormalization,  [21,  eqs.  (27)  and  (28)]  were 
used  with  a  C  parameter  of  1.1,  1.7,  and  2.0  for  Samples 
A,  B,  and  C,  respectively.  Each  of  these  parameters  affect  the 
peak  gain  position  and  overall  shape  of  the  spectra.  They  were 
selected  to  yield  the  best  overall  fit,  balancing  these  criteria. 
The  net  modal  gain  is  obtained  using  an  intrinsic  loss  which 
is  independent  of  wavelength,  and  a  confinement  factor  with 
a  linear  dependence  on  wavelength,  with  a  slope  equal  to 
1.66x10"^  nm”^  (estimated  from  the  wavelength-dependent 
refractive  index  and  waveguide  theory  for  the  fundamental 
mode)  and  the  vertical  intercept  treated  as  a  parameter  to  fit 
the  data. 

Fig.  5(a)  shows  a  plot  of  the  fitted  values  of  the  intrinsic 
loss  as  a  function  of  the  injected  current  for  each  device. 
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TABLE  IV 

Parameters  Used  in  and  Obtained  from  the  Theoretical  Fits 


Pwamcter 

Symbol 

Sample  A 

Sample  B 

Sample  C 

Fit  of  the  mod&l  g&in: 

iDterbftnd  rel&xAtion  time 

r»(/«) 

40 

60 

85 

CcrreUtioQ  time 

r.(/<) 

20 

20 

20 

Qu&si-Fermi  level  »epAr4tion 

AF(tV) 

O.SOl  (4mA) 

0.804  (5mA) 

0.806  (6mA) 

0.$13  (6mA) 

0.815  (7mA) 

0.816  (8mA) 

0.822  (8mA) 

0.823  (9mA) 

0.823  (10mA) 

0.828  (10mA) 

0.829  (12mA) 

0.827  (11.7mA) 

XV&nsform&tioD  of  tpont&neoui  emission: 

Scaling  coefficient 

Q(cm->nm-*(mM')->) 

6.11  X  10“^  ±0.22  X  10"' 

1.81  X  10-^  ±  0.08  X  10-^ 

8.2  X  10-*  ±  0.02  X  10-» 

Effective  index  of  refraction 

n 

3.3  ±  0.15 

3.3  ±0.15 

3.3  ±0.15 

Optical  confinement  factor 

r(-//.r.) 

0.15  ±0.02 

0.06  ±0.02 

0.11  ±  0.02 

Optical  confinement  (actor  per  well 

r. 

0.021  ±  0.0029 

0.012  ±0.004 

0.022  ±  0.004 

Fit  of  the  J,N  relationship: 

Coupling  coefficient 

1.S6  X  10-“  ±  0.01  X  10*“ 

1.64  X  10*“  ±  0.OS  X  10-“ 

1.21  X  10-“  ±  0.01  X  10-“ 

Monomolecular  coefficient 

Mr') 

1.1  X  10*  ±0.2  xlO* 

7x10' ±2x10' 

6.7x10' ±3x10' 

Auger  coefficient 

1.1  X  10->’±0.1S  X  10-” 

1.4  X  10-“  ±  0.3  X  10-” 

3  X  10-”  ±  2  X  10-” 

Injection  efficiency 

0.92  ±  O.IT 

0.95  ±  O.IS 

0.4$  ±0.15 

Characteristic  temperature  (measured) 

T.i'K) 

36 

53 

51.4 

The  observed  increase  in  the  loss  is  consistent  with  previous 
reports  [25]  and  can  be  explained  by  the  increase  in  free-carrier 
absorption  with  an  increasing  carrier  density.  The  fitted  optical 
confinement  factor  used  for  the  three  devices  is  shown  in  Fig. 
5(b).  Values  for  the  optical  confinement  factor  per  well  (Fy,) 
and  the  overall  confinement  factor  (F  =  N^T^v)  are  listed 
in  Table  IV.  These  fitted  values  are  in  good  agreement  with 
typical  values  for  InP-based  lasers  [26]  when  scaling  due  to 
well  width  is  taken  into  account.  Finally,  Fig,  5(c)  shows  a 
plot  of  peak  net  modal  gain  versus  the  total  current  density  for 
each  laser  system.  The  fits  of  the  gain  sp^tra  of  these  three 
material  systems  allow  us  to  obtain  the  density  of  carriers  N 
in  the  QW’s  corresponding  to  the  measured  injection  currents. 
For  each  system,  other  samples  from  the  same  wafer  exhibit 
similar  gain  spectra. 

The  same  model  used  to  create  the  gain  spectrum  was  also 
used  to  calculate  the  spontaneous  emission  rate  for  each  of  the 
material  systems  studied.  The  radiative  current  density  can  be 
related  to  the  spontaneous  emission  rate  as 

f  ~  (30) 

Jo 

where  c  is  the  charge  on  an  electron,  is  the  number  of 
wells,  Lw  is  the  width  of  a  QW,  rspon(A)  is  the  spontaneous 
emission  rate,  Beff  is  the  effective  radiative  recombination 
coefficient,  and  N  is  the  carrier  density.  Jrad  can  be  obtained 
from  our  calculated  spontaneous  emission  spectrum  using  (30) 
with 


In  order  to  validate  these  calculations  from  the  band  struc¬ 
ture,  it  is  important  to  be  able  to  obtain  the  radiative  current 
density  experimentally.  As  described  above,  we  measure  the 
spontaneous  emission  power  spectra  from  the  side  of  each 
laser  device.  The  spontaneous  emission  power  captured  by  the 
optical  fiber  can  be  related  to  the  spontaneous  emission  rate  by 

P(A)  =  A'i^r.po„(A)(iAl/  (32) 


where  P(A)  is  the  measured  power  in  the  fiber,  K\  is 
a  coupling  coefficient  which  includes  the  fraction  of  the 
spontaneous  emission  captured  by  the  optical  fiber,  h  is 
Planck’s  constant,  c  is  the  speed  of  light,  d\  is  the  specu-ometer 
resolution,  and  V  is  the  active  region  volume.  Since  K\  is  an 
unknown  quantity,  it  is  not  straightforward  to  obtain  a  value 
for  the  spontaneous  emission  rate  so  that  (30)  can  be  used  to 
extract  A  calibration  can  be  accomplished  based  on  the 
fundamental  relationship  between  the  material  gain  and  the 
spontaneous  emission  rate  [(23)  and  (24)),  rewritten  as 


P(A)  = 


A^ 


Sitcri^ 


1  —  exp 


fhc/X-AF 


keT 


)] 


rspon(A)  (33) 


where  n  is  the  effective  refractive  index,  AF  is  the  quasi- 
Fermi  level  separation  energy,  ks  is  Boltzmann’s  constant, 
and  T  is  the  temperature.  The  material  gain  per  well  can  be 
related  to  the  measured  net  modal  gain  as 


Gnet  =  rp(A)  -  Q,  (34) 


r.pon(A)  = —j^ff.pCA)  (31) 

where  psp(A)  is  given  in  (25),  A  =  2itc/w,  c  is  the  speed  of 
light,  and  n  is  the  refractive  index  of  the  QW.  The  calculated 
curves  for  of  each  sample  are  shown  in  Fig.  6. 


where  F  is  the  optica]  confinement  factor  and  o,  is  the  intrinsic 
loss  for  the  optical  waveguide  mode.  Combining  (32)-(34), 
we  obtain 

G„c.(A)  =  QA*^1  -  exp  jp(A)  -  a,  (35) 
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Fig.  5.  (a)  The  optical  confinement  factor  and  (b)  intrinsic  loss  used  in  the 
fits  of  the  modal  gain  for  Samples  A,  B.  and  C  in  Fig.  4.  (c)  The  peak  net 
modal  gain  versus  the  iota!  current  density. 
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Fig.  6.  The  calculated  total  radiative  current  density  through  all  wells  [see 
(30)]  as  a  function  of  carrier  density  (lines)  compared  to  the  extracted  total 
radiative  current  density  from  the  calibrated  spontaneous  emission  spectrum 
measurements. 


where 


r 

8irK\hc^n^  dW 


(36) 


Fig.  7.  (a)  The  measured  spontaneous  emission  spectra  for  the  ~0.78% 
compress! vely  strained  Ini-x-yGaxAIy  As  laser  system  (Sample  B  in  Table 
III)  and  (b)  the  measured  net  modal  gain  spectra  (solid  lines)  compared  to  the 
transformed  spontaneous  emission  spectra  (dashed  lines)  using  (35). 


This  equation  relates  the  measured  net  modal  gain  to  the 
measured  spontaneous  emission  power  using  the  three  pa¬ 
rameters  Q,  AF,  and  a,-.  AF  and  Oi  are  not  independent, 
however,  since  the  net  modal  gain  value  measured  at  the 
wavelength  corresponding  to  AF  is  exactly  the  transparency 
level  and  equal  to  a,.  (We  assume  here  that  the  intrinsic  loss  is 
independent  of  wavelength.)  Hence,  the  short-wavelength  side 
of  the  measured  net  modal  gain  spectrum  gives  the  relationship 
between  these  two  quantities. 

Using  Q  and  AF  as  fitting  parameters,  (35)  may  then  be 
used  to  transform  the  measured  spontaneous  emission  power 
into  the  net  modal  gain  spectrum.  By  choosing  the  best 
fit  between  the  measured  and  transformed  net  modal  gains, 
an  accurate  value  for  Q  can  be  obtained.  As  an  example 
of  this  procedure,  Fig.  7(a)  shows  the  spontaneous  emis¬ 
sion  power  spectra  measured  for  the  compressively  strained 
Ini-.x-yGaa;AlyAs  system.  These  spectra  were  measured  for 
each  of  the  currents  that  the  gain  was  measured  at.  The  best  fits 
between  the  transformed  spontaneous  emission  and  actually 
measured  net  modal  gain  spectra  are  shown  in  Fig.  7(b).  The 
extracted  values  of  Q  are  about  1.8x10’"'^  and  vary  by  less 
than  4%  among  all  of  the  currents  measured,  showing  excellent 
uniformity  in  the  measurement  of  the  SE  spectra.  For  AF, 
the  extracted  values  for  each  system  deviate  from  the  values 
obtained  for  the  fitting  of  the  experimental  gain  curves  by  less 
than  0.1%  at  low  currents  where  the  difference  between  tw'o 
adjacent  curves  in  Fig.  4  represent  a  shift  of  about  1%  in  AF. 
At  higher  currents,  the  deviation  gets  as  large  as  0.45%  w’hile 
the  separation  between  curves  is  about  0.5%. 

Once  accurate  values  of  Q  have  been  extracted,  the  mea¬ 
sured  spontaneous  emission  power  spectra  can  be  calibrated 
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in  the  following  way  to  get  the  actual  spontaneous  emission 
rate.  By  solving  (32),  r,pon(A)  can  be  expressed  as 


where 


r.pon(A)  =  TXP{X) 


(37) 

(38) 


Note  that  this  procedure  has  eliminated  the  need  to  derive 
the  coupling  coefficient  Ki  explicitly.  The  data  points  in  Fig. 
6  were  obtained  by  substituting  (37)  and  (38)  into  (30)  and 
numerically  integrating  the  calibrated  data.  The  value  for  n 
was  assumed  to  be  3.3  with  an  error  of  0.15  for  each  of  the 
material  systems  measured.  The  optical  confinement  factor  was 
assumed  to  be  constant  with  a  value  equal  to  the  average  of  the 
range  used  for  the  gain  fits,  as  shown  in  Fig.  5(b).  The  error 
range  is  taken  from  the  values  of  the  confinement  factor  in 
Fig.  5(b)  at  the  edges  of  the  transformed  spontaneous  emission 
spectra.  This  amounts  to  an  error  of  about  ±12%  for  the  value 
of  the  confinement  factor,  leading  to  a  total  error  of  slightly 
over  ±20%  for  the  radiative  current  density. 

From  the  calculated  curves  of  the  radiative  current  density 
Jfid  as  a  function  of  carrier  density  N  in  Fig.  6,  the  radiative 
recombination  coefficient  for  each  material  system  is 
obtained  from  (30)  and  is  presented  in  Fig.  8.  This  pa¬ 
rameter  is  independent  of  the  well  width/number  product, 
so  it  serves  as  a  fair  comparison  of  the  relative  probabil¬ 
ity  of  the  radiative  recombination  of  an  electron-hole  pair 
for  each  structure.  We  see  that,  at  low  carrier  densities, 
the  compressively  strained  Ini_,Ga,ASyPj_j,  system  has  the 
largest  radiative  recombination  coefficient  followed  by  the 
compressively  strained  Ini_*_yGaxAlyAs  system  and  then  the 
lattice-matched  Inj_*_yGajAlyAs  system.  It  is  not  surprising 
that  the  strained  samples  are  more  efficient  since  strain  has 
been  shown  to  decrease  the  in-plane  effective  mass  in  the 
top  (hh)  valence  band,  leading  to  a  more  symmetric  injection 
of  electrons  and  holes  in  the  conduction  and  valence  bands 
[27].  Also,  the  heavy-hole  and  light-hole  subband  separation 
becomes  larger  with  strain,  leading  to  less  wasted  carriers. 
At  higher  carrier  densities,  all  three  systems  show  a  decrease 
in  the  radiative  recombination  coefficient.  This  is  especially 
pronounced  in  the  compressive  Ina_*GaxASyPi_y  system. 
This  decrease  is  probably  due  to  the  significant  population 
of  higher  energy  valence  subbands  while  the  higher  energy 
conduction  subbands  are  virtually  empty.  For  example,  in  all  of 
these  systems,  the  second  conduction  subband  has  less  than  3% 
of  the  total  carrier  population.  Meanwhile,  for  the  compressive 
Ini_xGaxASyPi_y  system,  the  first  two  heavy-hole  subbands 
(which  are  the  lowest  subbands)  have  a  spacing  of  36  meV 
and  a  distribution  of  approximately  60.1%  of  carriers  in  the 
first  heavy-hole  subband  and  29.6%  in  the  second  heavy- 
hole  subband.  Also,  for  the  compressive  Inj-j-yGa^AlyAs 
sample,  there  is  a  spacing  of  65  meV  between  the  first  two 
heavy-hole  subbands  and  a  carrier  disU'ibution  of  about  80% 
and  19.4%  for  the  first  and  second  heavy-hole  subbands, 
respectively.  This  difference  is  mostly  due  to  the  thinner 
well  of  the  compressive  Inj_x  yGa^AlyAs  sample.  Holes 
in  the  second  heavy-hole  subband  are,  therefore,  unable  lo 
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Fig.  8.  The  radiative  recombination  coefficient  as  a  function  of  earner 
density  extracted  from  the  calculated  total  radiative  current  density. 


significantly  contribute  to  radiative  recombination  processes 
since  they  can  only  recombine  with  electrons  in  the  second 
conduction  subband,  which  is  sparsely  populated  for  each 
laser  system  studied.  The  decrease  seen  in  the  compressive 
Ini_xGa*AsyPi^y  system  is  more  significant  due  to  the 
smaller  spacing  between  valence  subbands.  At  the  higher 
carrier  densities,  this  spacing  leads  to  an  even  larger  increase 
in  the  population  of  higher  energy  subbands,  and  thus  more 
wasted  carriers. 

To  study  the  major  recombination  mechanisms  in  each  of 
these  lasers,  it  is  convenient  lo  use  the  measured  spontaneous 
power  which  has  been  monitored  by  an  optical  power  meter. 
In  this  case,  the  data  are  easily  measured  over  a  large  current 
range,  and  the  photodeteclor  performs  the  integration  of  the 
spontaneous  emission  power  spectrum.  Fig,  9  shows  a  log-log 
plot  of  the  square  root  of  spontaneous  emission  power  versus 
the  total  injected  current  for  each  device  studied.  This  format 
was  chosen  [28]  to  give  an  indication  of  the  relative  con¬ 
tributions  of  each  recombination  mechanism  (monomolecular, 
radiative,  and  Auger).  The  value  of  the  x  axis  can  be  expressed 
as 


y/P  =  K2  y/BeffN  (39) 


where  K2  is  another  unknown  coefficient  which  includes 
the  fraction  of  spontaneous  emission  power  captured  by  the 
optical  fiber.  We  see  that  this  is  almost  directly  proportional 
to  Nt  but  not  perfectly,  due  to  the  carrier  dependence  of 
Bcff*  as  discussed  above.  This  deviation  is  relatively  small, 
however,  since  we  can  see  clear  regions  of  slope  one  and  slope 
three  which  correspond,  respectively,  to  regimes  where  I  is 
proportional  to  N  (monomolecular  mechanisms  dominate)  and 
(where  Auger  mechanisms  dominate).  The  region  of  slope 
two  could  indicate  that  the  radiative  processes  are  dominant, 
but  this  is  not  clear  since  there  would  be  a  point  at  which 
the  slope  is  two  due  to  the  transition  between  the  N  and 
regimes,  even  if  radiative  processes  were  missing.  These  issues 
will  be  clarified  in  the  following  analysis. 

It  is  necessary  to  transform  the  measured  data  into  a 
relationship  between  total  current  density  and  carrier  density. 
Current  density  is  easily  obtained  from  the  current  as 


wL 


(40) 


where  w  is  the  stripe  width  and  L  is  the  cavity  length.  To 
obtain  N,  however,  the  unknown  coefficient  K2  must  be 


to 
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Fig.  9.  The  square  root  of  the  measured  integrated  spontaneous  emission 
power  y/P  (ploned  along  the  x  axis)  versus  the  measured  injection  current 
along  the  vertical  axis.  The  x  axis  is  roughly  proportional  to  carrier  density, 
and  the  slope  of  the  lines  on  the  log-log  plot  give  the  power  law  for  the 
injection  current  and  carrier  density  relation. 


determined.  The  fits  of  the  net  modal  gain  with  our  theory 
give  a  relation  between  I  and  N  over  a  limited  range  of 
the  measured  data  presently  under  consideration.  Using  the 
calculated  dependence  of  Besr  on  the  carrier  density  for  each 
system,  and  the  measured  ^/P  value  corresponding  to  the 
currents  for  the  fitted  gain  data,  (39)  yields  a  value  of  K2 
for  each  material  system.  The  average  between  the  values 
obtained  for  each  of  the  I-N  pairs  are  used  for  our  calculations 
and  are  listed  in  Table  IV.  All  values  are  within  2%  of  the 
average  value.  Since  K2  should  be  independent  (for  the  most 
part)  of  injection  current,  it  is  fair  to  extrapolate  the  value 
of  K2  from  this  small  range  of  currents  to  the  entire  range 
of  data  under  consideration.  Using  these  K2  values  and  the 
calculated  Bef[{N),  it  is  straightforward  to  obtain  a  value  for 
N  corresponding  to  the  measured  >/P  for  each  device. 

The  extracted  relationships  between  Jtot  and  N  are  shown 
as  dashed  lines  in  Fig.  10(a)  for  the  Ini^xUaxASyPi^y  system 
and  in  Fig.  10(b)  for  the  two  Ini_x-yGaxAlyAs  systems.  The 
solid  lines  in  these  plots  correspond  to  the  best  fit  results  using 
the  relation 

J^^,  =  ^IM:!L[AN  +  BMN)N^  +  CN^]  (41) 

T^inj 

where  7?,n;  is  the  injection  efficiency  accounting  for  the 
fraction  of  carriers  which  make  it  into  the  QW’s,  is 
the  number  of  wells,  A  is  the  monomolecular  recombination 
coefficient,  Beff(N)  is  the  calculated  radiative  recombination 
coefficient,  and  C  is  the  Auger  recombination  coefficient.  A  is 
a  fitting  parameter,  but  is  only  important  for  very  low  current 
densities.  Its  value  is  determined  by  fitting  the  behavior  of 
the  curve  at  these  low  current  densities.  Tyjnj  and  C  are  the 
major  parameters  which  control  the  behavior  of  the  curve.  For 
a  set  value  of  the  penalty  function  for  the  fit  is  very 
sensitive  to  the  C  coefficient,  so  the  value  can  be  very  well 
determined.  However,  a  range  of  T?inj  will  give  decent  fits,  so 
the  extracted  C  value  has  some  uncertainty.  The  best  fit  values 
of  77-, nj,  A,  and  C  are  presented  in  Table  IV.  The  uncertainties 
on  TTinj  and  C  are  determined  by  selecting  those  values  which 
raise  the  penalty  function  20%  over  the  best  fit  value.  The 
injection  efficiency  for  the  unstrained  Sample  C  appears  to 
have  a  relatively  low  injection  efficiency. 

It  should  be  noted  that  previous  measurements  of  the  Auger 
coefficient  in  the  Inj -xGajASyPi  ^  material  system  have 


(a) 


(b) 


Fig.  10.  The  extracted  current  density/carrier  density  relation  (dashed 
lines)  as  well  as  the  best  fits  from  the  theory  (solid  lines)  are  shown 
for  (a)  the  Ini.-xGax  ASyPi-j,  laser  system  (Sample  A)  and  (b)  the  two 
Ini-x-yGax  Aly  As  systems  (Samples  B  and  C). 


given  values  between  2x  10“^^cm®s“^  and  cm^s"^ 

[29H31).  As  discussed  previously,  the  fact  that  a  signifi¬ 
cant  amount  of  holes  occupy  states  where  they  may  only 
recombine  by  nonradiative  processes  could  explain  why  the 
extracted  Auger  coefficient  for  our  sample  is  significantly 
higher  than  this  range.  This  indicates  that  the  Auger  coefficient 
can  be  minimized  by  selecting  the  strain  and  well  width 
to  increase  the  separation  of  the  valence  subbands.  For  the 
Ini-x-yGaxAlj,As  material  system,  a  reported  value  of  the 
Auger  coefficient  [9]  is  3.6x10“^^  cm^s*"^  which  is  smaller 
than  our  value  for  the  compressive  sample.  Our  procedure, 
based  on  a  comprehensive  gain  model  and  a  fitting  of  the  gain 
spectra  to  determine  N  and  Beff,  should  provide  a  reliable 
estimate  of  the  C  coefficient. 

These  extracted  values  are  used  in  Fig.  11  to  com¬ 
pare  the  relative  contributions  of  Jmono*  and  JAug 

for  the  compressive  Ini^x-yGaxAlyAs  and  compressive 
Ini^xGaxASyPi.>y  systems.  Each  device  is  dominated  by  the 
Auger  recombination  well  before  the  laser  reaches  threshold. 
However,  due  to  an  Auger  coefficient  about  10  times  larger, 
the  Auger  recombination  dominates  at  a  much  lower  point  and 
becomes  much  larger  in  the  Ini-xGaxASyPi«y  system.  This 
leads  to  a  much  larger  threshold  current  density.  The  To  for 
each  of  these  lasers,  measured  from  the  shift  in  threshold  over 
a  temperature  range  from  15®  to  60®,  is  also  presented  in  Table 
IV.  These  relative  values  agree  with  the  trend  in  the  extracted 
Auger  coefficients,  indicating  that  Auger  recombination  is  a 
major  factor  reducing  To  for  these  laser  systems  [28],  [32]. 
Finally,  it  is  interesting  to  note  that  the  radiative  current  density 
never  becomes  the  dominant  recombination  mechanism.  At 
the  injection  level  where  Auger  recombination  becomes 
significant,  the  radiative  current  density  is  about  the  same 
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Fig.  II.  The  relative  contributions  of  the  monomolecular,  radiative, 
and  Auger  recombination  currents  for  (a)  —0.9%  compressively  strained 
lni-.xGaxAsyPi-y  (Sample  A),  (b)  —0.78%  compressively  strained 
Ini-x-yCaxAly  As  (Sample  B),  and  (c)  lattice-matched  Ini-x-uCax  ALAs 
(Sample  C).  ^  ^ 

as  the  monomolecular  current  density  for  each  system  studied. 
V.  Conclusion 

In  this  paper,  we  have  presented  a  comprehensive  model  for 
the  calculation  of  the  band  edges  of  both  Ini-iGaxASyPi^y 
and  Ini-je.yGaxAlyAs  QW  lasers  with  an  arbitrary  com¬ 
position.  Using  a  many-body  optical  gain  model,  we  have 
compared  the  net  modal  gain  for  both  material  systems  with 
calculations  from  the  realistic  band  structure.  These  compar¬ 
isons  allow  us  to  obtain  carrier  densities  for  each  of  the 
measured  curreni.s.  From  a  calibrated  measurement  of  the 
spontaneous  emission  spectrum,  we  were  able  to  obtain  the 
radiative  current  density.  These  results  are  all  in  excellent 
agreement  with  our  calculations.  From  a  fit  of  the  extracted 
relationship  between  carrier  density  and  total  current  density, 
we  were  able  to  extract  values  for  the  Auger  coefficient  of 


each  material  system  studied.  The  relative  values  of  these  pa¬ 
rameters  correspond  very  well  with  the  observed  temperature 
performance  of  each  device. 

From  this  study,  it  is  difficult  to  make  a  strong  statement  as 
to  which  of  the  two  material  systems  is  superior.  Obviously, 
more  studies  on  a  variety  of  active  region  designs  would  be 
needed  to  gauge  all  of  the  relevant  factors  and  draw  such 
a  conclusion.  Our  measurements  seem  to  indicate  that  the 
InGaAlAs  system  is  at  least  as  good  as  and  probably  better 
than  the  InGaAsP  systems  in  terms  of  Auger  recombination 
and  temperature  stability.  Also,  care  should  be  taken  to  design 
the  QW’s  such  that  there  is  a  large  separation  between  the 
lowest  and  highest  level  subbands,  in  order  to  reduce  the 
carriers  wasted  by  nonradiative  recombination  and  to  reduce 
the  temperature  sensitivity  of  the  lasers. 
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A  Wavelength-Tunable  Curved  Waveguide  DFB 
Laser  with  Integrated  Modulator 

A.  Hsu,  S.  L.  Chuang, 

W.  Fang,  L.  Adams,  G.  Nykolak,  and  T.  Tanbun-Ek 


Abstract —  A  distributed-feedback  (DFB)  laser  with  a 
curved  waveguide  in  the  active  region  integrated  with  an 
electroabsorption  modulator  is  studied  experimentally  and 
theoretically.  The  modulator  controls  the  lasing  wavelength 
of  the  DFB  laser  by  acting  as  an  optical  phase  shifter.  Our 
model,  which  is  based  on  the  transfer-matrix  method,  is 
used  to  simulate  this  multi-section  device  with  the  curved 
waveguide,  self-consistently  including  the  effects  of  spatial 
hole-burning  (SHB).  The  model  explains  the  features  and 
wavelength-tuning  behavior  of  the  spectrum  and  shows  good 
agreement  with  experimentally  measured  spectra.  We  also 
show  theoretically  that  the  curved  waveguide  suppresses  the 
longitudinal  photon  density  profile  of  the  device  compared 
with  a  straight  waveguide  case,  which  implies  reduced  SHB 
effects. 

Keywords —  Distributed-feedback,  modulator,  curved 
waveguide,  wavelength-tunable,  spatial  hole-burning 

I.  Introduction 

The  study  of  wavelength-tunable  laser  sources  is  an  im¬ 
portant  topic  of  research,  emphasized  by  the  rapidly 
expanding  development  of  wavelength-division  multiplexed 
(WDM)  optical  communications  systems,  which  can 
presently  transmit  data  at  rates  of  over  1  Tb/sec  over  a 
single  optical  fiber. 

Research  in  multi-section  wavelength-tunable  lasers  has 
been  investigated  extensively  for  use  in  WDM  and  coherent 
optical  communications  systems  and  some  of  the  progress 
in  this  field  is  well  summarized  in  Ref.[l],  [2].  Originally, 
the  goal  of  developing  distributed-feedback  (DFB)  lasers 
integrated  with  electroabsorption  modulators  was  to  im¬ 
prove  high-speed  performance  by  reducing  frequency  chirp¬ 
ing  of  DFB  lasers  by  external  modulation  [3]- [9].  Recently, 
these  devices  have  also  been  used  as  wavelength-tunable 
sources  [10],  [11].  The  use  of  a  chirped  grating  pitch  [12]- 
[14]  h^ls  been  investigated  as  a  way  to  acheive  a  distributed 
phase  shift  in  the  grating.  Several  groups  have  imple¬ 
mented  this  chirped  grating  pitch  using  a  bent  or  curved 
waveguide  to  acheive  wavelength  tuning  or  high  power  op¬ 
eration  and  have  shown  promising  results  [15]-[17].  In  pre¬ 
vious  work,  we  have  modeled  a  DFB  laser  and  an  inte¬ 
grated  electroabsorption  modulator-laser  (EML)  using  the 
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transfer-matrix  method  [18],  [19].  Multisection  DFB  lasers 
have  been  modeled  using  a  transmission-line  model  as  well 

The  wavelength-tunable  laser  studied  here,  a  possible 
candidate  for  WDM  systems,  consists  of  a  distributed-  ‘ 
feedback  (DFB)  laser  with  curved  waveguide,  integrated 
with  an  electroabsorption  modulator  which  acts  as  a 
wavelength-tuning  section.  The  structure  for  this  device 
is  identical  to  that  of  an  integrated  electroabsorption  mod¬ 
ulator  and  laser  (EML)  with  the  exception  that  the  facet 
coatings  on  the  laser  and  modulator  are  opposite  to  those 
of  the  conventional  EML.  A  voltage  bias  applied  to  the 
modulator  effectively  controls  the  leising  wavelength  of  the 
DFB  laser.  Previously  reported  results  for  this  device  have 
shown  a  wavelength  tuning  range  of  3.5  nm  and  high-speed 
optical  packet  switching  over  four  wavelength  channels  at 
a  rate  of  2.5  Gb/s  [11]. 

The  purpose  of  this  work  is  to  model  this  wavelength- 
tunable  multi-section  device,  including  spatial  hole  burning 
(SHB)  effects  self-consistently,  and  to  explain  the  charac¬ 
teristics  of  the  spectrum  and  the  wavelength-tuning  be¬ 
havior  of  this  device  by  showing  good  agreement  between 
theory  and  experiment.  We  first  describe  the  device  struc¬ 
ture  in  Section  II.  In  Section  III,  we  explain  our  model 
based  on  the  transfer-matrix  method  that  includes  spatial 
hole-burning  (SHB)  effects.  In  Section  IV,  we  compare 
measured  spectrum  with  theory  and  explain  the  spectral 
features  of  this  wavelength-tunable  laser.  This  includes  ex¬ 
planations  for  the  wavelength-tuning  behavior  of  this  de¬ 
vice  as  a  function  of  modulator  voltage  and  theoretical 
results  which  show  SHB  suppression  due  to  the  use  of  a 
curved  waveguide.  In  Section  V,  we  summarize  our  work 
and  present  our  conclusions. 

II.  Device  structure 

A  schematic  of  the  wavelength-tunable  laser,  which  con¬ 
sists  of  a  DFB  laser  with  curved  waveguide  integrated  with 
an  electroabsorption  modulator  separated  by  an  isolation 
region,  is  shown  in  Fig.  1.  The  length  of  the  DFB  section 
is  300  fim  long,  the  isolation  region  is  about  80  /im,  and 
the  modulator  is  250  ^m.  The  facet  reflectivity  at  the  HR 
coated  facet  on  the  modulator  side  is  997o,  and  the  AR 
coated  facet  reflectivity  on  the  DFB  laser  side  is  estimated 
to  be  1%  using  a  50  window  structure  on  the  AR  side  of 
the  DFB  section.  The  shape  of  this  curved  waveguide  ha.s 
been  modeled  previously  as  a  raised  sine  shape  [11].  How¬ 
ever,  this  waveguide  can  also  be  modelled  as  an  S-bencl 
shape  which  has  more  flexible  parameters  and  provides  a 
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more  general  waveguide  shape.  The  shape  of  the  S-bend 
curved  waveguide  is  given  by  the  following  equations  [16]: 

y{z)  =  2W  {otO<z<Lx 

!,(z)  =  + 

for  Li  z  K  Li  2Lb 

y{z)  =  0  for  Li^2Lb  <z<L  (1) 

where  L  is  the  total  length,  Li  is  the  length  of  the  straight 
section  before  the  S-bend,  Lb  is  equal  to  half  of  the  length 
of  the  S-bend  region,  and  W  is  the  width  of  half  of  the  S- 
bend.  For  this  device,  we  estimated  the  S-bend  parameters 
to  be  L  =  360  /xm,  Li  =  47  ptm,  Lb  =  133  /xm,  and  W  = 
10  /xm. 

In  both  devices,  the  active  region  in  the  DFB  and  the 
modulator  sections  consist  of  7  layers  of  compressively 
strained  InGaAsP  quantum  wells  and  were  fabricated  us¬ 
ing  a  selective  area  MOVPE  technique.  The  DFB  region 
has  an  estimated  compressive  strain  of  at  least  —1%  while 
the  compressive  strain  in  the  modulator  was  measured  to 
be  about  0.5%.  The  peak  wavelength  of  the  photolumi¬ 
nescence  was  measured  at  1.4891  /xm.  The  isolation  region 
between  the  DFB  and  modulator  sections  was  created  by 
shallow  trenching  a  depth  of  about  0.7  /xm  to  yield  an  elec¬ 
trical  isolation  with  a  typical  value  greater  than  50  kfi. 

III.  Theory 

A.  Longitudinal  Fields  in  a  Multi-section  Device 

We  model  the  longitudinal  properties  of  the  optical  fields 
in  the  curved  waveguide  DFB  laser  with  integrated  modula¬ 
tor  and  isolation  sections  using  the  transfer-matrix  method. 
Most  of  the  development  of  the  transfer-matrix  method 
used  in  our  model  can  be  found  in  Ref. [18],  [19]. 

The  laser,  isolation,  and  modulator  sections  can  each 

=zDFB  z=:isol 

be  represented  by  the  transfer  matrices  F  i  F  ,  and 

F  ,  respectively,  which  gives  the  slowly  varying  am¬ 
plitudes  of  the  forward  and  backward  propagating  waves, 
A  (2)  and  B  (2),  respectively,  in  each  section.  The  elements 
of  these  matrices  are  given  by  [18],  [19] 


Fiiil) 

=  r4r- 

i  —  m 

(2) 

Fx,{l) 

1  —  FpTp, 

(3) 

Fn{l) 

=  ’’p  (e’9'  e-‘«') 

1  — 

(4) 

F22{1) 

1  -  rpVm 

(5) 

q  =  -f  Kah^ba 

(6) 

where  (3  =  27rne  (A) /A  is  the  propagation  constant  where 
Tig  is  the  effective  index  of  the  fundamental  guided  mode, 
=  tt/A  where  A  is  the  grating  pitch,  and  the  detuning 
parameter  =  /?  -  /Jq. 

The  interface  between  sections  which  includes  differences 
in  refractive  index  between  each  section  is  also  represented 

by  a  transfer  matrix  T  and  is  derived  from  the  boundary 
counditions  that  tangential  electrical  and  magnetic  fields 
are  continuous  at  the  interface  at  a  given  z  =  /i  or 
given  by 


^(/+) 


srrt  nX 

=  T  (/) 


Mn  1 

Bin  J  ■ 


(9) 


The  forward  and  backward  propagating  waves  in  the 
whole  device  can  be  modeled  by  multiplying  the  matrices 
for  each  section  and  interface  together  as  given  by 


[  Bit:;’)  ] = [^■'<^"•'>1  [^"'<'’>1 

[r”‘(i.)]  ]  (10) 

where  Ltotai  is  the  length  of  the  entire  device,  Ldfs^  Li^Oi 
and  Lmod  are  the  lengths  of  the  DFB  laser,  isolation,  and 
modulator  sections,  respectively,  and  U  is  the  location  of 
section  i,  where  li  =  LjoFBt  h  =  Lbfb  +  Liso^  as  shown 
in  Fig.'l; 

We  model  the  amplified  spontaneous  emission  spectrum 
for  the  integrated  device  by  including  an  equivalent  sponta¬ 
neous  emission  source  [21]  into  the  transfer  matrix  method. 
We  place  an  an  equivalent  spontaneous-emission  surface 
current  density  at  a  position  z^  in  the  (active)  DFB 
section  of  the  cavity  and  solve  the  boundary  conditions  for 
the  electric  and  magnetic  fields  on  either  side  of  the  current 
sheet.  From  this,  we  obtain  the  relation  between  the  field 
amplitudes  near  the  source  position: 


where  Kah  and  n^a  are  the  coupling  coefficients  used  in  the 
coupled-mode  equations  [19].  Next,  we  match  boundary 
conditions  at  the  facets  and  obtain 
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and  ri  and  r2  are  the  complex  reflectivities  at  the  facets. 
The  matrix  F(Ltctai}  is  the  product  of  all  the  transfer  ma¬ 
trices  in  (10),  and  it  describes  the  total  three-section  device 

whereas  the  matrix  F(Ltotai  -  contains  only  those  ma¬ 
trices  between  and  Ltotai- 

The  amplified  spontaneous  emission  intensity  from  each 
facet  is  obtained  by  summing  the  square  of  the  electric  field 
envelopes  jB(0)  and  A[L)  due  to  a  single  source  at  over 
all  the  sources  J,  within  the  cavity  (0  <  r,  <  L)  as  given 
by 

/left  =  \B{0)fdz,  (15) 

/right  =  {l-\r2f)£\A{L)\-dz,.  (16) 

The  ASE  is  calculated  using  (15)  and  (16)  for  each  wave¬ 
length  to  get  the  entire  ASE  spectrum. 

The  photon  density  profile  is  obtained  from  the  electric 
field  envelopes  using 

5(r)cx:l>l(2)p  +  lB(z)l',  (17) 

Again,  the  total  photon  density  is  found  by  summing  the 
photon  density  S{z)  due  to  a  single  source  at  z,  over  all 
the  sources  Jg  within  the  cavity. 

B,  Spatial  Hole  Burning 

In  the  first  iteration  of  the  photon  density  calculation, 
the  gain  and  carrier  density  profiles  along  the  device  are 
considered  uniform.  If  the  applied  injection  current  is 
above  threshold,  the  calculated  photon  density  profile  will 
be  perturbed  as  a  function  of  longitudinal  cavity  position 
and  therefore  not  constant.  The  carrier  density  and  gain 
are  recalculated  using  the  carrier  density  rate  equation  at 
steady  state  given  by 

^  =  ^-AN-BN^-CN^-Vgg{N)S  (18) 
at  eV 

where  1  is  the  injected  current,  e  is  the  electron  unit  charge, 
V  is  the  volume  of  the  laser  active  region,  A  is  the  mono- 
carrier  lifetime  coefficient,  B  is  the  bimolecular  carrier  re¬ 
combination  coefficient,  C  is  the  Auger  carrier  recombina¬ 
tion  coefficient,  Vg  is  the  group  velocity,  g{N)  is  the  ma¬ 
terial  gain  as  a  function  of  carrier  density  A^  and  S  is 
the  photon  density.  In  the  model,  we  consider  a  linearized 
material  gain  given  by  ^^(A^)  =  ^^(A^  -  No)y  where  g'  is 
the  differential  gain  and  No  is  the  carrier  density  at  trans¬ 
parency.  This  results  in  a  nonuniform  carrier  density  and 
profile.  The  refractive  index  profile  also  becomes  nonuni¬ 
form  due  to  the  perturbed  carrier  density  profile  through 
the  relation: 

Ane  =  -^^AA'.  (19) 

47r 

where  Ane  is  the  change  in  effective  index,  A  is  wavelength, 
Of  is  the  linewidth  enhancement  factor,  and  gf  is  differ¬ 
ential  gain.  This  change  in  refractive  index  to  a  nonuni¬ 
form  carrier  density  profile  changes  the  threshold  gain  of 


potential  lasing  modes  and  accounts  for  the  decrease  in 
side  mode  suppression  ratio  and  mode-hopping  associated 
with  spatial  hole  burning.  The  photon  density  and  car¬ 
rier  density  are  continually  recalculated  using  (18)  until 
a  self-consistent  solution  is  found.  Only  a  single  longitu¬ 
dinal  mode  is  taken  into  account.  Our  model  has  been 
applied  to  the  study  of  spatial  hole  burning  in  conven¬ 
tional  DFB  lasers  with  excellent  agreement  with  experi¬ 
mental  data  [19].  It  should  be  noted  that  each  section  of 
the  device  is  futher  subdivided  into  multiple  sections,  and 
the  coupled-mode  equations  are  used  to  derive  the  transfer 
matrices  to  account  for  the  spatially-dependent  refractive 
index  and  gain  profiles  and  the  grating  pitch  parameter. 

C.  Curved  Waveguide  and  Dispersion 

The  curved  waveguide  in  the  DFB  laser  section  induces 
an  effective  variation  in  grating  pitch  Ae//  as  a  function  of 
longitudinal  position  and  is  included  by  [11]: 


A  - 

“  cos4.{z) 

(20) 

(21) 

where  y  and  z  are  defined  in  (1),  and  Ao  is  the  grating 
pitch  for  the  straight  waveguide  case.  Aejj  is  taken  into 
account  in  the  model  through  the  propagation  constant 
Po  =  whereas  in  our  previous  work  [19],  fdo  = 

tt/Ao.  In  addition,  the  use  of  a  curved  waveguide  slightly 
increases  the  effective  length  of  the  waveguide  in  the  DFB 
laser  section  through  which  the  fields  propagate.  This  is 
taken  into  account  in  the  model. 

The  effects  of  dispersion  are  also  considered  in  the  model. 
Over  a  narrow  range  of  A  near  a  reference  wavelength  Ao, 
the  effective  index  (A)  can  be  approximated  by  a  linear 
function  of  A  : 

ne(A)  =  n,(A,)-s(A-A,)  (22) 

where  the  slope  s  is 

S  =  (^g  -  (-^o))  (23) 

Ac 

and  Ug  is  the  group  index,  which  can  be  approximated  as 
a  constant  within  the  wavelength  range  of  interest.  We 
estimate  the  group  refractive  index  from  the  mode  spac¬ 
ing  outside  the  stop  band  in  the  experimental  spectrum, 
and  the  value  of  s  is  estimated  using  previous  experimen¬ 
tal  results  [21]  for  the  InGaAsP  system.  We  then  calculate 
Uo  (Ao)  and  solve  for  ne(A)  using  (22)  over  the  wavelengths 
of  interest. 

D.  Voltage-dependent  Refractive  Index  and  Gain  in  the 
Electoabsorption  Modulator 

Typically  in  EMLs,  the  bandgap  of  the  modulator  is 
slightly  larger  than  that  in  the  DFB  laser  section  with  a 
detuning  of  about  30-60  nm  [8],  [9]  so  that  for  zero  or  for¬ 
ward  modulator  voltage  bias,  there  is  no  attenuation  ex¬ 
cept  for  insertion  loss.  For  a  positive  voltage  bias,  carriers 
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are  injected  so  light  propagating  through  the  modulator 
experiences  gain,  and  the  refractive  index  of  the  modu¬ 
lator  decreases.  If  a  negative  voltage  bias  is  applied  to 
the  modulator,  absorption  is  increased  due  to  quantum- 
confined  Stark  effects  or  Franz-Keldysh  effects  which  in 
turn  increases  the  refractive  index  of  the  modulator.  We 
model  the  voltage-dependent  absorption  and  refractive  in¬ 
dex  in  the  modulator  through  a  simple  linear  additive  term 
given  by 

9mcd  ~  5^0  d"  "^Y  (2^) 


physics  of  the  tuning  mechansims  below.  Since  the  iso¬ 
lation  region  and  modulator  are  homogeneous,  meaning 
that  there  is  no  distributed  feedback  present,  the  round 
trip  through  these  sections  can  also  be  perceived  as  an  ef¬ 
fective  reflectivity  at  the  DFB-isolation  region  interface. 
Varying  the  modulator  voltage  will  change  the  magnitude 
and  phase  of  the  reflectivity  at  this  interface,  hence,  the 
modulator  acts  eis  a  phase-shifting  section.  The  voltage- 
dependent  magnitude  and  phase  of  the  interface  reflectivity 
is  approximately 

r,„t  =  (27) 


Tlmod  ~  ^0 


dn 

dVmod 


(25) 


where  gmoi  and  n^od  are  the  absorption  coeffcient  (1/cm) 
and  effective  index  in  the  modulator,  respectively,  and 
Tio  are  the  gain  and  effective  index  in  the  modulator  with 
no  voltage  bias,  dgIdVmod  and  dnfdVmod  are  empirical  con¬ 
stants,  and  Vmed  is  the  voltage  bias  at  the  modulator.  Note 
that  both  dgldVmod  and  dn/dVmod  are  positive.  Because 
the  change  in  refractive  index  in  the  forward  bias  case  is 
caused  by  carrier  injection  while  the  change  in  refractive 
index  in  the  reverse  bias  case  is  caused  by  the  electroab¬ 
sorption  effect,  the  value  of  dn/dVmod  is  different  for  each 
case.  In  both  cases,  the  change  in  refractive  index  in  the 
modulator  due  to  the  voltage  bias  is  the  mechanism  re¬ 
sponsible  for  the  phase-shifting  effect  as  will  be  explained 
in  Section  IV. 


E.  Leakage  Current  from  the  Modulator 

The  function  of  the  fabricated  isolation  region  is  to  pro¬ 
vide  electrical  isolation  between  the  modulator  and  laser 
sections  so  that  the  biasing  of  one  section  does  not  affect 
the  performance  of  the  other.  However,  the  isolation  region 
in  this  particular  device  does  not  provide  perfect  electrical 
isolation.  When  the  modulator  is  forward-biased,  carriers 
are  injected  into  the  modulator  and  a  small  leakage  current 
flows  from  the  modulator  into  the  laser  section.  This  leak¬ 
age  current  will  decrease  the  refractive  index  of  the  laser 
section  and  is  included  in  the  model  using  a  simple  empir¬ 
ical  relation. 


riDFB  —  rio—  ^mod  (26) 

(*ymod 

where  n^FB  is  the  modified  effective  index  of  the  laser  sec¬ 
tion,  no  is  the  original  effective  index  and  dnieahM^mod 
is  the  change  in  effective  index  due  to  the  leakage  current 
and  is  determined  empirically  by  the  observed  shift  in  las¬ 
ing  wavelength  in  the  measured  spectrum  as  a  function  of 

1  moii  • 

IV.  Results  and  Discussions 

A.  Wavelength-Tuning  Mechansims:  Phase-Shifting  and 
Current  Tuning 

Although  our  transfer-matrix  method  provides  the  full 
numerical  solution  taking  into  account  all  multiple  reflec¬ 
tion,  we  present  a  simple  physical  model  to  illustrate  the 


where  r,nt  and  are  the  magnitudes  of  the  field  reflectiv¬ 
ity  at  the  laser-isolation  region  and  HR  facet,  respectively, 
Liso  and  Lmod  are  the  lengths  of  the  isolation  region  and 
modulator,  respectively,  /?,ao  is  the  propagation  constant 
of  the  isolation  region,  and  ^mod  is  the  voltage-dependent 
propagation  constant  of  the  modulator  given  by 


^  27r  .^mod  fOQ\ 

Pmod  ^  *  2 

where  gmod  and  n^od  are  functions  of  modulator  voltage 
and  given  by  (24)  and  (25)  respectively.  Therefore,  the 
change  in  phase  at  the  interface  due  to  a  voltage-dependent 
change  in  refractive  index  is  given  by 


47rAneLmod 

A 


(29) 


where  An^  is  given  by  the  voltage-dependent  term  in  (25). 

It  is  well  known  that  in  an  HR/AR  coated  DFB  laser, 
the  lasing  wavelength  can  change  dramatically  for  different 
phases  at  the  HR  facet  [22].  By  integrating  a  phase-tuning 
section  into  the  HR-coated  side  of  the  DFB  laser,  we  use 
this  idea  to  achieve  voltage-controlled  wavelength  tuning 
over  the  stopband  of  the  DFB  laser. 

The  other  mechanism  behind  wavelength  tuning  in  this 
integrated  device  is  due  to  leakage  current  from  the  modu¬ 
lator.  An  increasing  positive  modualtor  voltage  will  inject 
carriers  into  the  laser  section,  decrease  the  refractive  index 
of  the  laser,  and  shift  the  spectrum  toward  longer  wave¬ 
lengths.  For  larger  negative  voltage  biases,  carriers  are 
absorbed  from  the  laser  and  the  spectrum  shifts  toward 
shorter  wavelengths.  Although  leakage  current  is  not  usu¬ 
ally  viewed  favorably,  this  effect,  when  combined  with  the 
phase-shifting  effect,  actually  allows  this  device  to  achieve 
a  larger  wavelength-tuning  range  than  phase-shifting  alone. 

The  wavelength-tuning  behavior  of  this  device  naturally 
leads  to  a  discussion  on  the  advantages  and  disadvantages 
of  current  tuning  and  phase-shifting  as  wavelength-tuning 
mechanisms  in  lasers.  The  idea  of  using  different  injection 
currents  and  operating  temperatures  to  produce  a  change 
in  the  wavelength  of  the  lasing  mode  has  long  been  known 
and  has  been  extended  to  a  variety  of  other  types  of  devices 
including  multi-section  lasers.  The  advantage  of  current 
tuning  is  that  one  can  achieve  a  continuous  tuning  range. 
The  main  disadvantage  is  that  for  wavelength  switching, 
as  used  in  traffic  routing  applications,  current  timing  is 
a  relatively  slow  process,  limited  by  the  carrier  lifetime, 
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and  cannot  support  wavelength  switching  at  rates  in  the 
multi-GHz  range.  In  order  to  increase  wavelength  switch¬ 
ing  rates,  several  groups  have  taken  advantage  of  the  faster 
Franz-Keldysh  effects  in  a  reverse-biased  electroabsoprtion 
modulator  to  produce  wavelength  switching.  Nakajima  et 
al.  demonstrated  wavelength  switching  at  10  Gb/s  over  a 
0.2  nm  range  [10].  This  wavelength  tuning  was  acheived  by 
absorption  of  carriers  by  the  modulator  and  is  essentially 
current  tuning  using  a  faster  process. 

Using  the  same  type  of  device  presented  here,  it  has 
been  shown  that  wavelength  switching  at  2.5  Gb/s  between 
4  wavelength  channels  spanning  a  2.1  nm  range  can  be 
achieved[ll].  The  wider  wavelength  range  achieved  here  is 
due  the  fact  that  both  phase-shifting  and  current  tuning 
from  quantum-confined  Stark  effects  were  used  to  produce 
the  wavelength  tuning. 

B,  Spectrum  of  the  Integrated  DFB  Laser- Modulator 

Our  model  was  first  used  to  explain  some  of  the  features 
present  in  the  output  spectrum  of  this  curved  waveguide 
DFB  laser  with  integrated  modulator.  The  spectrum  of 
a  conventional  HR/AR  coated  DFB  laser  was  measured 
above  threshold  as  shown  in  Fig.  2(a).  The  DFB  laser 
spectrum  typically  has  a  lasing  mode  with  a  side  mode 
suppression  ratio  (SMSR)  greater  than  30dB  on  one  side 
of  the  stopband  and  side  modes  with  amplitudes  which  de¬ 
crease  gradually  for  shorter  and  longer  wavelengths  away 
from  the  stopband.  The  measured  spectrum  of  the  inte¬ 
grated  device  with  no  modulator  voltage  bias  is  shown  in 
Fig.  2(b)  and  shows  several  interesting  features  which  differ 
from  the  conventional  DFB  laser.  The  first  feature  is  the 
absence  of  a  visible  stopband.  The  second  feature  is  that 
the  modes  alternate  in  amplitude,  whereas  the  mode  am¬ 
plitudes  of  the  conventional  DFB  laser  gradually  decrease 
away  from  the  stopband.  The  third  feature  is  that  compet¬ 
ing  side  modes  can  be  seen  on  both  sides  of  the  lasing  mode, 
whereas  in  the  conventional  DFB  laser,  there  is  usually  only 
one  competing  mode  located  on  the  other  side  of  the  stop- 
band.  The  theoretically  calculated  spectrum  shown  in  Fig. 
2(c)  exhibits  features  which  are  very  similar  to  the  mea¬ 
sured  spectrum  shown  in  Fig.  2(b).  The  parameters  used 
in  simulations  of  the  integrated  wavelength-tunable  device 
are  shown  in  Table  I. 

The  spectrum  of  the  integrated  device  can  be  explained 
by  looking  at  the  threshold  gain  of  the  modes  as  a  func¬ 
tion  of  wavelength  for  several  cases.  First  we  consider  the 
effect  of  using  a  two-section  DFB  laser-modulator  device 
compared  to  the  single-section  DFB  laser  case.  Shown 
in  Fig.  3(a)  are  the  modes  for  a  single-section  DFB  laser 
with  straight  waveguide  case.  Shown  in  Fig.  3(b)  are 
the  modes  for  the  two-section  integrated  laser-modulator. 
For  the  single-section  DFB  laser,  the  modes  show  a  well- 
defined  stopband.  However,  for  the  two-section  DFB  laser- 
modulator  device,  the  modes  seem  to  separate  into  two 
sets  near  the  stopband,  where  the  set  of  modes  with  higher 
threshold  gain  correspond  to  the  coupled-cavity  effects 
from  the  modulator  region,  and  the  modes  with  lower 
threshold  gain  correspond  to  the  DFB  laser  region.  This 


accounts  for  the  alternating  mode  amplitude  behavior  away 
from  the  stopband  toward  shorter  and  longer  wavelengths 
as  seen  in  the  spectrum. 

In  a  straight  waveguide,  the  grating  pitch  is  constant 
as  a  function  of  longitudinal  cavity  postion  which  produces 
two  well-defined  Bragg  frequency  modes  located  oppositely 
across  the  stopband.  In  the  curved  waveguide,  the  grating 
pitch  is  chirped  as  a  function  of  cavity  position  which  pro¬ 
duces  a  mode-flattening  effect  around  the  stopband  which 
increases  as  a  function  of  S-bend  width  W.  This  effect  is 
illustrated  by  simulating  a  single-section  AR/AR  coated 
DFB  laser  in  order  to  isolate  the  effects  of  using  the  curved 
waveguide  and  more  clearly  illustrate  the  effect.  Figs.  4(a)- 
(c)  where  the  threshold  gain  of  the  theoretically  calculated 
modes  are  plotted  versus  wavlength  for  increeising  values  of 
W.  The  mode-flattening  effect  comes  from  the  fact  that  us¬ 
ing  a  curved  waveguide  or  chirped  grating  pitch,  the  Bragg 
wavelength  is  no  longer  defined  as  a  single  frequency  but  a 
range  of  frequencies.  Therefore,  the  distributed-feedback  is 
weaker  for  one  particular  frequency,  and  the  SMSR  as  de¬ 
termined  by  the  threshold  gain  decreases  as  a  result.  This 
mode-flattening  effect  around  the  stopband  accounts  for 
the  presence  of  multiple  competing  side  modes  in  the  spec¬ 
trum  of  the  integrated  device  as  well  as  the  absence  of  a 
stopband.  Although  a  reduced  SMSR  is  not  desirable,  this 
mode-flattening  effect  may  extend  the  wavelength-tuning 
range  of  this  type  of  device  by  providing  more  competing 
modes  for  lasing,  which  may  be  selected  by  the  modula¬ 
tor  section.  Far  outside  the  stopband,  the  mode  behavior 
was  observed  to  be  determined  by  the  Fabry- Perot  cavity 
in  both  experiment  and  theory  as  expected.  Simulations  of 
the  modes  for  an  HR/AR  coated  DFB  laser  shows  a  similar 
mode-flattening  effect  due  to  the  curved  waveguide. 

C,  Wavelength  Tuning  of  the  Lasing  Mode 

The  model  was  used  to  simulate  the  wavelength-tuning 
behavior  in  the  integrated  DFB  laser-modulator  for  several 
values  of  modulator  voltage.  Fig.  5(a)  shows  the  measured 
spectrum  of  the  integrated  device  for  modulator  voltage 
values  of  0,  +0.47,  +0.78,  +0.87,  +0.92,  and  +1,03  V.  As 
explained  earlier,  the  wavelength  tuning  occurs  due  to  a 
combination  of  phase-shifting  and  current  tuning.  In  gen¬ 
eral,  the  wavelength  tuning  behavior  of  this  device  can  be 
described  in  the  following  way.  As  the  modulator  volt¬ 
age  is  increased  from  0  V,  the  lasing  mode  remains  at  one 
wavelength  until  +0.47  V,  where  the  lasing  mode  switches 
abruptly  to  a  second  wavelength.  As  the  modulator  is  in¬ 
creased  further,  the  lasing  mode  remains  at  this  second 
wavelength  until  +0.78  V,  where  the  lasing  mode  switches 
to  another  wavelength,  and  this  behavior  continues  un¬ 
til  +1.03  V.  Fig.  5(b)  shows  the  theoretically  calculated 
spectrum  for  modulator  values  of  0,  +0.65,  +0.75,  +0.85, 
+0.95,  and  +1.10  Y,  which  shows  reasonably  good  agree¬ 
ment  using  the  linear  model  for  voltage-dependent  refrac¬ 
tive  index. 

For  a  positive  voltage  bias,  the  modulator  gain  and  re¬ 
fractive  index  could  be  modeled  more  accurately  by  making 
them  functions  of  carrier  density  in  the  modulator  which 
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is  a  function  of  injected  current  in  the  modulator.  Also, 
wavelength-tuning  of  the  lasing  mode  for  positive  bias  volt¬ 
ages  is  also  associated  with  changes  in  threshold  current, 
linewidth,  and  output  power.  This  was  not  included  here, 
but  we  refer  to  Ref. [24]  for  a  more  detailed  treatment. 

Wavelength-tuning  of  the  lasing  mode  was  also  observed 
for  negative  modulator  bias  voltages,  although  for  a  more 
limited  range  for  this  particular  device.  While  in  the  pos¬ 
itive  voltage  bias  case,  the  refractive  index  decreases  due 
to  the  injection  of  carriers,  for  this  case,  the  refractive  in¬ 
dex  increases  due  to  absorption  of  carriers  from  quantum- 
confined  Stark  effects.  Fig.  6(a)  shows  the  measured  spec¬ 
trum  for  modulator  voltage  values  of  0,  -0.45,  and  -0.78 
V.  Fig.  6(b)  shows  the  theoretically  calculated  spectrum 
for  modulator  values  of  0,  -0.50,  and  -0.90  V,  which  shows 
good  agreement  with  the  measured  spectrum.  The  lasing 
wavelength  as  a  function  of  positive  and  negative  modula¬ 
tor  voltages  is  shown  for  both  the  experimentally  measured 
and  theoretically  calculated  spectrum  in  Fig.  7.  The  closed 
circles  are  measured  data  points.  The  lasing  wavelengths 
between  data  points,  shown  by  the  lines,  were  observed  but 
not  measured  and  are  accurate  to  within  ±1  angstrom.  The 
open  triangles  are  the  theoretically  calculated  values  and 
show  reasonably  good  agreement  with  the  measured  spec¬ 
trum  values. 

D,  Reduced  Spatial  Hole  Burning  due  to  the  Curved  Waveg- 
uide 

Spatial  hole  burriing  (SHE)  is  an  important  phenomena 
which  has  been  studied  extensively  [25]-[32],  particularly 
in  DFB  lasers,  because  it  can  severely  limit  the  maximum 
output  power  and  cause  mode-hopping  or  multimode  be¬ 
havior.  SHE  is  caused  by  a  nonuniform  distribution  of  the 
photon  density  along  the  longitudinal  cavity  which  can  be 
caused  by  the  distributed  feedback  in  DFE  lasers  as  well  as 
asymmetries  in  HR/AR  coated  DFE  lasers.  A  greater  pho¬ 
ton  density  in  one  region  of  a  DFB  laser  will  cause  greater 
carrier  depletion  in  that  region  so  the  carrier  density  pro¬ 
file  also  becomes  nonuniform  along  the  cavity.  This  in  turn 
produces  a  nonuniform  refractive  index  profile  which  affects 
the  threshold  gain  conditions  of  the  modes  and  can  result  in 
mode-hopping  or  multimode  lasing  at  moderate  injection 
currents  above  threshold.  The  model  was  used  to  calcu¬ 
late  the  longtudinal  photon  density  profile  for  the  straight 
and  curved  waveguide  cases  as  shown  in  Fig.  8.  The  cor¬ 
responding  facet  coatings  and  device  sections  are  labelled 
with  cavity  position.  The  solid  lines  are  for  Vmod-^^  and 
dashed  lines  are  for  V;nod=H-0.65  V.  Considering  only  the 
laser  section,  the  photon  density  profile  show^s  a  larger  per¬ 
turbation  in  the  straight  waveguide  case  than  the  curved 
waveguide  case.  One  way  to  evaluate  the  amount  of  spatial 
hole  burning  present  in  a  DFB  laser  is  to  calculate  a  flat¬ 
ness  parameter  [33].  We  use  a  modified  flatness  parameter 
given  by 

F=jJ{IM„-h„)‘dz  (30) 


where  L  is  the  length  of  the  DFB  laser  section,  /  (^)norm 
the  normalized  local  field  intensity,  and  lavg  is  the  average 
normalized  intensity.  The  integration  is  performed  only 
over  the  DFB  laser  cavity  length.  Note  that  a  greater  F 
value  indicates  larger  SHB  effects.  For  the  Vmod  ==  0  V  case, 
the  calculated  value  of  F  was  0.78  for  the  straight  waveg¬ 
uide  and  0.21  for  the  curved  waveguide.  Very  similar  results 
were  obtained  for  the  other  values  of  Vmod-  These  results 
suggest  that  the  curved  waveguide  exhibits  less  SHB  than 
the  straight  waveguide  case.  This  is  significant  because  a 
device  with  less  SHB  can  remain  at  stable  operation  at 
higher  output  powers. 

The  photon  density  profile  of  the  curved  waveguide  as 
shown  in  Fig.  8  exhibits  an  oscillatory  behavior  [34],  [35]. 
This  can  be  explained  by  the  fact  that  this  curved  waveg¬ 
uide  produces  multiple  distributed  phase  shifts  due  to  the 
chirped  pitch  grating.  Other  groups  have  modeled  multi¬ 
ple  phase-shift  devices  which  show  multiple  peaks  in  the 
photon  density  profile  [36]  and  also  distributed  phase-shift 
devices  which  exhibit  a  rounding  and  flattening  of  a  pho¬ 
ton  density  profile  peak  when  compared  to  an  abrupt  phase 
shift  [16]. 

Several  designs  can  be  implemented  to  improve  the  per¬ 
formance  of  this  device. 

1)  The  present  device  mode-shifts  between  only  two 
modes.  By  optimizing  the  curved  waveguide  shape  and 
increasing  the  mode-flattening  effect,  we  should  be  able 
to  obtain  mode-shifting  between  three  or  more  modes, 
thereby  increasing  the  wavelength-tuning  range.  However,- 
the  tradeoff  would  be  a  decrease  in  side-mode  supression 
ratio. 

2)  The  concept  of  an  amplitude-modulated  grating  [37] 
can  be  used  to  extend  the  phase-shifting  range.  Although 
originally  designed  to  flatten  the  photon  profile  and  reduce 
SHB  effects,  this  idea  may  also  be  used  to  enhance  the 
phase-shifting  range. 

V.  Conclusions 

In  conclusion,  we  have  theoretically  modeled  a  DFB 
laser  with  a  curved  waveguide  integrated  with  an  electroab¬ 
sorption  modulator  which  acts  as  a  wavelength-tuning  sec¬ 
tion.  The  characteristics  of  the  output  spectrum  have  been 
identified,  and  the  mechanisms  behind  the  wavelength  tun¬ 
ing,  namely,  phase-shifting  and  current  tuning  have  been 
explained  using  simulation  results  which  show  good  agree¬ 
ment  with  experiment.  The  model  also  shows  that  the 
curved  waveguide  suppresses  the  photon  density  profile  and 
reduces  SHB  compared  to  the  straight  waveguide  case.  The 
development  of  wavelength-tunable  lasers  with  greater  tun¬ 
ing  range  and  higher  power  operation  should  be  achievable 
using  a  curved  waveguide. 

Acknowledgments 

This  work  was  supported  by  ONR  under  Grant  N00014- 
96-1-0303  and  N00014-96-1-0902.  The  authors  would  like 
to  acknowledge  R.  People,  C.  Bethea,  A.  I\l.  Sergent,  and 
W.  T.  Tsang  for  their  technical  expertise  and  discussions. 


IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS 


7 


References 

{l]  K.  Kobayashi  and  I.  Mito,  “Single  frequency  and  tunable  laser 
diodes  "  IEEE  J.  Lightwave  Tech.,  vol.  6,  pp.  1623-1633,  1988. 

[2]  H.  Kawaguchi,  “Progress  in  optical  functional  devices  using  two- 
section  laser  diodes/amplifiers,”  lEE  ProceedingS’J,  vol,  140,  pp. 
3-15,  1993. 

[3]  Y,  Kawamura,  K.  Wakita,  Y.  Yoshikuni,  Y.  Itaya,  and  H.  Asahi, 
“Monolithic  integration  of  a  DFB  laser  and  an  MQW  optical 
modulator  in  the  1.5  range,”  IEEE  J.  Quantum  Electron,, 
vol.  23,  pp.  915-918,  1987. 

[4]  M.  Suzuki,  Y.  Noda,  H,  Tanaka,  S.  Akiba,  Y.  Kushiro,  and  H. 
Issikhi,  “Monolithic  integration  of  an  InGaAsP/InP  distributed 
feedback  laser  and  electroabsorption  modulator  and  vapor  phase 
epitaxy,”  IEEE  J,  Quantum  Electron,,  vol.  23,  pp.  915-918, 
1987. 

[5]  M.  Aoki,  M.  Suzuki,  H.  Sano,  T.  Kawano,  T.  Ido,  T.  Tani- 
watari,  K.  Uomi,  and  A.  Takai,  “InGaAs/InGaAsP  MQW  elec¬ 
troabsorption  modulator  integrated  with  a  DFB  Iztser  fabricated 
by  band-gap  energy  control  selective  area  MOCVD,”  IEEE  J, 
Quantum  Electron,,  vol.  29,  pp.  2088-2096,  1993. 

[6]  T.  Kato,  T.  Sasaki,  K.  Komatsu,  and  I.  Mito,  “DFB- 
LD/modulator  integrated  light  source  by  bandgap  energy  con¬ 
trolled  selective  MOVPE,”  Electron.  Lett.,  vol.  28,  pp.  153-154, 
1992. 

[7]  T.  Tanbun-Ek,  Y.  K.  Chen,  J.  A.  Grenko,  E.  K.  Byrne,  J.  E. 
Johnson,  R.  A.  Logan,  A,  Tate,  A.  M.  Sergent,  K.  W.  Wecht,  P. 
F.  Scriortino,  Jr.,  S.  N.  G.  Chu,  “Integrated  DFB-DBR  laser 
modulator  grown  by  selective  area  metalorganic  vapor  phase 
epitaxy  growth  technique,”  J.  Crystal  Growth,  vol.  145,  pp, 
902-906,  1994. 

[8]  A.  Ramdane,  F.  Devaux,  N,  Souli,  D.  Delprat,  and  A.  Ougaz- 
zaden,  “Monolithic  integration  of  multiple-quantum-well  lasers 
and  modulators  for  high-speed  transmission,”  IEEE  J.  Selected 
Topics  Quantum  Electron.,  vol.  2,  pp.  326-335, 1996. 

[9]  H.  Takeuchi,  K.  Tsuzuki,  K.  Sato,  M,  Yamamoto,  Y.  Itaya,  A. 
Sano,  M.  Yoneyama,  and  T.  Otsuji,  “Very  high-speed  light- 
source  module  up  to  40  Gb/s  containing  an  MQW  electroabsorp¬ 
tion  modulator  integrated  with  a  DFB  laser,”  IEEE  J.  Selected 
Topics  Quantum  Electron.,  vol.  3,  pp.  336-343,  1997, 

IlO]  H.  Nakajima,  J.  Charil,  S.  Slempkes,  A.  Gloukhian,  B.  Pierre,  J. 
Landreau,  and  Y,  Raffle,  “Very  high-speed  wavelength  switching 
capability  Franz-Keldysh  electroabsorption  DFB  lasers,”  Optical 
Fiber  Conj.  Tech.  Digest,  p.  276,  1996. 

[11]  T.  Tanbun-Ek,  L.  Adams,  G.  Nykolak,  C.  Bethea,  R.  People,  A. 
Sergent,  P.  Wisk,  P.  Sciortino,  Jr.,  S.  Chu,  T.  Fullowan,  and  W, 
T.  Tsang,  “Broad-band  tunable  electroabsorption  modulated 
laser  for  WDM  application,”  IEEE  J.  Selected  Topics  Quantum 
Electron.,  vol.  3,  pp.  960-966, 1997. 

[12]  A.  Katzir,  A.  C.  Livanos,  and  A.  Yariv,  “Chirped-grating  output 
couplers  in  dielectric  waveguides,”  Appl.  Phys.  Lett.,  vol.  30, 
pp.  225-226,  1977. 

[13]  P.  Zhou  and  G.  S.  Lee,  “Chirped  grating  A/4  shifted  DFB  laser 
with  uniform  longitudinal  field  distribution,”  Electron.  Lett., 
vol.  26,  pp.  16660-1661,  1990. 

[14]  M.  Okai,  T.  Tsuchiya,  K.  Uomi,  N.  Chinone,  and  T.  Harada, 
“Corrugation-pitch  modulated  MQW-DFB  lasers  with  narrow 
spectral  linewidth,”  IEEE  J.  Quantum  Electron.,  vol.  27,  pp. 
1767-1772,  1991. 

[15]  W.  T.  Tsang,  R.  M.  Kapre,  R.  A.  Logan,  and  T.  Tanbun-Ek, 
“Control  of  lasing  wavelength  in  distributed  feedback  lasers  by 
angling  the  active  stripe  with  respect  to  the  grating,”  IEEE 
Photon.  Technol.  Lett.,  vol.  5,  pp.  10-12,  1993. 

[16]  H.  Hillmer,  A.  Grabmaier,  S.  Hansmann,  H.-L.  Zhu,  H. 
Burkhard,  and  K.  Magari,  “Tailored  DFB  lasers  properties  by 
individually  chirped  gratings  using  bent  waveguide,”  IEEE  J. 
Selected  Topics  Quantum  Electron.,  vol.  1,  pp.  356-362,  1995. 

[17]  J.  Salzman,  H,  Olesen,  A.  Moller-Larson,  O.  Albrektsen,  J.  Han- 
berg,  J.  Norregaard,  B.  Jonsson,  and  B.  Tromborg,  “Distributed 
feedback  lasers  with  an  S-bent  waveguide  for  high-power  single¬ 
mode  operation,”  IEEE  J.  Selected  Topics  Quantum  Electron., 
vol.  1,  pp.  346-355,  1995. 

[is]  W.  Fang,  S.  L.  Chuang,  T.  Tanbun-Ek,  and  Y.  K.  Chen,  “Mod¬ 
eling  and  experiment  of  1.55  /im  integrated  electroabsorption 
modulator  with  distributed-feedback  laser,”  Proc.  SPIE;  Opto¬ 
electronic  Integrated  Circiiifs,  Eds.  Y.  S.  Park  and  R.  V.  Ra- 
maswamy,  vol.  3006,  pp.  207-215,1997. 

[19]  W.  Fang,  A.  Hsu,  S.  L.  Chuang.  T.  Tanbun-Ek,  and  A.  M. 
Sergent,  “Measurement  and  modeling  of  distributed-feedback 


Lasers  with  spatial  hole  burning,”  IEEE  J.  Selected  Topics 
Quantum  Electron.,  vol.  3,  pp.  547-554,  1997. 

[20]  C.  H.  Chen  and  G.  GrifTel,  “Static,  dynamic,  and  noise  analysis 
of  multisection  DFB  lasers  using  frequency-domain  transmission 
line  model,”  IEEE  J.  Quantum  Electron.,  vol.  34,  pp.  1533-1544, 
1998. 

[21]  C.  S.  Chang,  S.  L.  Chuang,  J.  R.  Minch,  W.  Fang,  Y.  K. 
Chen,  and  T.  Tanbun-Ek,  “Amplified  spontaneous  emission 
spectroscopy  in  strained  quantum-well  lasers,”  IEEE  J.  Selected 
Topics  Quantum  Electron.,  vol.  1,  pp.  1100-1107,  1995. 

[22]  S,  Chinn,  “Effects  of  mirror  reflectivity  in  a  distributed-feedback 
laser,”  IEEE  J.  Quantum  Electron.,  vol.  9,  pp.  574-580,  1973. 

[23]  M.  Kitamura,  M.  Yamaguchi,  K.  Emura,  I.  Mito,  and  K. 
Kobayashi,  “Lasing  Mode  and  Spectral  Linewidth  Control  by 
Phase  Tunable  Distributed  Feedback  Laser  Diodes  with  Double 
Channel  Planar  Buried  Hcterostructure  (DFB-DC-PBH  LD’s),” 
IEEE  J.  Quantum  Electron.,  vol.  21,  pp.  415-417,  1985. 

[24]  H.  Olesen,  X.  Pan,  and  B.  Tromborg,  “Theoretical  Analysis  of 
Tuning  Properties  for  a  Phase-Tunable  DFB  Laser,”  IEEE  J. 
Quantum  Electron.,  vol.  24,  pp.  2367-2375,  1987. 

[25]  M.  Aoki,  K.  Uomi,  T.  Tsuchiya,  S.  Sasaki,  M.  Okai,  and  N.  Chi¬ 
none,  “Quantum  size  effect  on  longitudinal  spatial  hole  burning 
in  MQW  A/4-shifted  DFB  laisers,”  IEEE  J.  Quantum  Electron., 
vol.  27,  pp.  1782-1789,  1991. 

[26]  F.  Girardin,  G.  H.  Duan,  and  A.  Talneau,  “Modeling  and  mea¬ 
surement  of  spatial-hole-buming  applied  to  amplitude  modu¬ 
lated  coupling  distributed  feedback  lasers,”  IEEE  J.  Quantum 
Electron.,  vol.  3,  pp.  834-841,  1995. 

[27]  M.  R.  Phillips,  T.  E.  Darcie,  and  E.  J.  Flynn,  “Experimental 
measure  of  dynamic  spatial-hole  burning  in  DFB  lasers,”  IEEE 
Photon.  Tech.  Lett.,  vol.  4,  pp.  1201-1203,  1992. 

[28]  S.  Hansmann,  H.  Walter,  H.  Hillmer,  and  H.  Burkhard,  “Static 
and  dynamic  properties  of  InGaAsP-InP  distributed  feedback 
lasers  -  A  detailed  comparison  between  experiment  and  theory,” 
IEEE  J.  Quantum  Electron.,  vol.  30,  pp.  2477-2484,  1994. 

[29]  I.  Orfanos,  T.  Sphicopoulos,  A.  Tsigopoulos,  and  C.  Caroubalos, 
“A  tractable  above-threshold  model  for  the  design  of  DFB  and 
phase-shifted  DFB  lasers,”  IEEE  J.  Quantum  Electron.,  vol. 
27,  pp.  946-956,1991. 

[30]  M.  G.  Davis  and  R.  F.  O’Dowd,  “A  transfer  matrix  method 
based  large-signal  dynamic  model  for  multielectrode  DFB 
lasers,”  IEEE  J.  Quantum  Electron.,  vol.  30,  pp.  2458-2466, 
1994. 

[31]  P.  Correc,  “Stability  of  phase-shifted  DFB  lasers  against  hole 
burning,”  IEEE  J.  Quantum  Electron,  vol. ,30,  pp.  2467-2476, 
1994, 

[32]  S.  F.  Yu,  “A  quasi-three-dimensional  large-signal  dynamic  model 
of  distributed  feedback  lasers,”  IEEE  J.  Quantum  Electron.,  vol. 
32,  pp.  424-432,1996. 

[33]  H.  Soda,  Y.  Kotaki,  H.  Sudo,  H.  Ishikawa,  S.  Yamakoshi,  and 
H.  Imai,  “Stability  in  single  longitudinal  mode  operation  in 
GalnAsP/InP  phase  adjusted  DFB  lasers,”  IEEE  J.  Quantum 
Electron.,  vol.  23,  pp.  804-814,  1987. 

[34]  A.  Hsu,  W.  Fang,  S.  L.  Chuang,  T.  Tanbun-Ek,  C.  Bethea,  and 
R.  People,  “Integrated  tunable  laser  with  mode  selection  mod¬ 
ulator,”  Photonics  West,  San  Jose,  CA,  January  24-30,  1998. 

[35]  A.  Hsu,  W.  Fang,  and  S.  L.  Chuang,  “Modeling  of  normal 
and  backward  integrated  electroabsorption  modulator  and  laser” 
Integrated  Photonics  Research,  Victoria,  British  Columbia, 
Canada,  March  29-April  3,  1998. 

[36]  H.  Ghafouri-Shiraz  and  B.  S.  K.  Lo,  Distributed  Feedback  Laser 
Diodes,  New  York,  NY:  Wiley,  1996. 

[37]  G.  Morthier,  K.  David,  P.  Vankwikelberge,  and  R.  Baets,  “A 
new  DFB-Iaser  diode  with  reduced  spatial  hole  burning,”  IEEE 
Photon.  Ttchnoi  Lett.,  vol.  2,  pp.  3SS-390,  1990. 


IEEE  JOURNAL  OF  QUANTUM  ELECTRONICS 


8 


Fig.  1.  Top:  A  side  view  schematic  of  the  two-section  wavelengtli- 
tunable  laser  which  consists  of  a  DFB  laser  section  at  the  AR 
side  and  a  modulator  at  the  HR  side  of  the  device  which  tunes 
the  lasing  wavelength,  Ltotal  is  the  length  of  the  entire  device, 
Ldfb.  Liso,  and  Lmod  are  the  lengths  of  the  DFB  laser,  iso¬ 
lation,  and  modulator  sections,  respectively,  l\  =  and 

h  =  L£)FB  +  Liso-  Bottom:  A  top  view  schematic  illustrating 
the  curved  waveguide  of  the  DFB  laser  section. 


L.  E.  Adams. 


Fig.  2.  (a)  Measured  spectrum  of  a  conventional  HR-AR  coated 

DFB  laser.  The  spectrum  typically  has  a  lasing  mode  with  good 
G.  Nykolak.  SMSR  on  one  side  of  the  stopband  and  side  modes  with  am¬ 

plitudes  which  decrezuse  gradually  for  shorter  and  longer  wave¬ 
lengths  away  from  the  stopband,  (b)  Measured  spectrum  for  the 
integrated  wavelength-tunable  laser.  This  spectrum  differs  from 
the  conventional  DFB  laser  in  that  the  stopband  is  less  visible, 
T.  Tanbun-Ek.  the  modes  alternate  amplitudes  away  from  the  stopband,  and 

competing  side  modes  are  present  on  either  side  of  the  lasing 
mode,  (c)  Theoretically  calculated  spectrum  of  the  integrated 
wavelength-tunable  laser.  This  spectrum  shows  similar  features 
to  the  experiementally  measured  spectrum  in  (b). 


Fig.  3.  (a)  Modes  for  a  one-section  DFB  laser,  where  the  modes  show 
a  well-defined  stopband,  (b)  Modes  for  the  two-section  integrated 
laser-modulator:  The  modes  show  a  splitting  near  the  stopband, 
the  upper  modes  corresponding  to  the  modulator  region  and  the 
lower  modes  corresponding  to  the  DFB  laser  region. 


A.  Hsu. 


S.  L.  Chiiang. 


W.  Fang. 


Fig.  4.  The  theoretically  calculated  threshold  gain  of  the  modes 
for  an  AR/AR  coated  DFB  laser  as  a  function  of  wavelength  for 
(a)  W  =  O^m,  (b)  W  =  5;im,  and  (c)  W  =  15;4m.  In  general, 
increasing  the  width  of  curved  waveguide  W  produces  a  mode¬ 
flattening  effect  around  the  stopband,  which  can  enhance  mode 
tunability. 


Fig.  5.  (a)  Measured  spectrum  of  the  integrated  device  at  modulator 
voltage  values  of  0,  4-0.47,  +0.78,  4-0.87,  4'0.92,  and  4“103.  (b) 
Theoretically  calculated  spectrum  of  the  integrated  devices  at 
modulator  voltage  values  of  0,  +0.65,  +0.75,  +0.85,  +0.95,  and 
+  1.10. 


Fig.  6.  (a)  Measured  spectrum  of  the  integrated  device  at  modulator 
voltage  values  of  0,  -0.45,  and  -0.7S.  (b)  Theoretically  calculated 
spectrum  of  the  integrated  devices  at  modulator  voltage  values 
of  0,  -0.50,  and  -0.90. 
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Fig.  7.  Measured  and  theoretically  calculated  values  of  the  las¬ 
ing  wavelength  as  a  function  of  modulator  voltage.  The  closed 
circles  represent  measured  l2ising  wavelengths.  The  lines  repre¬ 
sent  observed  lasing  wavelengths  and  are  accurate  to  within  ±1 
angstrom.  The  open  triangles  represent  the  theoretically  calcu¬ 
lated  values. 


Fig.  8.  The  theoretically  calculated  photon  density  as  a  function 
of  longitudinal  cavity  position  for  the  straight  waveguide  (W=0 
^m)  and  curved  waveguide  (W=10  tim)  cases.  The  solid  lines 
are  for  Knod==0,  and  dashed  lines  are  for  a  different  lasing  mode 
at  V^ojf=-|-0.65.  The  curved  waveguide  shows  less  perturbation 
in  the  photon  density  profile  and  therefore  less  SHE  effects  com¬ 
pared  to  the  straight  waveguide  case. 
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Table  I 


Parameter  List  Used  for  Simulations 

power  reflectivity  at  HR  facet  Rfm 

0.99 

power  reflectivity  at  AR  facet  R^r 

0.01 

phase  at  HR  facet  <f>HR 

30" 

phase  at  AR  facet  <f>/iR 

0" 

coupling  coefficient-length  product  kL 

3.0 

length  of  laser  section  Ldfb  (/^^) 

360 

length  of  isolation  section  L/so  i^rn) 

80 

length  of  modulator  section  La/oo  1 

260 

grating  pitch  A© 

0.2308 

S-bend  parameter  Li  (//m) 

47 

S-bend  paramter  Lb 

133 

S-bend  parameter  W  (fim) 

10 

group  index  rip 

3.6 

linewidth  enhancement  factor  Og 

4 

optical  confinement  factor  T 

0.03 

transparent  carrier  density  No  (cm~^) 

1.5  X  10-1® 

differential  gain  (cm“^) 

2.7  X  10-1® 

dnuak/dVmod 

0.002169 

dg/dVmod  (cm“^ 

8.0 

dnIdVmcd  for  Vmod  >  0 

0.013135 

dn/dVmod  for  Vmod  <  0  {V~'^) 

0.01129 

group  velocity  Vp  [m/s) 

8.33  X  10® 

dispersion  parameter  s 

0.21930 

(C) 


2  Ca)  ,( ^  ("c:  > 


1542  1544  154$  1548  1550  1552  1554  1556 
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Abstract — We  present  both  theoretical  and  experimental  studies 
of  four-wave  mixing  in  a  distributed-feedback  (DFB)  laser.  The 
measured  amplified  spontaneous  emission  and  probe  gain  spectra 
are  also  used  to  extract  the  DFB  laser  parameters.  We  use 
mean-field  theory  to  investigate  both  neariy  degenerate  and  non- 
d^enerate  four-wave  mixings.  For  small  pump-probe  detunings, 
the  dependence  of  multiple  conjugate  waves  on  the  pump  and 
probe  powers  is  measured.  For  larger  detunings,  the  gain  of 
the  probe  signal  and  the  conversion  efficiency  are  measured 
in  great  detail  up  to  a  detuning  of  ±550  GH^  The  measured 
conversion  efficiencies  as  a  function  of  the  pump-probe  detuning 
and  the  power  scaling  laws  are  compart  with  those  of  the 
theoretical  results.  We  also  demonstrate  the  DFB  cavity  effects 
on  the  conversion  efficiency  spectrum. 

Index  Terms —  Amplified  spontaneous  emission,  distributed- 
feedback  lasers,  optical  frequency  conversion,  optical  mixing, 
wavelength-division  multiplexing. 

I.  iNTOODUCnON 

OUR-WAVE  MIXING  (FWM)  is  one  of  the  top  contend¬ 
ing  technologies  for  wavelength  conversion  applications 
in  future  multiple  wavelength  telecommunication  systems  [1]. 
Wavelength  conversion  using  four-wave  mixing  in  semicon¬ 
ductor  optical  amplifiers  has  been  demonstrated  at  high  bit 
rates  [2]-[5]  and  for  fi^uency  shifts  of  several  THz  [6].  Four- 
wave  mixing  is  one  of  the  few  technologies  which  converts 
amplitude  and  Ifrequency  modulation  equally  well  and  has  the 
ability  to  convert  multiple  signals  in  parallel. 

Most  four-wave  mixing  experiments  for  wavelength  con¬ 
version  reported  in  the  literature  use  a  semiconductor  optical 
amplifier  (SOA)  as  the  mixing  device.  A  strong  external 
pump  signal,  and  a  weaker  external  probe  signal,  are  typically 
injected  into  the  SOA  and  the  new  signal,  called  the  conjugate 
wave,  is  extracted  from  the  device.  A  simpler  method  to  obtain 
four-wave  mixing  is  to  use  a  distributed-feedback  (DFB)  laser 
as  the  mixing  device.  It  is  well  known  that  DFB  lasers  provide 
a  single,  stable  lasing  mode.  This  feature  provides  a  mixing 
medium  with  a  very  intense  pump  beam  already  in  the  cavity. 
All  that  is  needed  is  an  external  probe  signal  to  be  inject^ 
into  the  cavity  to  create  four-wave  mixing.  Thus,  only  two 
laser  devices  are  needed  to  perform  the  experiment. 

Studies  of  four-wave  mixing  in  DFB  lasers  for  small  de¬ 
tunings  have  recently  been  reported  [7],  [8].  Also,  four-wave 
mixing  with  higher  detunings  (greater  than  1  THz)  has  been 
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observed  for  a  conventional  DFB  laser  [9],  although  restricted 
to  a  few  cavity  resonances,  and  a  quarter-wavelength  shifted 
DFB  laser  [10].  Recently,  we  observed  higher-order  mixing 
peaks  for  small  detunings  of  the  probe  beam  from  the  pump 
and  demonstrated  the  dependence  of  the  scaling  law  of  the 
power  of  the  conjugate  waves  on  the  powers  of  the  pump  and 
probe  signals  [1 1].  We  also  studied  the  efficiency  of  the  four- 
wave  mixing  wavelength  conversion  process  up  to  detunings 
of  ±550  GHz. 

In  this  paper,  we  present  new  experimental  data  for  the 
probe  gain  and  conversion  efficiency  spectra  and  compare 
with  the  results  of  a  general  theory  for  four-wave  mixing 
inside  a  cavity  of  a  laser  operating  above  threshold.  The 
probe  gain  and  the  amplified  spontaneous  emission  spectrum 
from  the  DFB  lasers  are  measured  and  compared  with  our 
theoretical  model  with  excellent  agreement.  Our  experimental 
results  for  the  four-wave  mixing  are  also  well  explained  by  the 
theoretical  model.  The  general  FWM  theory  includes  carrier 
density  modulation  (CDM),  spectral-hole  burning  (SHB),  and 
carrier  heating  (CH)  and  can  thus  be  used  for  both  nearly 
degenerate  and  nondegenerate  detuning  cases.  Two  special 
cases  are  derived  to  account  for  our  experimental  conditions. 
For  the  nearly  degenerate  case,  we  focus  on  the  multiple 
conjugate  waves  and  their  dependences  on  both  the  powers  of 
the  pump  and  probe  waves.  The  contributions  from  the  SHB 
and  CH  nonlinearities  are  negligible  in  this  case.  We  keep 
the  general  power  level  of  both  the  pump  and  probe  waves. 
For  the  conversion  efficiency  measurements,  we  consider  a 
probe  wave  power  much  smaller  than  the  pump  power.  An 
undepleted  pump  power  approximation  is  then  used  to  consider 
the  FWM  conversion  efficiency  for  the  whole  range  of  the 
detuning  frequency. 

In  Section  n,  we  first  derive  an  expression  for  the  gain  of 
a  signal  injected  into  a  DFB  cavity  for  comparison  with  the 
probe  gain  spectrum.  We  then  present  a  general  theory  for 
four-wave  mixing  in  a  DFB  laser.  We  detail  our  experimental 
measurements  in  Section  H.  In  Section  HI,  we  discuss  our 
experimental  results  and  compare  with  the  theory.  Finally,  in 
Section  IV,  we  end  with  a  brief  conclusion. 

n.  Theory 

A.  Probe  Gain 

Light  traveling  in  a  DFB  cavity  can  be  described  with  an 
electrical  field  of  the  form 

E{2)  =  yU{x)  [X(2)e*^-*  +  (1) 
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where  U{x)  is  the  transverse  waveguide  mode  profile,  /?o  = 
TTiTT /A  is  the  Bragg  wave  number  where  A  is  the  period  of  the 
grating  and  m  is  the  order  of  the  grating,  and  A{z)  and  B{z) 
are  the  coefficients  representing  the  slowly  varying  amplitudes 
of  the  forward  and  backward  propagating  waves,  respectively. 
These  quantities  obey  the  well-known  coupled-mode  equations 
[12]  whose  solutions  can  be  written  as 


where 


)  =  Aoe'^^ -i- Botme-^^ 

(2) 

)=Aorpe‘«*-l-B<,e-‘’* 

(3) 

«6o 

(4) 

"9  + A/3 

-«o6 

(5) 

9  + a3 

=  i  (A;0)^  -|-  Kab^ba 

(6) 

=  9r  +  *9; 

(7) 

and  A/3  =  ^  where  /3  =  -  i(G/2)  is  the  propaga¬ 

tion  constant  of  the  guided  (fundamental)  optical  mode,  with 
Tieff  as  the  effective  refractive  index,  ko  as  the  wave  number 
in  free-space  and  G  is  the  net  modal  gain.  K^h  and  Kta  are 
the  coupling  coefficients  of  the  DFB  grating.  Ao  and  Bo  are 
constants. 

Assume  that  the  facet  at  z  =  0  has  a  field  reflection 
coefficient  rj  and  the  facet  at  z  =  L  has  a  reflection  coefficient 
r2.  For  an  incident  (probe)  wave  at  z  =  0  with  an  electric  field 
Eia,  the  boundary  conditions  are 


Pump 


Fig.  1.  A  plot  showing  pump,  probe,  and  multiple  conjugate  waves  used  in 
the  four  wave-mixing  theory. 


B.  Four-Wave  Mixing  Theory 

We  first  consider  a  general  model  for  the  four-wave  mixing 
for  both  nearly  degenerate  and  nondegenerate  cases.  The 
three  mechanisms;  carrier-density  modulation,  spectral-hole 
burning,  and  carrier  heating  are  treated  with  equal  footing. 
The  electrical  field  inside  the  mixing  laser  cavity  is  expressed 
as  =  £(t)  exp(-ia>x.t)  where  ui,  is  the  frequency  of 
the  mixing  laser  (pump  beam).  Similarly,  the  electrical  field 
of  the  injection  laser  beam  (probe  beam)  is  written  as  = 
£i„(t)  exp  +  n)t]  where  is  the  detuning  frequency 

between  the  pump  and  probe  beams.  Based  on  the  mean-field 
theory  [7],  [8],  the  rate  equations  for  the  carrier  density  N{t) 
and  field  E{t)  are 


j4(0)  =riJ5(0)  +  tia^in 
r2e2‘^»'^A(I)=B(I)  (9) 

£out=t2iA(I)  (10) 


where  ti2  is  the  field  transmission  coefficient  from  the  air  into 
the  semiconductor  and  <21  is  *6  transmission  coefficient  from 
the  semiconductor  into  the  air.  We  derive  the  expression  for 
the  power  gain  of  the  incident  electric  field  to  be 


*  in 

|£outP 
■  l^inP 
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where  the  coefficients  ci  and  C2  are 
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(15) 


where  J  is  the  injection  current  density,  q  is  the  magnitude  of 
a  unit  electron  charge,  d  is  the  thickness  of  the  active  region, 
Vg  is  the  group  velocity,  h  is  Planck’s  constant,  Tir  is  the 
refractive  index,  and  Tp  are  the  earner  and  photon  lifetimes, 
respectively,  fee  is  the  coupling  constant  which  accounts  for  the 
coupling  of  the  probe  beam  into  the  mixing  cavity,  g  is  the 
material  gain,  fe*  is  the  wave  number  of  the  pump  beam,  and 
fe  is  the  wave  number  of  the  mixing  device.  In  general,  both 
g  and  fe  are  functions  of  the  carrier  density  N  and  the  optical 
field  intensity  S  =  |£p. 

Let  us  first  consider  that  the  probe  beam  is  absent  (i.e., 
Ein  =  0).  By  setting  the  time  derivative  to  zero  and  adding 
a  subscript  “e”  to  the  variables,  we  obtain  the  threshold 
condition 


BVgQo  = 


1 


(16) 


and  the  field  intensity  of  the  pump  beam 


Note  that  the  power  gain  expression  (1 1)  is  derived  without 
the  presence  of  the  pump  wave.  The  calculated  power  gain 
spectrum  of  the  probe  wave  using  (11)  will  be  compared  with 
the  measured  data  to  extract  the  parameters  for  the  mixing 
DFB  laser. 


where  the  threshold  current  density  is  Jth  =  qdNelTei^e)-  ^^ 
we  include  die  probe  beam,  we  can  express  the  carrier  density 
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N,  the  optical  field  E,  and  the  optical  field  intensity  5  as 
linear  combinations  of  harmonics  with  frequencies  equal  to 
mil,  where  m  is  an  integer  [13] 


a/2 


E=  Y 

(18) 

m=  — A/i 

A/j  +A/2 

S=  Y  Sme-'”'^^ 

(19) 

m=-(A/i-f  A/2) 

A/1+A/2 

N=Ne-i-  Y  ANme-'”'^* 

(20) 

m=~(A/i+A/2) 


where  Mi  and  M2  are  integers  chosen  according  to  the  de¬ 
sired  spectrum  to  be  modeled.  Fig.  1  illustrates  the  frequency 
relations  of  the  optical  field  of  the  pump,  probe,  and  multiple 
conjugate  waves  used  in  the  four  wave-mixing  theory.  The 
harmonics  of  the  field  intensity  5^  are  related  to  Em  as 

5„  =  53  (21) 

where  the  superscript  means  the  complex  conjugate.  Note 
that  both  S  and  AN  =  iV  ~  are  real  quantities.  The 
materia]  gain  g  including  the  nonlinear  gain  suppression  due 
to  the  Spectral-hole  burning  (SHB)  and  carrier-heating  can  be 
written  as  [14] 


9{N,  5)=po(A^e  +  AiV) 


1  -  ^  H^S^e-' 


—imQt 


m 

dgo 


—  imOt 


9o{K  +  AiV)  ^  (iV  -  iVe) 


=  9ot+ 9'  ^  AA^me"’ 


— imOt 


(22) 


(23) 


where  9o{N)  is  the  linear  part  of  the  material  gain,  is  the 
differential  gain,  Cah  and  gch  represent  the  relative  strengths 
of  the  spectral-hole  burning  and  carrier-heating  processes, 
respectively,  and  the  response  functions  of  the  two  processes 
are  [14],  [15] 


jysA  ^  * _ 

^  1  —  tmnT2 

^  (1  -  imfiri)  (1  -  imO,T2) 


(24) 

(25) 


where  T2  is  the  time  constant  of  the  carrier-carrier  scattering 
which  is  responsible  for  the  intraband  relaxation,  and  ti  is  the 
time  constant  of  the  carrier-LO  phonon  scattering  which  is 
responsible  for  the  heating  or  cooling  of  the  carriers.  Typical 
values  [16]  for  the  time  constants  tj  and  T2  are  650  fs  and  100 
fs,  respectively.  A  similar  expression  can  be  used  to  describe 
the  wave  number  of  the  cavity  k.  However,  for  simplicity,  we 
neglect  the  nonlinear  effects  in  k. 


where  is  the  linewidth  enhancement  factor.  Also  note  that 
Qe,  ^  5oc(l  ^sh  \E^e\  )  9ch  (27) 


By  equating  both  sides  of  the  rate  equations  with  the  harmonic 
expansions,  we  obtain  the  equations  for  the  normalized  carrier 
density  [ANm  =  {g'lgc)^Nm] 


(1  -  im^l)  AN 
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and  the  nonnalized  field  (Em  =  Em/Eg) 
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where 


p2  _  2  fiUL 

*  Con?  Vg9'T[, 


(30) 


I  a  =  {Dcon^VgE^  is  the  saturation  intensity,  t'  is  the 
dfferential  carrier  lifetime,  H  =  fir',  E{^  =  kcT'^Em/Eg, 

—  ^m/Egi  €gh  “  ^ahEgt  and  ”  (.9ch/9oe')Eg, 

Equations  (21),  (28),  and  (29)  form  a  system  of  nonlinear 
equations  which  can  be  solved  to  obtain  the  field  harmonics 
Em  inside  the  cavity  as  a  function  of  the  detuning  fi,  injection 
probe  beam  intensity  |J5in  p,  and  the  injection  current  density 
J  (i.e.,  |£e|^)-  The  general  model  is  then  further  simplified  to 
consider  two  of  our  experimental  conditions. 

1)  Multiple-Wave  Generation  From  Nearly  Degenerate 
FWM:  For  the  nearly  degenerate  FWM,  we  may  neglect  both 
the  SHB  and  CH  (i.e.,  €,h  =  0  and  pcfc  =  0).  The  material 
gain  is  then  a  linear  function  of  only  the  carrier  density.  We 
can  then  numerically  integrate  the  rate  equations  (14)  and 
(15)  with  the  initial  conditions 


E{Q)=Eg  (31) 

N{0)=Ng.  (32) 


k  =  k{Ne  +  AN)  After  some  initial  oscillations,  tiie  optical  field  E{t)  reaches 

g'  —  steady  state  and  becomes  a  periodic  function  of  time.  The 

«  Ale  -  Qe  —  Y  **”  ‘  (26)  haiTOonics  Em&  can  then  be  obtained  by  finding  the  coefficients 

of  the  Fourier  series  of  the  periodic  part  of  the  field. 
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2)  Weak  Probe  Beam  Injection:  Wth  a  weak  probe  beam 
injection,  i.e.,  [Einl  <  we  consider  only  three  waves: 
probe  pump  (£o)»  and  conjugate  (-E-i),  and  assume 
that  £o  =  Ee  and  AA^o  =  0. 

From  (21),  we  have  the  following  approximations 


5o«|£e|2  __  (33) 

(34) 

£_i  =£J  (35) 

5-2  =S5  «  0.  (36) 

Using  these  ^proximations  in  (28)  we  have 

A77i  =  Ai^lEi  +  Eelli )  (37) 

AN. I  =ANl  (38) 
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The  above  equations  for  ANi  and  AN^i  can  then  be  put 
into  (29)  for  m  =  ±1  to  form  a  system  of  linear  equations  of 
variables  Ei  and  Elj. 
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m.  Experiments  and  Comparison  with  Theory 
A.  Experimental  Setup 

For  our  experiments’,  we  use  an  index-coupled  DFB  laser 
with  a  lasing  wavelength  of  about  1.535  /xm  as  the  mixing 
device.  Fig.  2  shows  a  diagram  of  the  experimental  setup. 
Another  DFB  laser  similar  to  the  mixing  laser  is  used  to 
provide  the  probe  signal.  The  probe  signal  is  coupled  into 
a  fiber  and  passed  through  an  optical  isolator  to  eliminate 
injection  of  light  from  the  mixing  DFB  laser  into  the  probe 
DFB  laser.  The  probe  light  is  then  passed  through  a  20-dB 
fiber  coupler  to  monitor  the  power  of  the  probe  signal  injected 
into  the  mixing  DFB  laser.  Light  emitted  from  the  other  facet 


Detector 

20dB  Fiber  r-\ 

□0-®— — oao — 

optical  Isolator  DFB  Laser 

(Four-Wave  Mixing  Device) 

Fig.  2.  A  diagram  showing  the  experimental  setup  for  four-wave  mixing 
using  two  DFB  lasers.  The  probe  beam  is  coupled  from  the  first  laser  into  an 
optical  fiber  and  passed  through  an  optical  isolator  to  remove  unwanted  optical 
feedback.  The  probe  power  injected  into  the  second  DFB  laser  is  monitored 
using  a  20-dB  fiber  directional  coupler.  The  output  pump,  probe  and  conjugate 
waves  are  coupled  into  another  fiber  at  the  other  laser  facet  and  then  measured 
with  an  HP  (Optical  Spectrum  Analyzer. 

of  the  mixing  DFB  laser  is  coupled  into  another  fiber  and 
monitored  using  an  HP  optical  spectrum  analyzer  to  obtain  the 
output  powers  of  the  pump,  probe,  and  conjugate  signals.  The 
operating  temperatures  of  both  probe  and  mixing  DFB  lasers 
are  maintained  using  thermoelectric  controllers.  Tuning  of  the 
probe  signal  is  accomplished  by  changing  the  operating  bias 
current  and  temperature  of  the  probe  DFB  laser.  A  single  probe 
DFB  laser  has  a  typical  span  of  2.5-3.0  nm  with  a  temperature 
range  of  35  and  a  current  range  of  laser  threshold  to  60  mA. 
A  larger  tuning  range  is  achieved  by  using  several  different 
probe  DFB  lasers, 

B,  Characterization  of  Mixing  Laser 

In  order  to  obtain  the  gain,  refractive  index,  and  DFB  cavity 
parameters  in  our  models,  we  first  characterized  the  mixing 
DFB  laser  {Ith  =  8.4  mA)  using  the  amplified  spontaneous 
emission  (ASE)  spectrum  of  the  laser  just  below  threshold 
and  at  the  operating  bias  (40  mA,  20  °C)  used  in  the  four- 
wave  mixing  experiments.  The  details  of  this  procedure  are 
described  in  [17].  The  ASE  spectrum  shown  in  Fig.  3(a)  (40 
mA,  20  °C)  allows  us  to  obtain  the  gain  spectrum  from  the 
ratio  of  the  maximum  and  minimum  envelopes  based  on  the 
Hakki-Paoli  method  [18],  [19].  Fig,  3(b)  shows  the  extracted 
net  round-trip  gain  spectrum  at  both  8.2  mA  {Uh  =  8.4  mA) 
and  40  mA.  Far  above  the  bandedge,  these  two  spectra  are 
the  same  since  the  carrier  concentration  in  the  active  region 
is  more  or  less  pinned  above  threshold.  Near  the  bandedge, 
the  two  diverge,  however,  because  of  an  increase  in  the 
background  noise  which  leds  to  an  artificial  decrease  in  the 
extracted  Hakki-Paoli  gain.  To  obtain  the  true  intrinsic  loss 
level  and  the  quasi-Fermi  level  separation,  we  must  use  the 
plateau  of  the  net  gain  spectrum  at  8.2  mA.  Unfonunately,  near 
the  stopband  region,  the  strong  dependence  of  the  reflectivity 
of  the  DFB  grating  causes  the  extracted  data  to  deviate  sharply 
from  the  actual  gain  spectrum.  The  actual  gain  must  be 
interpolated  from  the  surrounding  data  [17].  Also,  an  effective 
refractive  index  spectrum  is  extracted  to  match  the  peaks  in 
the  ASE  spectrum  far  away  from  the  stopband  region.  The 
effective  refractive  index  in  the  vicinity  of  the  stopband  is 
interpolated  from  this  data.  The  cavity  length  of  the  laser  was 
measured  to  be  about  255  /xm.  Using  these  parameters,  the 
coupling  coefficient,  the  facet  reflectivities,  and  the  Bragg 
wavelength  can  be  found  by  performing  a  fit  of  the  ASE 
spectrum  as  done  in  reference  [17].  This  fit  is  not  shown 
because  the  same  parameters  are  used  later  in  this  section 
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Fig.  3.  (a)  The  measured  an^lified  spontaneous  emission  spectrum  of  the 
mixing  DFB  laser  well  above  threshold  (40  mA,  20  ®C).  Ttns  spectrum  is 
used  to  characterize  this  laser,  (b)  The  extracted  net  round-trip  gain  of  the 
mixing  DFB  laser  obtained  from  the  maxima  (/max)  and  the  minima  (/min) 
of  the  ASE  spectrum  such  as  given  in  part  (a)  Circles  show  die  net  round-trip 
gain  well  above  threshold  at  /  =  40  mA,  triangles  show  the  net  round-trip 
gain  just  below  the  laser  threshold  at  /  =  8.2  mA  {1th  =  8.4  mA). 


to  fit  the  probe  gain  spectrum  which  is  nearly  identical  to  the 
ASE  spectrum. 

C.  Power  Dependence  and  Multiple  Conjugate  Waves 

For  our  four-wave  mixing  measurements,  we  studied  both 
pump-probe  detunings  in  the  nearly  degenerate  and  the  non¬ 
degenerate  regimes.  For  small  detunings  we  observed  multiple 
conjugate  signals  for  a  single  injected  probe  beam.  Fig.  4 
shows  a  typical  spectrum  demonstrating  our  observations. 
Note  that  the  linewidths  of  the  spectra  shown  in  Fig.  4  are 
limited  by  the  minimum  resolution  of  the  optical  spectrum 
analyzer.  The  true  linewidths  are  much  narrower  than  those 
shown.  By  varying  the  coupling  of  the  emitted  light  from  the 
probe  DFB  laser  to  the  optical  fiber,  we  were  able  to  control 
the  power  of  the  probe  signal  injected  into  the  mixing  DFB 
laser.  This  enabled  us  to  study  the  dependence  of  higher  order 
conjugate  signals  on  the  power  of  the  probe.  Fig.  5(a)  shows 
the  measured  powers  of  four  of  the  conjugate  waves  as  a 
function  of  the  measured  probe  power.  The  slop>es  of  the  lines 
shown  on  the  log-log  plot  give  the  power  law  for  each  of 
these  conjugate  signals. 

The  dependence  of  the  conjugate  waves  on  the  pump 
power  is  much  more  difficult  to  measure  because  the  lasing 
wavelength  of  the  mixing  DFB  laser  will  change  as  the 
injection  current  is  changed.  So,  it  is  difficult  to  change  the 
pump  power  without  changing  the  pump-probe  detuning  at  a 


Fig.  4.  A  measured  optical  spectrum  showing  multiple  conjugate  signals  for 
a  pump  and  probe  with  a  small  detuning  for  nearly  degenerate  four-wave 
mixing. 

fixed  temperature.  In  our  measurements,  the  probe  signal  was 
tuned  to  a  fixed  detuning  from  the  pump  for  each  injection 
current  used  in  the  DFB  mixing  laser.  As  before,  the  power  of 
the  probe  signal  was  varied,  and  the  power  of  the  output  pump, 
probe,  and  conjugate  signals  were  measured.  The  dependence 
of  the  conjugate  powers  on  the  pump  power  was  obtained 
by  selecting  a  pump  power  common  to  the  range  measured 
for  each  bias  level  of  the  mixing  DFB  laser  and  interpolating 
the  value  of  the  conjugate  power  at  each  bias  level.  The 
interpolated  conjugate  powers  and  the  measured  pump  powers 
are  shown  in  Fig.  5(b). 

Our  model  for  the  nearly  degenerate  case  is  applied  to 
compare  with  our  experimental  data  shown  in  Rg.  5(a)  and 
(b).  We  change  th^normalized  injection  field  (£■„)  and  the 
injection  current  (Et)  and  calculate  the  norm^zed  power 
of  the  multiple  conjugate  waves.  The  results  are  shown  in 
Fig.  5(c)  and  (d).  A  linewidth  enhancement  factor  of  Oe  =  5 
is  used  in  the  calculation.  The  differential  carrier  lifetime  is 
treated  as  a  fitting  parameter  and  was  found  to  be  about  1 10  ps. 
The  saturation  intensity  Is  is  then  found  to  be  2.35  MW/cm^. 
With  an  emitting  area  of  1  pmx  0.02  /zm,  the  corresponding 
normalization  power  for  Fig.  5  (c)  and  (d)  is  0.47  mW.  The 
slopes  of  the  log-log  plot  of  the  theoretical  probe  and  pump 
power  laws  are  close  to  the  experimental  data.  The  relative 
magnitudes  of  the  conjugate  waves  show  some  disagreement 
between  experimental  and  theoretical  results.  A  better  estimate 
of  the  saturation  intensity  should  improve  the  comparisons. 

D.  Probe  Gain  and  Conjugate  Conversion  Efficiency 

Measurements  were  also  made  for  large  detunings  of  the 
pump  and  probe.  Fig.  6  shows  the  measured  output  spectrum 
for  several  detunings  of  the  pump  and  the  probe.  In  order 
to  maintain  a  constant  pump  power,  the  temperature  and 
operating  bias  current  of  the  mixing  laser  were  held  constant 
llie  probe  input  power  was  monitored  through  the  fiber 
coupler  and  was  typically  controlled  to  be  between  10  and 
100  times  smaller  than  the  pump  power  so  that  the  pump 
signal  of  the  mixing  DFB  laser  would  not  be  affected.  The 
probe  beam  was  tuned  from  about  550  GHz  below  to  about 
550  GHz  above  the  pump  frequency. 

Fig.  7  shows  the  ratio  of  the  measured  output  probe  power 
to  the  measured  input  probe  power.  We  have  not  included  the 
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(a) 


Wavelength  (nm) 


Fig.  6.  Measured  DFB  laser  spectra  with  several  probe  detunings  for  non- 
degenerate  four-wave  mixing. 


(c) 


Fig.  5.  Plots  of  the  dependence  of  the  measured  multiple  conjugate- wave 
ou5)ut  powers  on  (a)  the  probe  power  and  (b)  the  pump  power.  Tlieorctical 
results  for  the  normalized  conjugate-wave  intensity  as  a  function  of  (c)  the 
normalized  probe  intensity  and  (d)  the  normalized  pump  intensity  are  also 
shown.  All  of  the  powers  are  normalized  by  0.47  mW.  The  slopes  of  the 
log-log  plots  give  the  power  dependence  on  the  pump  or  probe  for  each 
conjugate  wave  shown. 

coupling  losses  from  the  fiber  to  the  mixing  DFB  laser  and 
from  the  laser  to  the  fiber,  but  this  factor  will  just  scale  our 


Fig.  7.  A  plot  of  the  gain  experienced  by  the  probe  through  the  DFB  cavity. 
Symbols  show  the  measured  data  and  the  solid  line  shows  the  best  fit  using 
parameters  extracted  from  the  laser  ASE  spectrum. 


Detuning  (GHz) 


Fig.  8.  Measured  four-wave  mixing  conversion  efficiency  (symbols)  as  a 
function  of  the  probe  detuning.  The  solid  curve  gives  a  fit  using  the  presented 
model  based  on  the  mean-field  theory. 

measured  ratio  by  a  constant.  To  match  this  spectra,  we  use  the 
same  parameters  as  before  from  our  fit  of  the  ASE  spectrum 
in  (11).  Fig.  7  shows  the  best  theoretical  fit  of  the  probe  gain 
yielding  =  rf  =  0.9,  =  320®,  i22  =  ^2  ~  ^2  = 

30®,  /c  =  66.7  cm“^  and  Ab  =  1.5334  /zm.  Fig.  8  shows  the 
conversion  efficiency  of  the  four-wave  mixing  defined  to  be 
the  ratio  of  the  measured  output  conjugate  power  to  the  input 
probe  power.  Again,  this  ratio  needs  to  be  modified  by  the 
same  coupling  loss  factor  as  in  the  probe  gain  in  order  to  give 
the  true  conversion  efficiency. 

The  weak-probe  wave  theory  is  used  to  study  the  conversion 
efficiency  of  FWM.  The  theoretical  conversion  efficiency  is 
calculated  as  the  product  of  a  coupling  constant  and  t^  ratio 
between  the  normalized  power  of  the  coi^gate  wave  (\E^i\^) 
and  that  of  the  injected  probe  wave  The  coupling 
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Fig.  9.  A  plot  showing  the  cavity  effects  of  the  conversion  efficiency.  The 
solid  line  corresponds  to  the  measured  ASE  power  spectrum  of  the  mixing 
DFB  laser  at  the  conjugate  signal  frequency. 

constant  is  used  to  account  for  the  coupling  of  the  probe  wave 
from  fiber  into  the  cavity  of  the  mixing  laser  and  the  coupling 
of  the  conjugate  wave  from  the  cavity  to  the  collecting  fiber. 
We  treated  the  coupling  constant  as  a  fitting  parameter  to 
fit  the  overall  shape  of  the  measured  conversion  efficiency. 
The  theoretical  results  of  the  conversion  efficiency  are  shown 
as  the  solid  line  in  Fig.  8.  The  coupling  constant  used  is 
0.01.  The  time  constants  used  are  650  fs  and  100  fs  for  ri 
and  r2,  respectively.  The  normalized  strength  of  the  CH  and 
SHE  nonlinearities  are  Cch  =  0.004  and  tsh  =  0.016.  The 
peaks  of  the  solid  line  in  Fig.  8  near  detunings  of  ±15  GHz 
are  due  to  the  resonance  of  the  carrier  density  modulation. 
These  peaks  are  the  same  as  the  resonance  peaks  which 
occur  in  direct  modulation  of  semiconductor  lasers.  When  the 
detuning  is  much  larger  then  the  CDM  resonance  frequency, 
the  contribution  from  CDM  to  the  FWM  process  decreases 
quickly  and  then  nonlinear  processes  become  dominant.  The 
CH  (ti)  is  the  dominant  nonlinear  mechanism  in  the  FWM 
process  for  the  detuning  range  of  our  experiment  data,  while 
the  contribution  from  the  SHE  (t2)  is  almost  negligible.  The 
conversion  efficiency  from  the  CH  contribution  has  a  much 
larger  bandwidth  then  the  CDM  process.  The  experimental 
data  show  some  resonance  peaks  which  are  not  predicted  by 
the  theory.  These  resonance  peaks  are  due  to  the  DFE  cavity 
resonance  effects  which  are  neglected  in  the  mean-field  theory 
and  are  explained  below. 

There  are  several  mechanisms  which  can  contribute  to  the 
overall  cavity  effects  in  the  four-wave  mixing  efficiency  spec¬ 
trum.  The  enhancement  of  the  probe  and  conjugate  beams  by 
the  cavity  is  the  most  obvious  effect.  In  Fig.  9,  the  conversion 
efficiency  is  superimposed  with  a  curve  showing  the  power  of 
the  ASE  spectrum  of  the  mixing  DFE  laser  at  40  mA  and  20 
®C  at  the  wavelength  of  the  conjugate  signal.  This  curve  gives 
a  very  good  indication  of  the  enhancement  that  the  cavity 
gives  to  the  conjugate  wave.  One  can  see  that  some  of  the 
general  trends  of  the  cavity  effects  in  the  conversion  efficiency 
correlate  very  well  to  this  curve.  It  is  difficult  to  correlate  a 
similar  curve  showing  the  enhancement  of  the  probe  wave  (the 
probe  gain  in  Fig.  6)  to  the  conversion  efficiency  because  large 
dips,  which  are  present  in  the  conversion  efficiency  spectrum, 
seem  to  occur  near  the  peaks  in  the  probe  gain  curve.  These 
dips  are  probably  caused  by  more  subtle  effects,  such  as  the 


phase  distribution  of  the  created  conjugate  waves  in  the  laser 
cavity.  These  effects  need  further  investigation. 

rv.  Conclusion 

In  this  paper,  we  have  presented  both  theoretical  and  ex¬ 
perimental  investigations  of  wavelength  conversion  in  a  DFE 
laser  using  four-wave  mixing.  For  detunings  in  the  nearly 
degenerate  regime,  we  observed  multiple  conjugate  signals. 
We  have  characterized  the  dependence  of  the  power  of  these 
conjugate  signals  on  both  the  probe  and  pump  powers.  Our 
measured  data  agree  veiy  well  with  the  theoretical  results 
based  on  the  mean-field  theory.  For  larger  detunings  in  the 
nondegenerate  regime,  we  present  a  detailed  measurement  of 
the  four-wave  mixing  conversion  efficiency.  We  point  out  the 
noticeable  DFB  cavity  effect  on  the  spectrum  and  correlate  the 
effects  with  the  cavity  enhancement  of  the  conjugate  wave. 
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The  piezoelectric  effects  on  the  optical  gain  of  wurtzite  GaN/AIGaN  QW  lasers  taking  into  account 
the  many-body  effects  are  presented.  The  self-consistent  model  with  piezoelectric  field  effect  shows 
that  band  structures  and  optical  gain  are  significantly  affected  by  the  piezoelectric  field  at  relatively 
low  carrier  densities.  The  peak  gain  is  redshifted  and  smaller  when  compared  to  the  flat-band  model 
without  piezoelectric  field  effect.  Only  gain  peaks  corresponding  to  Cl-HHl  and  Cl-LHl 
transitions  are  observed  in  the  investigated  range  and  transitions  for  C1-HH2  and  C1-LH2  are 
negligible  due  to  the  large  subband  energy  spacing  at  low  carrier  densities  and  small  matrix 
elements  at  high  carrier  densities.  At  high  carrier  densities,  the  self-consistent  model  shows  band 
structures  and  optical  properties  similar  to  the  flat-band  model  due  to  the  screening  effects.  ©  1998 
American  Institute  of  Physics.  [80003-6951(98)00924-3] 


The  wide-band-gap  semiconductors  including  GaN, 
AIN,  InN,  and  their  ternary  compounds  have  great  potential 
for  applications  in  high  power  and  optoelectronic  devices. 
Blue-green  light-emitting  diodes  based  on  wurtzite  GaN 
have  already  been  demonstrated,  and  room-temperature 
continuous-wave  operation  of  strained  InGaN/AlGaN  multi¬ 
quantum  well  (MQW)  laser  diodes  with  a  lifetime  of  many 
hours  has  recently  been  reported.*  With  the  current  progress 
in  these  blue-green  lasers,  considerable  efforts  have  been  de¬ 
voted  to  the  calculation  of  the  optical  properties  with  many- 
body  effects,  because  Coulomb  interaction  effects  among 
carriers  are  expected  to  be  important  due  to  the  large  exciton 
binding  energies  in  the  wide-band-gap  semiconductors.^"** 
Also,  recent  experimental  studies  show  that  III-V  nitride  het¬ 
erostructures  with  the  wurtzite  crystal  structure  pseudomor- 
phically  grown  in  the  (OOOl)-orientation  have  large  strain- 
induced  piezoelectric  fields.^  This  internal  field  can  greatly 
influence  the  optical  properties,^  and  therefore  the  piezoelec¬ 
tric  effect  in  addition  to  many-body  effects  should  be  taken 
into  account  for  the  design  of  high  performance  devices 
based  on  nitrides.  However,  most  calculations  of  the  elec¬ 
tronic  and  optical  properties  for  these  wurtzite  QW  structures 
have  been  calculated  mainly  under  the  flat-band  (FB)  model. 

In  this  letter,  we  report  new  results  on  the  piezoelectric 
effect  in  wurtzite  GaN/AIGaN  QW  lasers  taking  into  account 
many-body  effects.  Recently,  Chow  et  al^  treated  carrier- 
carrier  (C-C)  scattering  effects  on  the  optical  gain  at  the  level 
of  quantum  kinetic  theory  in  the  Markovian  limit  and 
showed  that  both  diagonal  and  nondiagonal  Coulomb  corre¬ 
lations  give  rise  to  important  modifications  of  the  optical 
gain.  Here,  the  line  shape  with  only  the  diagonal  C-C  scat¬ 
tering  contribution  given  by  Asada^  is  used  in  calculating  the 
optical  gain  as  an  approximation.  However,  the  line  shape 
function  differs  from  a  Lorentzian  in  that  the  broadening 
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factor  r  is  given  as  a  function  of  photon  energy  hw.  The 
band  structure  is  determined  based  on  a  self-consistent  (SC) 
solution. 

We  consider  a  GaN  QW  of  a  width  grown  along  the 
(0001)  (z-axis)  sandwiched  between  two  Alo.3Gao.7N  barri¬ 
ers.  The  piezoelectric  field  of  the  barrier  region  is  assumed  to 
be  zero.  The  valence  band  structure  is  calculated  by  using  the 
block-diagonalized  3X3  Hamiltonian®  based  on  the  k^p 
method.  The  material  parameters  for  GaN  and  AIN  are  taken 
from  Ref.  8,  except  for  those  listed  in  Table  I.  The  SC  band 
structures  and  wave  functions  are  obtained  by  iteratively 
solving  the  Schrodinger  equation  for  electrons,  the  block- 
diagonalized  3x3  Hamiltonian  for  holes,  and  Poisson’s 
equation.^  The  total  potential  profiles  for  the  electrons  and 
holes  are 

Vc(z)  =  Vc».(z)+k|Fz-|€|<^(z), 

V'„(z)  =  V„H,(z)  +  |elFz-lc|<^(2), 

where  VrM,(z)  and  V„h,(z)  are  the  square-well  potentials  for 
the  conduction  band  and  valence  band,  respectively,  F  is  the 
piezoelectric  held,  and  ^(z)  is  the  screening  potential  in- 


TABLE  I.  Parameters  used  for  the  calculations. 


Parameter 

GaN 

Defonnation  potentials  (eV) 

D, 

“1.70* 

D2 

6.30* 

D3 

s.oo* 

Da 

-4.00* 

Dielectric  constant(static)  e 

Elastic  stiffness  constants  (X  lO'*  dyn/cm^) 

10.0*’ 

C’li 

39.0* 

C12 

14.5' 

Piezoelectric  constant  (cm/V) 

-1.7X10"'“ 

•Reference  14, 
**Refercnce  15. 
^Reference  16. 
‘^Reference  5. 
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duced  by  the  charged  carriers,  which  satisfies  Poisson’s 
equation.  The  strain-induced  piezoelectric  field  is  given  by 


where  d-^  is  the  piezoelectric  constant,  c  is  the  static  dielec¬ 


tric  constant,  c,j  are  the  elastic  stiffness  constants,  and  is 
the  in-plane  component  of  the  strain  tensor.  The  wurtzite 
GaN  is  well  known  to  have  a  large  piezoelectric  field  along 
the  (OOOl)-direction  because  it  is  polar  along  the  c-axis. 

Using  a  spontaneous-emission  transformation  method 
and  a  Fade  approximation,^  the  optical  gain  g(oi),  taking 
into  account  the  many-body  effects,  is  given  by 


, _ rct,(<:||3w) _ 


(1-  Re  - fea»)/(r ,hi»)) 

(1-  Re  qj2+(Im  q* 


where  w  is  the  angular  frequency,  fiQ  is  the  vacuum  perme¬ 
ability,  e  is  the  dielectric  constant,  is  the  in-plane  wave 
vector,  Lj  is  the  well  thickness,  is  the  momentum 

matrix  element’^  in  the  strained  QW,/}^  and/^  are  the  Fermi 
function  for  the  conduction  band  states  and  the  valence  band 
states,  and  is  Planck’s  constant.  The  indices  /  and  m  de¬ 
note  the  electron  states  in  the  conduction  band  and  the  heavy 
hole  (light  hole)  subband  states  in  the  valence  band.  Also, 
£,„(it||)  =  £[(it||)-£:;(^l|)  +  £g+ A£sx+A£ch  «  the  renor¬ 
malized  transition  energy^  between  electrons  and  holes, 
where  Eg  is  the  band-gap  of  the  material,  and  A£sx  and 
AEch  are  the  screened  exchange  and  Coulomb-hole  contri¬ 
butions  to  the  band-gap  renormalization.  q{ki,h(u)  is  the 
excitonic  or  Coulomb  enhancement  of  the  interband  transi¬ 
tion  probability.^  In  the  above,  — 

-l-r„„(<:||,fiw),  where  T,,  and  are  broadening  factors 
related  to  the  decay  rate  of  the  electron  and  hole  wave  func¬ 
tions.  Here,  only  the  hole-hole  scattering  is  considered  be¬ 
cause  it  is  dominant  among  the  various  processes  such  as  the 
carrier-carrier  and  carrier-phonon  scattering.’ 

Figure  1  shows  the  potential  profiles  and  the  wave  func¬ 
tions  (Cl,  HHl,  and  HH2)  at  zone  center  for  the  self- 
consistent  model  with  the  piezoelectric  field  of  a  5  nm 
GaN/Alo.jGaoTN  QW  laser  at  (a)  N=2X  lO”  cm"^  and  (b) 
A^=10X10'’  cm"^,  respectively.  The  dotted  line  corre¬ 
sponds  to  potential  profiles  for  the  flat-band  model  without  a 
piezoelectric  field.  The  GaN  QW  is  under  0.724%  compres¬ 
sive  strain  and  the  strain-induced  piezoelectric  field  is  esti¬ 
mated  to  be  about  1.33  MV/cm.  This  large  field  induces  a 
spatial  separation  of  the  electron  and  hole  wave  functions 
and  creates  a  screening  electric  field  in  the  well  to  compen¬ 
sate  the  piezoelectric  field.  However,  the  SC  model  for  rela¬ 
tively  low  carrier  densities  shows  asymmetric  electron  and 
hole  wave  functions  relative  to  the  center  of  the  well  and  is 
spatially  separated  as  shown  in  Fig.  1  (a),  although  there  ex¬ 
ists  a  large  screening  electric  field  in  the  well.  Also,  it  is 
observed  that  the  ground  state  transition  energy  Cl -HHl  is 
redshifted  compared  to  the  FB  model  due  to  the  quantum- 
confined  Stark  effect.  For  example,  transition  energies  for 
the  FB  and  SC  models  are  3.55  and  3.50  eV  at 
N=2X  10'^  cm"^,  respectively.  On  the  other  hand,  as  the 


carrier  density  increases,  the  spatial  separation  of  the  elec¬ 
tron  and  hole  wave  functions  is  gradually  reduced.  The  SC 
model  at  N=  lOX  lO”  cm"^  shows  potential  profiles  and 
ground  state  transition  energy  Cl -HHl  similar  to  the  FB 
model.  Thus,  the  piezoelectric  effect  on  the  electronic  and 
optical  properties  is  dominant  in  a  range  of  relatively  low 
carrier  densities. 

Figure  2  shows  the  valence  band  structures  for  the  self- 
consistent  model  with  the  piezoelectric  field  of  a  5  nm 
GaN/Alo.3Gao.7N  QW  laser  at  (a)  2  X  and  (b) 

N=  10X10^’  cnl"^  respectively.  The  dotted  line  corre¬ 
sponds  to  the  band  structure  for  the  flat-band  model  without 
piezoelectric  field.  The  SC  model  at  A^=2xl0’^  cm'^ 
shows  a  larger  energy  spacing  between  the  first  two  sub¬ 
bands  (HHl  and  LHl)  and  higher  subbands  (HH2  and  LH2) 
than  those  of  the  FB  case.  The  increase  of  the  subband  en¬ 
ergy  spacing  reduces  the  carrier  population  in  the  higher  sub¬ 
bands.  However,  this  effect  is  compensated  by  the  decrease 
of  the  optical  matrix  element  due  to  the  spatial  separation  of 
the  wave  functions.  The  energy  spacing  decreases  with  in¬ 
creasing  carrier  density  and  the  SC  model  shows  band  struc¬ 
tures  similar  to  those  of  the  FB  model  at  a  carrier  density  of 


Length  (A)  Length  (A) 

FIG.  1.  The  potential  profiles  and  the  wave  functions  (Cl,  HHl,  and  HH2) 
at  zone  center  for  the  self-consistent  model  (solid  line)  with  the  piezoelectric 
field  of  a  5  nm  GaN/AlojGao 7N  QW  laser  at  (a)  N=2xl0‘’  cm"'’  and 
N=10X10’^  cm"^  respectively.  The  dotted  line  corresponds  to  potential 
profiles  for  the  flat-band  model  without  a  piezoelectric  field. 
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FIG.  2.  The  valence  band  structures  for  the  self-consistent  model  with  the 
piezoelectric  field  of  a  5  nm  GaN/Alo.3Gao.7N  QW  laser  at  (a) 
N  =  2X10'’  cm^^  and  (b)  N=\0'X\0^^  cm“^,  respectively.  The  dotted 
line  corresponds  to  the  band  structure  for  the  flat-band  model  without  piezo- 
elecuic  field. 

lOX  10*^  cm"*^.  This  is  due  to  the  increase  of  the  screening 
electric  field  to  compensate  the  piezoelectric  effect.  In  Fig.  3, 
we  plot  the  normalized  optical  matrix  element  at  zone  center 
using  the  self-consistent  model  as  a  function  of  carrier  den¬ 
sity  in  a  5  nm  GaN/Alo.3Gao7N  QW  laser.  For  comparison, 
using  the  flat-band  model,  we  find  that  the  normalized  opti¬ 
cal  matrix  elements  at  the  zone  center  are  0.249  and  0.248 
for  Cl-HHl  and  Cl-LHl,  respectively,  and  zero  for  Cl -HH2 
and  C1-LH2.  In  contrast,  the  SC  model  for  C1-HH2  and 
C1-LH2  has  matrix  elements  of  a  finite  magnitude  compared 
to  zero  for  the  FB  model.  Also,  the  matrix  elements  for  Cl - 
HH2  and  C1-LH2  are  larger  than  those  of  Cl-HHl  and  Cl- 
LHl  at  low  carrier  densities.  However,  transitions  for  Cl- 
HH2  and  C1-LH2  are  reduced  largely  due  to  the  increase  of 
the  subband  energy  spacing  as  shown  in  Fig.  2(a). 

Figures  4(a)  and  4(b)  show  the  optical  gain  spectra  and 
peak  gain  for  FB  and  SC  models  as  a  function  of  carrier 
density  in  a  5  nm  GaN/Alo.3Gao.7N  QW  laser.  The  optical 
gain  spectra  in  Fig.  4(a)  for  both  cases  have  only  one  peak 
which  corresponds  to  Cl-HHl  and  Cl-LHl  transitions.  The 
optical  gain  for  C1-HH2  and  C1-LH2  transitions  is  always 
smaller  than  that  of  the  Cl-HHl  and  Cl-LHl  transitions 
because  of  the  large  subband  energy  spacing  at  low  carrier 
densities  and  small  matrix  elements  at  high  carrier  densities. 


FIG.  3.  The  normalized  optical  matrix  clement  at  zone  center  using  the 
self-consistent  model  as  a  function  of  the  carrier  density  in  a  5  nm 
GaN/Alo aGao 7N  QW  laser.  Cl-HHl  and  Cl-LHl  are  on  top  of  each  other, 
so  are  C1-HH2  and  C1-LH2. 


10000 

"  :  S«M-consi8t«nt  (wtm  PZ 

6aN/AI„6a„N  (b) 

L,*  5  nm 

8000 

6000 

>7*^N(X10’*cm  *)  s . 

Rat-b^ 

/ 

4000 

f 

//\ 

2000 

0 

mmm 

/  X  Self-consistent 

p I  / . . i7.  i./r..  t  _ _ _ _ _ _ _ _ — . — ^ — 1 

0.34  0.35  0.36  0.37  0.38  0.39  0.40  0123456769  10 


Wavelength  {\irr\)  Carrier  density  (xlO^’cm*®) 


FIG.  4.  The  optical  gain  spectra  (a)  and  peak  gains  (b)  for  FB  and  SC 
models  as  a  function  of  carrier  density  in  a  5  nm  GaN/Alp jGap  7N  QW  laser. 


The  peak  gain  for  the  SC  model  at  low  carrier  densities  is 
redshifted  compared  with  the  FB  model. 

On  the  other  hand,  the  peak  wavelength  for  the  SC  cal¬ 
culation  approaches  that  of  the  FB  model  with  increasing 
carrier  density.  This  is  due  to  the  increase  of  the  screening 
electric  field  with  carrier  density.  As  the  carrier  density  in¬ 
creases,  the  peak  gain  in  Fig.  4(b)  for  the  SC  case  increases 
linearly  while  that  of  the  FB  case(dashed  line)  saturates 
gradually.  Thus,  it  is  expected  that  peak  gains  for  both  mod¬ 
els  coincide  with  each  other  at  a  high  carrier  density.  The  SC 
model  predicts  a  larger  threshold  carrier  density  than  that 
estimated  by  the  FB  model  in  the  investigated  range. 

In  summary,  the  piezoelectric  effect  on  the  optical  gain 
of  wurtzite  GaN/AlGaN  QW  lasers  taking  into  account  the 
many-body  effects  is  presented.  Our  calculations  show  that 
the  piezoelectric  effect  on  the  electronic  and  optical  proper¬ 
ties  plays  an  important  role  up  to  a  carrier  density  of  10 
XIO’^  cm‘^ 

This  work  was  supported  by  ONR  Grant  No.  N00014- 
96-1-0303.  One  of  the  authors  (S.H.P.)  was  supported  by 
Catholic  University  of  Taegu-Hyosung. 
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Abstract — A  theory  for  the  electronic  band  stmcture  and  the 
free-carrier  optical  gain  of  wnrtzite-strained  quantum-well  lasers 
is  presented.  We  take  into  account  the  strain-induced  band-edge 
shifts  and  the  realistic  band  structures  of  tbe  GaN  wurtzite  crys¬ 
tals.  The  effective-mass  Hamiltonian,  the  basis  functions,  tbe  va¬ 
lence  band  structures,  the  interband  momentum  matrix  elements, 
and  the  optical  gain  are  presented  with  analytical  expressions  and 
results  for  GaN-AlGaN  strained  quantum-weU  lasers. 
This  theoretical  model  provides  a  foundation  for  investigating  the 
electronic  and  optical  properties  of  wnrtzite-strained  quantum- 
well  lasers. 


I.  Introduction 

Recently,  wumite-strained  quantum  wells  have  been 
studied  intensively  in  an  effort  to  realize  GaN  blue- 
green  lasers.  With  the  current  progress  in  E-VI  and  III-V 
blue-green  lasers  [1],  [2],  theoretical  models  for  the  electronic 
and  optical  properties  of  these  quanmm-well  (QW)  structures 
become  very  important.  Although  optically  pumped  stimulated 
emission  from  GaN  was  first  observed  more  than  25  years  ago 
[3]  and  many  researchers  have  been  intensively  working  on 
the  electrical  injection  laser  diodes,  stimulated  emission  [4] 
and  laser  action  [2]  by  current  injection  have  only  recently 
been  achieved. 

There  have  been  more  reports  on  the  experimental  data  and 
theoretical  models  of  zinc-blende  compounds  than  those  for 
the  wurtzite  semiconductors.  Therefore,  wurtzite  semiconduc¬ 
tor  band  structures  and  optical  properties  are  less  understood, 
and  many  band-structure  parameters  are  not  available.  Fur¬ 
thermore,  a  GaN  layer  is  usually  grown  on  AlGaN  barrier  and 
cladding  layers,  which  are  grown  on  a  sapphire  substrate.  A 
significant  amount  of  lattice  mismatch  between  the  substrate 
and  the  AlGaN  layer  results  in  a  large  dislocation  or  defect 
density.  The  defect  density  may  broaden  the  gain  spectrum 
and  can  be  simply  accounted  for  in  the  linewidth  broadening 
parameter.  Strain  also  exists  in  the  GaN-AlGaN  layers.  As 
shown  in  Fig.  1  and  Table  I,  the  lattice  constant  of  GaN  is 
3.1892  A  and  that  of  AIN  is  3.112  A  [5],  [6].  The  corre¬ 
sponding  bandgap  energies  for  Al^Gai-xN  and  Ini_xGaxN 
ternary  compounds  with  bowing  effects  as  a  function  of  the 
lattice  constant  are  shown  as  solid  curves.  The  band  gap  of 
Al*Ini_xN  and  its  bowing  parameter  are  less  understood. 
From  the  only  available  data  [7],  there  seems  to  be  a  significant 
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Fig.  1.  Bandgap  energy  of  GaN,  AIN,  InN,  and  their  ternary  compounds 
or  possible  quaternary  compounds  as  a  function  of  the  lattice  constants.  The 
bowing  effects  for  the  AlGaN  and  InGaN  are  included  in  the  plot  of  tbe 
bandgap  energies  (solid  curves).  Because  not  enough  experimental  dau  for  the 
ternary  AUnN  conpounds  are  available,  a  dashed  line  is  shown  for  illustration. 

amount  of  bowing  effects.  We  use  a  dashed  line  to  show  the 
possible  range  of  band  gap  and  lattice  constant.  We  would 
expect  that  a  quaternary  compound  using  AlInGaN  will  have 
its  band  gap  and  lattice  constant  within  the  triangular  area 
in  Fig.  1.  This  means  that,  in  general,  various  amounts  of 
tensile  and  compressive  strains  are  possible  using  different 
combinations  of  ternary  and  quaternary  compounds  in  the 
wells  and  barriers.  However,  many  band-structure  parameters 
such  as  the  conduction-band  and  valence-band  deformation 
potentials  and  the  valence-band  effective-mass  parameters  are 
either  absent  or  cannot  be  measured  because  of  a  lack  of 
high-quality  pseudomotphically  strained  QW  samples. 

Previously,  the  effective-mass  Hamiltonian  for  wurtzite 
crystals  was  derived  using  an  invariant  method  [8]-[10]-  Th^ 
necessary  band-structure  parameters  to  describe  tiie  heavy- 
hole  (HH),  light-hole  (LH),  and  crystal-field  split-hole  (CH) 
bands  have  only  been  recently  obtained  by  fitting  the  band 
stmcture  of  the  effective-mass  Hamiltonian  to  that  of  a  first- 
principles  calculation  [11]  and  then  applied  to  the  calculations 
of  QW  band  stmcture  and  optical  gain  [12],  [13].  We  have 
recently  derived  the  effective-mass  Hamiltonian  based  on 
tiie  k  •  p  method  [14]  and  defined  important  band-stmcture 
parameters  in  a  similar  way  as  those  defined  using  the 
Luttinger-Kohn  method  [15]  for  the  zinc-blende  crystals.  Our 
model  is  attractive  since  the  Hamiltonian  of  a  wurtzite-strained 
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TABLE  I 

The  Material  Parameters 


Parameter 

GaN 

AIN 

Lattice  constants*^^  (A) 

a 

3.1892 

3.112 

c 

5.1850 

4.982 

Energy  parameters*"*^ 

E,{tV)  at  300K 

3.44 

6.28 

Ai  =  Acf  (meV) 

16^(22*) 

-58.5* 

(meV) 

12^  (11') 

A,  -  As  -  Ak/3  (meV) 

4 

6.80' 

ConductioD'band  effective  masses' 

mil  mo 

0.20 

0.33 

ml!  mo 

0.18 

0.25 

Valence-band  effective-mass  parameters' 

M 

-6.56 

-3.95 

-0.91 

-0.27 

5.65 

3.68 

Ai 

-2.83 

•1.84 

As 

-3.13 

-1.95 

Ao 

-4.86 

-2.91 

Deformation  potentials  (eV) 

0  (Interband) 

-8.16^ 

Ct  =  0.5o 

-4.08^ 

D, 

0.7^ 

D2 

2.1^ 

Da 

L4f 

D4 

-0.7^ 

Elastic  stiffness  constants*’*  (10”dyn/cm^) 

Ci3 

15.8 

12.0 

C33 

26.7 

39.5 

•The  linear  interpolation  formual  P(A]zGai-xN)  =  xP(AlN)  +{1— x)PGaN) 
is  used  for  most  of  the  parameters  for  AltGai_*N  barriers  except  for 
the  bandgap  energy,  which  has  a  bowing  parameter  [5]  of  0.98  eV. 
^^(AlxGai-zNXeV)  =  x(AlN)  +(1  ~  x)es(GaN)-0.98x{l  -  x), 

°  [51  *  [6j,  ^  [33],  ^  [27],  «  [11],  f  [14]. 

semiconductor  resembles  that  of  a  zinc-blende  semiconductor 

[16] -[18],  except  that  different  energy-splitting  parameters 
and  basis  functions  are  necessary.  The  six-by-six  Hamiltonian 
can  also  be  diagonalized  [19]  in  a  similar  fashion  as  that  for 
a  zinc-blende  structure. 

The  optical  gain  spectrum  for  a  strained  zinc-blende  QW 
laser  was  first  calculated  using  a  four-by-four  Hamiltonian 

[17] ,  [20]  ignoring  spin-orbit  coupling  and  then  using  a  six- 
by-six  Hamiltonian  including  the  spin-oibit  coupling  [18]. 
It  has  been  shown  [21]  recently  that  our  theoretical  results 
for  the  optical  gain,  the  refractive  index,  and  the  linewidth- 
cnhancement  factor  spectra  of  long-wavelength  (1.55  /xm) 
semiconductor  lasers  agree  very  well  with  our  experimental 
data. 

In  this  paper,  we  present  new  results  including  analytical 
expressions  for  the  interband  momentum  matrix  elements  and 


numerical  results  for  the  optical  gain  of  strained  wurtzite  QW 
lasers.  The  Hamiltonian  is  taken  from  a  block-diagonalized 
Hamiltonian  consisting  of  two  three-by-three  matrices,  which 
we  have  recently  derived  [19].  A  GaN-AlxGai-xN  strained 
QW  structure  is  considered  in  our  model  calculations.  We 
present  the  formulation  of  the  free-carrier  theory  for  the 
optical  gain  of  a  strained  QW  structure  in  Section  n.  The 
theory  for  electronic  properties,  such  as  the  conduction-band 
and  valence-band  structures,  the  strain-induced  band-edge 
shifts,  and  the  model  for  a  wurtzite  QW  heterostruture,  is 
then  presented  in  Section  HI.  We  then  show  in  Section  IV 
analytical  expressions  for  the  interband  momentum-matrix 
elements  and  the  final  expressions  for  the  optical  spontaneous 
emission  spectrum,  which  is  used  to  calculate  the  optical  gain 
coefficient.  Numerical  results  for  GaN-AlxGai_xN  strained 
QW  lasers  are  then  presented  in  Section  V.  The  conclusions 
are  given  in  Section  VI. 


n.  General  Formulation  of  the  Optical  Gain  Spectrum 

In  this  section,  we  present  the  general  formula  for  the  optical 
gain  of  a  strained  QW  laser  in  terms  of  the  electronic  properties 
such  as  the  conduction  band  and  the  valence  band  energies  and 
the  interband  transition  matrix  elements.  The  optical  material 
gain  coefficient  of  a  quantum  well  can  be  calculated  either 
from  Fermi’s  golden  rule  [17],  including  a  Lorentzian  line- 
shape  function  or  calculated  from  the  spontaneous  emission 
spectrum  [17],  [21],  [22]  using  a  fundamental  relation  between 
the  spontaneous  emission  and  gain  spectra.  Both  methods  are 
essentially  the  same.  However,  we  adopt  the  latter  because 
it  gives  more  physical  results  below  the  band  gap  than  the 
former  model  and  the  optical  gain  has  a  transparent  energy 
at  the  quasi-Fermi  level  separation.  The  spontaneous  emission 
rate  per  unit  volume  per  energy  interval  at  an  optical  energy 
hw  is  given  by  [21],  [22] 


'sp 


{M  = 


nrC€omlu)L^  ^  5,^ 

T*  u  0  *  17=1,1  <rs=C;, Ln,m 


(1) 


(2) 


where  q  is  the  magnitude  of  the  electron  charge,  mo  is  the 
electron  rest  mass  in  free  space,  c  and  cq  are  the  velocity 
of  light  and  permittivity  in  free  space,  respectively,  r\  is  the 
electron  spin,  which  can  be  either  T  or  i,  e  is  the  polarization 
vector  of  the  optical  electric  field,  rir  and  Lz  are  the  refractive 
index  and  well  width  of  the  quantum  well,  and  ^7  is  the  half 
linewidth  of  the  Lorentzian  function  (7  =  (0.1  ps)'^  is  used 
in  this  p^r  for  sin?>licity).  MX^(fct)  = 
is  the  momentum  matrix  element  for  transitions  between  the 
conduction-band  state  ’*’^4,(2)  and  the  valence-band  state 
p  is  the  momentum  operator.  For  convenience, 
gtp(hu)  is  defined  to  have  the  same  dimension  (1/cm)  as  the 
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material  gain  coefficient.  The  expression 


=  Kikt)  -  (3) 


is  the  interband  energy  of  the  conduction  band  and  the  valence 
band.  The  Fertni-Dirac  distributions  for  the  electrons  in  the 
conduction  and  valence  bands  are 


mt)  = 


f:M  = 


l+exp(£i^) 
l  +  exp(^^^) 


(4) 

(5) 


where  Fc  and  are  the  quasi-Fermi  levels  for  electrons  and 
holes. 

The  material  gain  can  be  determined  from  the  spontaneous 
emission  rate  [21],  [22]  using 


g{hu})  =  9lj,{hu}) 


exp 


( 


huj^AF\ 
ksT  ) 


(6) 


for  e  =  TE  {—x  or  y,  which  is  perpendicular  to  the  c  axis)  or 
€  =  TM  (=2,  which  is  parallel  to  the  c  axis)  polarization.  The 
modal  gain  =  Tg  is  obtained  by  multiplying  the  material 
gain  by  the  optical  confinement  factor  F.  The  separation  of 
the  quasi-Fermi  levels,  AF,  depends  on  the  injected  carrier 
densities 


AF  =  Fc-F^..  (7) 

The  quasi-Fermi  levels  will  be  discussed  in  Section  FV. 

ni.  The  Electronic  Band 
Structures  and  Wave  Functions 

Consider  a  strained  GaN-iayer  wurtzite  crystal  pseudomor- 
phically  grown  along  the  c  axis  (z  axis)  on  a  thick  AlxGai_xN 
layer.  TTie  strain  tensor  f  in  the  well  region  has  the  following 
elements: 

do  ^  a 

€xx  ~  ^yy  ~ 

2Ci3 

€z2  —  ~  ^  Cxz 

^xy  ^  ^yz  “  ^zx  “  0  (8) 

where  ao  and  a  are  the  lattice  constants  of  the  AlxGai-^N 
barrier  and  the  GaN  well  layers,  respectively.  C13  and  C33 
are  the  stiffness  constants  of  the  GaN  well  layer.  The  general 
strain-stress  relation  for  the  hexagonal  crystal  can  be  found 
in  [23]. 

The  conduction  bands  can  be  characterized  by  a  parabolic- 
band  model,  and  the  effective-mass  Hamiltonian  can  be  written 
as 

where  the  wave  vector  kt  =  =  -idjdz,  and  m* 

and  m*  are  the  electron  elfective  masses  perpendicular  (t) 
and  parallel  (z)  to  the  growth  direction,  respectively.  The 
hydrostatic  energy  shift  in  the  conduction  band,  F C£(^)»  is  zero 


(a)  Unstrained  (b)  Comprcssively  strained 

wurtzite  GaN  wurtzite  GaN 


Fig.  2.  Schematic  diagram  showing  the  bulk  band  structures  of  (a)  an 
unstrained  GaN  wurtzite  semiconductor  and  (b)  a  comprcssively  strained  GaN 
layer. 


Alj^Gaj  GaN 


T 

AE,. 


^\z) 


-  -E^’^GaN) 


Eg+Ai+A2 


AEy 


-  E5’\GaN)=0 


•^\z) 


Fig.  3.  The  conduction-  and  valence-band  energy  profiles  of  a 
GaN-Alo.3Gao.7N 'quantum  well  with  a  well  width  Lu-  before  adding 
the  strain-induced  band-edge  shifts.  The  reference  energy  is  chosen  to  be 
Ej  of  an  unstrained  GaN  material. 


in  the  barrier  regions  and  is  due  to  a  compressive  strain  in  the 
well  region 


Pc€  *—  ^cz^zz  0.ct{^xx  “b  ^yy)  (10) 

where  dez  and  det  are  the  conduction-band  deformation  po¬ 
tentials  along  the  c  axis  and  pcipendicular  to  the  c  axis, 
respectively.  We  take  =  flcf(=  ^c)  for  simplicity  and 
assume  a  value  of  ^4.08  eV  from  a  recent  fit  [14]  to  the 
available  experimental  data.  Since  the  interband  deformation 
potential  is  negative,  the  band  gap  increases  with  an  increasing 
compressive  strain,  as  shown  in  Fig.  2.  The  conduction-band 
edge  is  shifted  upward  for  compressive  strain,  as  in  a  thin  GaN 
layer  sandwiched  between  two  thick  AlGaN  layers.  E^iz)  in 
(9)  is  the  QW  potential  profile  of  the  unstrained  conduction- 
band  edge  and  is  given  by  E^{z)  =  ^^(GaN)  for  z  in  the 
well,  and  =  F°(AlxGai-*N)  for  z  in  the  barrier,  as  shown  in 
Fig.  3.  We  choose  the  reference  energy  at  E^{G^N)  =  0, 
which  is  the  band  edge  of  the  CH  band  of  an  unstrained 
wurtzite  GaN  before  including  the  spin-orbit  interaction  [14], 
E^  =<  z\Ho\z  >.  Note  that  the  heavy-hole,  light-hole,  and 
conduction-band  edges  are  referenced  to  the  same  energy 
value,  F2(GaN),  following  the  convention  for  the  wurtzite 
structure.  In  this  convention,  ^^(GaN)  —  FS(GaN)  =  Fp  -f 
Ai  +  A2* 
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The  wave  function  of  the  nth  conduction  subband  can  be 
written  as 

=  (11) 

where  kt  is  the  transverse  wave  vector,  rt  is  the  two- 
dimensional  position  vector  and  A  is  the  area  of  the  quantum 
well  in  the  x-y  plane,  (^>n(z)  is  the  envelope  function  for  the 
^-dependent  part  of  the  nth  conduction  subband,  and  7?  is  the 
electron  spin  (t  or  [). 

The  conduction-band-edge  energy  =  E^{kt  =  0) 
is  determined  by  solving  the  effective-mass  equation  for 
electrons 

H^{kt  =  0,  =  -id/dz)4>n{z)  =  E^(pn{z)  (12) 

with  the  envelope  function  normalized  as 

J  =  1.  (13) 

We  use 

+  (14) 

to  obtain  the  conduction-band  structure,  where  ^  is  the 
electron  effective  mass  in  the  well  region  parallel  to  the 
quantum  well. 

The  bulk  valence-band  structure  can  be  determined  by 
finding  the  eigenvalues  of 

det[Hij{k) {k)]^0  (15) 

where  k  is  a  real  wave  vector  and  the  envelope  functions  are 
plane  waves  with  the  same  wave  vector  k.  The  Hamiltonian 
for  the  valence-band  structure  has  been  derived  recently  by 
the  k  •  p  method  [14],  and  the  matrix  elements  are  defined 
consistently  using  our  chosen  basis  functions.  We  found  that 
the  six-by-six  Hamiltonian  can  be  block-diagonalized  using  a 
unitary  transformation  [19] 

liv  /i.\  _  /i^\ 

He.eW-[  0 

where  ^zxz  three-by-three  matrices  defined  as 

'F  Kt  -iHt  " 

=  Kt  G  A-iHt 
iHt  A^iHt  A 
f  F  Kt  iHt  * 

Kt  G  A  +  iHt  .  (17) 

-iHt  A -iHt  A 

Note  that  the  relations  hold:  H^  =  {H^Y  =  {H^Y  and 
the  reference  energy,  £'5,  which  should  appear  in  all  the 
diagonal  terms  of  the  matrices,  has  been  set  to  zero  for  a 
bulk  semiconductor.  The  matrix  elements  contain  the  general 
expressions  for  a  strained  wurtzite  semiconductor 


A,  =  Di«  zz  •f  £2(6  XX  ^yy) 

6  ^  6k  Of 

+  Aik^) 

6f  =  Dzt 

ZZ  +  Dj^i^xx  +  €yy) 

Kt  =  ^Azkt 

ZTTIq 

±.2 

Hi  =  ^Asktk, 

2mo 

A  =  V2A3  (18) 

where  kt  =  +  k^  is  the  magnitude  of  the  wave  vector 

in  the  kx-ky  plane.  Note  that  the  bases  for  the  Hamiltonian 
in  (16)  are  defined  as 

|1)  =  0*1^11  T)  +  a|yi_i  1) 

m  =  m-i  T)+/?*|yu  i) 
i3)  =  riyio  T)+mo  1) 

|4)  =  a*|yn  T>-aiyi-i  1) 

\5)  =  m-i  T)-rini  1) 

|6)  =  -r|yio  ])  +  0\Yio  1) 

Yn  =  ^\X  +  iY) 

Yio  =  \Z) 

yi_i  =  ;^|j-iy).  (19) 


1  + 


and  4>  =  tan“^(A:a:/fcy)  is  the  azimuthal  angle  in  the  kx- 
ky  plane  and  covers  a  range  between  0  and  27r.  The  wave 
functions  Yimil  -  l,m  =  -1,0,1)  in  (19)  are  the  spherical 
harmonics  of  the  p  state  of  hydrogen  atom  [17]  in  cartesian 
coordinates  following  the  phase  convention  of  [24] .  The 
i4iS  are  the  band-structure  parameters,  which  play  the  same 
role  as  the  Luttinger  parameters,  71,72^  and  73,  of  zinc- 
blende  crystals.  Ai  is  the  crystal-field  split  energy  and  is 
caused  by  the  anisotropy  of  the  wurtzite  symmetry.  It  is 
the  difference  between  the  eigenvalues  of  the  unperturbed 
Hamiltonian,  Hq,  of  the  AT-like  (or  y-like)  states  and  the  Z- 
like  states  before  the  spin-orbit  coupling  effect  is  included, 
{X\Ho\X)  =  {Y\Ho\Y)  =  £^+Ai,  and  {Z\Ho\Z)  =  E^.  £. 
is  the  reference  energy.  A2  and  A3  account  for  spin-orbit 
interactions.  It  should  be  pointed  out  that  under  the  cubic 
approximation  [8],  [11],  the  following  relations  hold  among 
the  parameters  >ltS  and  AiS: 


Ai  —  A2  =  —.A3  =  2^4,  .A3  -h  4A5  =  y/2A^ 

A  -  £2  =  -Dz  =  2jD4,  A2  =  A3.  (21) 


E  =  Ai  +  A2  -f*  A  -|-  5 
G  “  Ai  —  A2  -l-  A  4-  ^ 
A  =  Ajt  4“  Ae 


Therefore,  five  band-structure  parameters  such  as 
A.i,A2,>l5,  Ai,  and  A2,  plus  two  deformation  potentials  Di 
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and  D2y  are  necessary  for  the  calculation  of  the  valence- 
band  structures.  Sometimes,  the  spin-orbit  split  energy 
Aso  =  3A2  =  3A3  is  used. 

For  a  QW  structure,  the  corresponding  material  parameters 
in  the  barrier  and  the  well  regions  have  to  be  used  in 
the  Hamiltonian.  The  band-edge  discontinuities  have  to  be 
added  properly  in  the  effective-mass  equation.  The  hole  wave 
function  in  a  quantum  well  can  be  written  as 


+  flW(z;fc*)|3))  (22) 

^{gW{z-,kt)\4)+gl^Hz;kt)\5) 

+  £Hz;kt)\6))  (23) 


where  gm-i  =  1.2.  -.e  are  the  envelope  functions  of  the 

mth  valence  subbands,  which  follow  the  normalization  rules 


^  J  dz\gli,Hz:kt)f  =  h  J  dz\g!^Hz:kt)f‘  =  1. 

(24) 

The  valence  subband  structure  for  the  upper  Hamiltonian 
E^{kt)  can  be  determined  by  solving 


for  i  =  1,2,3.  Similarly,  the  band  structure  for  the 
lower  Hamiltonian,  E^{kt),  can  be  obtained  using  Hlj  + 
6ijE^{z),  E^Xz)  is  the  QW  potential  profile  of  the  unstrained 
v^ence-band  reference  energy,  and  is  given  by  E^{z)  = 
£2(GaN)  for  z  in  the  well,  and  =  £2(^xGai-xN)  for  z  in 
the  barrier.  We  choose  the  reference  energy  at  ^^(GaN)  =  0. 
The  band-edge  discontinuity  for  the  valence  band  is  AjE^  == 
QvJ^Eg  where  AEg  is  the  bandgap  difference  bewteen 
the  unstrained  Al^Gai-xN  and  GaN.  The  partition  ratio  is 
assumed  to  be  Qv  =  0-33  for  the  valence  band,  which  follows 
the  experimental  results  for  the  GaN-AlN  heterojunction  [25]. 
Thus,  we  have  the  QW  potential  E^{z)  =  0  in  the  well,  and 
£^2(2:)  =  -QvAEg  in  the  barriers.  The  bandgjq)  energy  of 
AlxGai^xN  has  a  bowing  parameter  b  =0.98  eV  and  is  given 
in  Table  I  together  with  odier  material  parameters,  which  are 
linearly  interpolated  from  the  corresponding  GaN  and  AIN 
binary  wurtzite  semiconductors.  Note  that  the  strain-induced 
band  edge  shifts  are  taken  into  account  in  the  diagonal  tenxis 
of  the  Hamiltonian  matrices. 

Since  and  differ  only  by  a  change  of  signs  in  front 
of  all  the  terms  with  a  single  derivative  with  respect  to  z 
(i.e.,  djdz  in  Ht\  it  is  easy  to  see  that  for  a  quantum  well 
with  a  reflection  symmetry  E^{z)  =  E^{-z),  the  upper  and 
lower  Hamiltonians  have  the  same  valence-band  structures, 
and  the  wave  functions  are  related  by 
gm\z\kt)  =  9m\-z;kt),9m\z\kt)  =  and 

9m\z\kt)  =  9m\-z;kt).  At  the  zone  center  {kt  =  0),  the 
HH  subbands  are  decoupled  from  the  LH  and  CH  subbands, 


while  there  is  always  a  coupling  between  the  LH  and  CH 
bands. 

rv.  iNTERBAND  MOMENTUM-MaTRIX  ELEMENTS 
FOR  Gain  and  Spontaneous  Emission  Spectra 

Using  the  expressions  for  the  basis  |uf )  in  (19),  we  obtain 

+2cos(2(^){())„|p^^)((?!'nl5^^)}  (26) 

for  7]  =  both  T  and  J.,  where  we  have  used  the  following 
momentum  matrix  elements  from  the  basis  functions: 

TE(e  =  x)  polarization 

(STb.|i)  =  -^(5|Pxi:^> 

{S  T  bx|2)  =  -^{S\p.\X) 

(5T1p.|3>  =  0 

(ST|Px|4)  =  -^(5|Px|X> 

(5  T  1Px|5>  =  -^iS\P.\X) 

{S  T  bx|6)  =  0  (27) 


TM(e  =  z)  polarization 

(5Tbx|l)  =  0 

(S  T  bx|2>  =  0 
(5Tbx|3>  =  r(5bx|2) 

(5Tbx|4)  =  0 
{S  T  bx|5)  =  0 

(5Tbxl6)  =  -;3-(51Pxl2>  (28) 

for  p  =  t-  Simlar  expressions  can  be  derived  for  the  electron 
spin  p  =  I  in  the  conduction  band. 

Since  the  only  ^dependent  term  in  the  integrand  of  (2) 
appears  in  the  matrix  elements,  the  integration  of  the  cos(2i^) 
term  over  ^  vanishes.  We  then  obtain  two  scalar  polarization- 
dependent  matrix  elements,  which  result  from  the  <j>  integration 
of  the  momentum  matrix  elements: 

TE-polarization  (e  =  i  or  p  ±  c  axis): 

+ <^„  Is®)’} 

for  a  =  (7 

= Is®)’ + <^,iss>)’} 

for  (7  =  1.  (29) 


for  £7  =  t/ 

for  (7  =  L. 

(30) 
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We  define  two  equivalent  energy  parameters,  Ep.  and  Epxt  in 
terms  of  the  band-edge  momentum  matrix  elements 

^^h^  =  \{S\p,\X)f  (31) 


for  the  interband  transition  in  a  similar  way  to  that  of  the  zinc- 
blende  crystals  [26],  for  which  only  one  parameter  is  required. 
These  energy  parameters  can  be  expressed  in  terms  of  the 
band-structure  parameters  by  generalizing  Kane’s  model  [17], 
[26]  to  those  for  the  wurtzite  crystals,  taking  into  account  the 
hexagonal  symmetry  [14].  We  obtain 

(£g  +  Ai+A2)(JE:a  +  2A2)-2Ai 

(£g  +  2A2) 

£g  [(£g  +  Ai  +  A2)(£s  +  2A2)  -  2Al] 

(£s  +  Ai  +  A2)(£s  +  A2)-A|  • 

(32) 

Therefore,  the  numerical  values  for  the  interband  matrix  el¬ 
ements  can  be  calculated  using  the  parameters  listed  in  Ta¬ 
ble  I.  It  should  be  pointed  out  that  no  experimental  data 
for  the  magnitudes  of  the  interband  optical  matrix  elements 
are  available  due  to  the  lack  of  a  good  calibration  scheme 
and  probably  high-quality  GaN  samples;  the  above  analytical 
model  does  provide  a  resonable  estimation  for  the  interband 
matrix  elements.  For  GaAs  materials,  a  similar  formula  [17] 
for  the  energy  parameter  Epz  =  Epx  of  the  zinc-blende 
crystals  following  Kane’s  moidel  gives  a  value  of  about  21.4 
eV,  which  is  about  17%  smaller  than  the  reported  experimental 
value  of  25.7  eV.  We  would  expect  a  similar  percentage 
of  errors  in  the  numerical  values  when  using  the  above 
expressions  for  wurtzite  crystals. 

The  final  expression  for  gspi^u)  in  the  spontaneous  emis¬ 
sion  rate  is 


^  —  T1,TT1 


(33) 


The  total  spontaneous  emission  rate  per  energy  interval  per 
unit  volume,  rsp(ftw),  is  given  by 


Up{hu)  = 


(2gfp+gs%) 

TT^hc^  3 


(34) 


which  uses  the  fact  that  the  momentum-matrix  element  of  the 
total  spontaneous  emission  is  the  angular  average  of  two  TE- 
polarization  components  and  one  TM-polarization  component 

|Msp|2  =  i(2|M,|2  +  lM.p).  (35) 


For  a  given  injection  level,  the  electron  and  hole  concentra¬ 
tions  are  related  by  the  charge  neutrality  condition 


n-¥  Na-P’^  Nr>  (36) 


where  Na  and  Nd  are  the  ionized  acceptor  and  donor  concen¬ 
trations,  respectively,  in  the  active  region  and  are  assumed  to 


be  zero  for  an  undoped  active  layer.  The  electron  concentration 
n  and  the  hole  concentration  p  are  related  to  the  quasi-Fermi 
levels,  Fc  and  Ft,,  in  (4)  and  (5),  by 


ksTml 

wh^Lw 


E>“[i+ 


,\F,-EimykBT-\ 


c=U,L  m  '' 


(37) 

(38) 


V.  Numerical  Results  and  Discussions 

We  consider  first  a  layer  of  GaN  under  either  a  compressive 
or  tensile  strain,  depending  on  the  barrier  materials,  with 
the  convention  that  the  in-plane  strain  Cxx  is  negative  for 
compression  and  positive  for  tension.  The  conduction  band- 
edge  energy  in  the  GaN  bulk  layer  is 

£c  =  £°  +  Ai+A2  +  £g  +  Pc,-  (39) 


The  valence  band-edge  energies  can  be  obtained  by  finding  the 
eigenvalues  of  the  three-by-three  Hamiltonians,  (17),  evaluated 
at  the  zone  center  (k  =  0) 


Fhh  =  -f  Ai  +  A2  +  +  Ac, 


+ 


+  24S 


~  A2  +  0, 


')  +2AI 


3- 


(40) 


The  band  edge  energies  of  the  conduction  band,  Ec,  the  HH, 
the  LH,  and  the  CH  bands  are  shown  in  Fig.  4(a)  and  (b).  With 
increasing  compressive  strain,  the  bandgap  energy  Ec  -  Fhh 
increases,  where 


Fc  -  Fhh  =  Fp  -f  Fc<  ~  (Ac  -f  6^) 

=  Fp  -f  (acz  -  Di-  Ds)€zz 
+  (Uct  —  F2  -  D4){€xx  +  6yy).  (41) 


It  should  be  pointed  out  that  most  measurements  of  the 
bandgap  energy  as  a  function  of  an  externally  applied  hydro¬ 
static  pressure  can  only  provide  the  total  inteiband  deformation 
potential,  (acz  -Di-  Dz)  or  (uct  -  F2  -  F4),  and  it  is  not 
possible  to  dintiguish  these  two  numbers,  or  individual  param¬ 
eters.  The  hydrostatic  pressure  also  leads  to  the  same  sign  for 
^nd  tzzy  whilc  a  biaxial  strain  in  the  GaN-AlGaN 
QW  structure  has  a  sign  difference  between  the  in-plane  strain 
(txx  =  Cyy)  and  the  strain  along  the  growth  direction  (c^^).  At 
zero  strain,  Fhh  =  Ai  -f  A2  =  0.020  eV  using  the  parameters 
listed  in  Table  I.  These  parameters  were  obtained  by  fitting  to 
the  experimental  data  collected  in  [27].  With  increasing  tensile 
strain,  which  would  require  an  AUnN  or  AUnGaN  barrier 
material  (see  Fig.  1),  the  band  gap  of  the  tensile-strained  GaN 
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Fig.  4.  The  band  edges  of  (a)  the  conduction  band  and  (b)  the  HH,  the  LH, 
and  the  CH  bands  as  a  function  of  the  compressive  and  tensile  strains  for  a 
GaN  crystal.  The  reference  energy  is  set  at  of  an  unstrained  GaN  material. 


will  have  a  smaller  band  gap  than  the  unstrained  GaN.  The 
strong  coupling  of  the  LH  and  CH  bands  also  splits  the  energy 
at  a  value  of  strain  near  0.004  (or  0.4%).  Since  the  band- 
structure  parameters  such  as  the  deformation  potentials  and 
effective-mass  parameters  for  the  valence  band  are  available 
only  for  AlxGai-xN  compounds,  we  present  numerical  results 
for  the  GaN-AlGaN  QW  structures  in  the  following  numerical 
examples.  This  system  has  only  compressive  strain  instead  of 
tensile  strain  in  the  GaN  well  layers. 

The  valence  subband  dispersion  curves  are  shown  in 
Fig.  5(a)  and  (b)  for  two  QW  structures  with  well  widths 
iu,  =  26  A  and  50  A,  respectively.  We  choose  62  A-26  A-^2 
A  and  50  A-50  A-50  A  as  the  barricr-weU-bairier  widths  and 
use  an  interval  of  2  A  in  a  finite-difference  method  [28]  to 
find  the  electron  and  hole  wave  functions  when  we  solve  the 
coupled  differential  equations.  The  label  of  the  HH,  the  LH, 
and  the  CH  subbands  is  determined  by  the  wave  functions  and 
their  components  at  the  band  edges,  kt  =  0.  For  HH  subbands, 
the  wave  function  component,  or  is  dominant.  For 
LH  subbands,  or  and  for  CH  subbands,  g^^^  or  g^^\ 
are  dominant.  Fig.  5(a)  and  (b)  also  shows  that  a  thick  well  has 
closer  spacings  between  the  subband  energies  than  those  of  the 
thin  well.  For  the  quasi-two-dimensional  nature  of  the  density 
of  states  to  play  a  significant  role  on  the  optical  gain  of  a  QW 
laser,  a  thin  well  is  desirable,  especially  because  the  electron 
and  hole  effective  masses  of  GaN  are  very  heavy.  However,  a 


k  (l/A) 


(a) 


k,(l/A) 

(b) 


Fig.  5.  The  valence  subband  structures  of  a  strained  GaN-Alo.3Gao.7N 
quantum  well  with  well  widths:  (a)  Lu-  =  26  A  and  (b)  Lyj  =  50  A.  The 
reference  energy  is  set  at  £2  of  an  unstrained  GaN  material. 

thin  well  reduces  the  modal  gain;  therefore,  multiple  quanmm 
wells  are  required. 

We  show  in  Fig.  6(a)  and  (b)  the  dispersion  of  the 
momentum-matrix  elements  as  a  function  of  the  in-plane 
wavenumber  kt  on  the  plane  for  the  first  conduction 
subband  (Cl)  to  the  first  heavy-hole  (HHl),  the  light-hole 
(LHl),  and  the  crystal-field  split-hole  (CHI)  transitions  and 
their  polarization  dependences  for  Lw  =  26  A.  Both  the  Cl- 
HHl  and  Cl-LHl  transitions  have  strong  TE  components  near 
the  zone  center,  while  the  Cl-CHl  transition  has  a  dominant 
TM  component,  which  can  be  larger  than  the  other  transitions 
for  the  TE  polarizations.  However,  since  the  compressive 
strain  in  the  quantum  well  makes  the  top  few  valence  subbands 
HH-  or  LH-like,  the  dominant  transitions  for  the  compressive 
GaN-AlGaN  transitions  will  be  TE  polarization.  This  result 
is  different  from  that  of  a  zinc-blende  structure,  for  which 
the  conduction  to  LH  subband  transition  has  a  dominant  TM 
component,  and  the  LH  subband  can  be  shifted  downward 
away  from  the  top  HH  subband  with  a  compressive  strain. 
For  a  wurtzite  structure,  the  HH  and  LH  subbands  tend  to 
shift  in  the  same  direction  by  almost  the  same  amount  of 
energy,  resulting  in  competition  for  the  hole  population  in 
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(a) 


(b) 

Fig.  6.  The  normalized  momentum-matrix  elements  for  the  first  conduction 
(Cl)  subband  to  the  HHl,  the  LHl,  and  the  CHI  subband  transitions  for  (a) 
TE  polarization  [on  the  x-y  plane,  which  is  perpendicular  to  the  c  axis  (=? 
axis)]  and  (b)  TM  polarization  (along  the  c  axis)  for  a  well  width  Lu-  =  26  A. 

each  subband.  Both  the  HH  and  LH  bands  contribute  to  the 
same  (TE)  polarization,  yet  at  different  energies  result  in  a 
broader  gain  spectrum  than  those  of  the  zinc-blende  structures. 
The  momentum-matrix  elements  in  Fig.  6(a)  and  (b)  also  show 
that  the  coupling  between  the  LHl  and  CHI  subbands  is  weak, 
while  the  coupling  between  the  LHl  and  HH2  becomes  strong 
at  kt  near  0.07  as  can  be  checked  from  the  valence  band 
structure  in  Fig.  5(a)  and  (27>-(28). 

Fig.  7(a)  and  (b)  show  the  material  gain  spectrum,  g{fuJ), 
and  the  spontaneous  emission  spectrum,  which  is  proportional 
to  ^r|p(^),  of  a  single  QW  structure  for  well  widths  of 
26  A  and  50  A,  respectively.  The  solid  curves  are  for  the  TE 
polarization  and  the  dashed  curves  for  the  TM  polarization.  A 
narrow  well  has  a  larger  separation  of  subbands,  resulting  in 
sharper  gain  and  spontaneous  emission  spectra  than  those  of 
a  wider  well.  We  plot  in  Fig.  8(a)  and  (b)  the  gain  spectrum 
for  tiiree  carrier  concentrations  ranging  from  1  x  10^^  cm*"^ 
to  3  X  10^^  cm'^  for  TE  (solid  curves)  and  TM  (dashed 
curves)  polarizations.  It  is  clearly  seen  that  the  TM  polarized 
gain  is  very  smaU,  as  expected  from  the  valence  subband 
structures,  since  the  HH  and  LH  subbands  contribute  mostly 
to  TE  polarization. 

The  spontaneous  emission  spectrum  is  useful  since  a  lot 
of  die  experimental  data  have  been  obtained  using  a  surface 
emission  configuration  via  optical  pumping  of  GaN  films 


Fig.  7.  Spontaneous  emission  specmim  (right  axis  in  arbitrary  units),  which 
is  proportional  to  (hi*;),  and  optical  gain  spectrum  (left  axis),  <7(hi*;),  for 
a  strained  GaN-Alo.3Gao.7N  quantum  well  with  a  well  width  of  (a)  =  26 

A  and  (b)  Lw  =  50  A  at  the  same  carrier  concentration  n  =  3  x  10^®  cm“^. 
Solid  curves  are  for  the  TE  polarization,  which  is  on  the  x-y  plane.  Dashed 
curves  are  for  the  TM  polarization. 

instead  of  an  edge-emission  configuration  due  to  the  difficulty 
of  cleaving  GaN  layers  to  form  mirrors.  The  spontaneous  emis¬ 
sion  speefrum  is  also  used  to  calculate  the  gain  spectrum  in  our 
model  based  on  (2)  and  (6).  Previously,  optical  gain  of  GaN 
zinc-blende  crystal  was  calculated  using  a  simple  parabolic 
band  model  [29]  and  using  a  many-body  formulation  [30] 
with  valence-band  mixing.  Our  results  for  wurtzite  structures 
have  the  same  order  of  magnitude  compared  with  the  previous 
calculations  [13],  [29],  [30]  for  both  crystal  structures  since  the 
effective  masses  are  about  the  same. 

The  peak  gain  coefficient  is  plotted  as  a  function  of  the 
carrier  concentration  per  quantum  well  in  Fig.  9(a).  The  trans¬ 
parent  carrier  densities  are  larger  than  1  or  1.5  x  10^^  cm“^ 
and  the  gain  increases  very  fast  once  the  carrier  densities  are 
above  the  transparent  values.  It  is  interesting  to  note  that  the 
earliest  experimental  report  [3]  of  stimulated  emission  from 
GaN  crystals  using  optical  pumping  has  a  gain  coefficient 
ranging  from  1.3  x  10^  to  1  x  10^  cm"“^  depending  on 
the  pump  intensity.  Although  the  transparent  carrier  denisty 
is  very  large,  ffie  gain  increases  quickly  with  increasing 
carrier  density  above  the  transparent-density  value.  The  gain 
coefficient  of  a  thinner  well  has  a  larger  transparent  carrier 
density,  yet  a  larger  differential  gain,  than  those  of  a  wider 
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Fig.  8.  Optical  gain  spectra  for  a  strained  GaN-Alo  .3Gao  .7N  quantum  well 
with  a  well  width  of  (a)  Lu-  =  26  A  and  (b)  Lu-  =  50  A  at  carrier 
concentrations  7?  =  1  x  10^®,  2x  10^^,  and  3  x  10^®  cm“^.  Solid  curves 
are  for  the  TE  polarization.  Dashed  curves  are  for  the  TM  polarization. 


well.  Both  the  carrier  density,  n,  and  the  material  gain  per 
well,  p,  are  proportional  to  or  the  density  of  states. 

The  modal  gain,  which  determines  the  threshold  condition  of 
a  laser  diode,  is  proportional  to  the  material  gain  multiplied 
by  the  optical  confinement  factor,  which  is  proportional  to  the 
well  width.  Since  the  injection  current  density  is  in  amperes 
per  unit  area,  we  replot  the  gain-weUwidth  product,  (p-Liy), 
versus  the  surface  carrier  concentration  in  Fig.  9(b)  for  another 
comparison.  The  surface  carrier  concentration,  Us  =  will 
be  the  number  of  carriers  injected  into  each  well  per  unit  area. 
The  plot  in  Fig.  9(b)  shows  an  interesting  improvement  in 
the  transparent  surface  carrier  density  and  differential  gain  by 
using  a  thin  well  rather  than  a  wide  well  in  the  active  region 
of  a  QW  structure. 

Using  the  electron  and  hole  effective  masses  [5],  [6], 
mc  *=  0.20mo  and  nth  =  O.SOmo,  the  dielectric  constant  [31], 
c  =  9.3860,  where  60  is  the  permittivity  of  free  space,  we 
obtain  an  exciton  binding  energy 


ntre^ 

2(47r€ft)2 


(42) 


equal  to  25  meV,  and  the  exciton  Bohr  radius 


as  = 


4%€fi^ 

VUrt^ 


(43) 


(b) 

Fig.  9.  (a)  Peak  optical  gain  coefficient  as  a  function  of  the  injected  carrier 
density,  n,  for  a  strained  GaN/Alo.3Gao.7N  quantum  well  with  a  well  width 
Lu  =  26  A,  and  i a-  =  50  A  and  0))  Gain-weilwidth  product,  which  is 
the  material  gain  g  in  (a)  multiplied  by  the  well  width  Lw  and  is  proportional 
to  the  modal  gain  coefficient  of  a  double  heterojunction  laser,  is  plotted  as  a 
function  of  the  injected  surface  carrier  density,  n*  =  nLv;,  The  polarization 
is  on  the  x-y  plane  (TE). 


about  31  A.  Here,  ntr  =  nt^nthlim^  +  m^)  is  the  reduced 
effective  mass.  We  then  estimate  the  carrier  density  above 
which  the  exciton  phase  becomes  unstable  with  respect  to  an 
electron-hole  plasma  for  a  two-dimensional  structure,  Uc  = 
l/(7ra|)  =  3.3x  10^^  cm"^^.  This  value  is  smaller  than  that 
of  ZnSe  [32],  5.7  x  10^^  cm“^.  From  Fig.  9(b),  we  find  a 
significant  amount  of  gain  occurs  when  die  surface  carrier 
concentration  n*  is  larger  than  4  x  10^^  cm“^.  Therefore, 
the  electron-hole  plasma  model  should  be  appropriate  in  the 
estimation  of  the  optical  gain.  It  should  be  noted  that  bandgap 
shrinkage  for  a  strained  wuitzite  QW  laser  and  an  enhancement 
of  the  matrix  elements  due  to  the  screened  Coulomb  interaction 
in  an  electron-hole  plasma  will  increase  the  gain  value  by 
about  30%  [30]  and  shift  the  peak  gain  energy  by  an  amount 
of  the  order  of  20-50  meV,  wWch  is  determined  by  the  carrier 
density.  On  the  other  hand,  strain  may  be  relaxed  in  the  GaN 
layer  because  of  dislocation,  and  the  amount  of  the  increase  in 
bandgap  energy  in  the  above  theory  due  to  compressive  strain 
may  be  ignored  in  strain-relaxed  samples. 

VI.  Conclusion 

We  have  presented  a  comprehensive  theoretical  model  for 
calculating  ^e  electronic  band  structures  and  optical  gain  for 
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a  wurtzite-strained  QW  structure.  Our  model  is  based  on  the 
effective-mass  Hamiltonian  derived  recently  using  the  k  •  p 
method  for  wurtzite  crystals.  We  also  presented  new  analytical 
expressions  and  numerical  results  for  the  QW  band  structure, 
the  basis  transformation,  the  interband  optical-matrix  elements, 
and  their  applications  to  the  calculations  of  the  optical  gain 
coefficient.  Our  model  provides  a  foundation  for  calculating 
the  optical  gain  and  spontaneous  emission  properties  of  a 
strained  wurtzite  GaN  QW  laser.  For  a  GaN  layer  with  much 
heavier  electron  and  hole  effective  masses  than  those  of  a 
zinc-blende  structure,  the  necessity  of  using  thin  multiple 
quantum  wells  is  pointed  out  and  demonstrated  in  our  model 
calculations. 
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